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An  experimental  and  analytical  investigation  of  saturated 
and  subcooled  flow  boiling  under  microgravity  conditions  was 
conducted  to  study  gravity  effects  on  two-phase  flow  boiling 
heat  transfer.  Optical  measurement  of  Freon  113  two-phase 
flow  patterns  in  a 0.5m  long,  5mmx5mm  vertical  channel 
uniformly  heated  with  transparent  conductive  films  identified 
bubbly,  slug,  and  annular  flow  under  microgravity  conditions. 
A significantly  different  void  distribution  was  observed  in 
bubbly  flow.  Velocity  and  temperature  distributions  in  the 
radial  direction  in  bubbly  flow  were  calculated  by  Sato  and 
Sadatomi's  method.  The  distributions  were  changed  by  the 
increased  void  fraction  and  coalescence  of  bubbles  under 
microgravity  conditions.  Analysis  of  experimental  results 
showed  that  the  bubble-induced  turbulence,  accounting  for  a 
large  fraction  of  heat  transfer  in  normal  gravity,  became 
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smaller  under  microgravity  conditions.  Effects  of  gravity  on 
bubble  agitation  were  less  pronounced  in  turbulent  flow  than 
in  laminar  flow  because  of  the  inherent  turbulence.  A 
substantial  increase  in  void  fraction  resulted  in  an  increase 
in  the  accelerational  and  frictional  pressure  drops  in  the 
low-quality  range,  suggesting  that  the  pressure  drop  in  the 
subcooled  region  is  also  significant  in  microgravity. 

Motivated  by  the  observed  changes  in  flow  and  heat 
transfer  rates  in  the  subcooled  region  under  microgravity 
conditions,  the  author  designed,  fabricated,  and  tested  an 
additional  forced  convective  boiling  rig  aboard  NASA's  low- 
gravity  Learjet.  A 0.5m  long,  11.3mm  diameter  vertical  copper 
tube  was  heated  uniformly  with  a constant  heat  flux  of 
23.9kW/m^.  Freon  113  in  the  once-through  flow  line  was  pumped 
by  air-driven  pistons  to  reduce  the  system  transient  time.  An 
increase  in  wall  superheats  observed  near  the  test  section 
exit  during  the  18  seconds  of  microgravity  implies  a decrease 
in  bulk  convection  in  the  absence  of  buoyancy  force.  The  heat 
transfer  analysis  showed  that  the  reduced  gravitational 
acceleration  affected  the  subcooled  two-phase  heat  transfer 
through  an  increase  in  void  fraction  and  stagnation  of 
coalesced  bubbles.  The  contribution  of  evaporation  heat 
transfer  was  small  in  the  low  heat-flux  range.  The  overall 
two-phase  flow  heat  transfer  coefficient  in  vertical  upward 
flow  was  decreased  under  microgravity  conditions. 
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CHAPTER  1 
INTRODUCTION 

1 . 1 Background 

Heat  transport  is  an  essential  issue  associated  with  all 
large-power  systems  used  in  space  where  environment  is 
characterized  by  low-gravity  and  vacuum.  The  heat  transport 
process  in  space  is  divided  into  two  categories;  a)  heat 
transport  from  spacecraft  to  environment,  and  b)  heat 
transport  inside  the  spacecraft.  For  heat  rejection  from  the 
spacecraft,  thermal  radiation  is  the  only  method  available  in 
space.  This  problem  has  been  intensively  discussed  since  the 
early  stage  of  space  development  and  several  new  versions  of 
radiator  systems  have  been  proposed[l] [2] . 

Heat  transport  and  fluid  behavior  inside  the  spacecraft 
are  practically  identical  to  ground-based  systems  when  single- 
phase flow  media  are  employed.  However,  as  heat  transport 
requirements  grow  to  the  30  kilowatt-meters  level  and 
beyond[3],  the  existing  single-phase  transport  systems  become 
inefficient  due  to  significant  temperature  drop  along  the  flow 
direction.  Temperature  drop  is  highly  undesirable  for  the 
radiator-dependent  systems  since  the  emissive  power  decreases 
proportional  to  the  4th  power  of  radiator  temperature. 
Although  passive  flow  recirculation  devices  such  as  heat-pipes 
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are  the  alternatives,  they  are  not  applicable  to  the  transport 
capacities  in  excess  of  50  kilowatt-meters[3] . The 
circulating  two-phase  flow  involving  phase  change  is  capable 
of  transporting  thermal  energy  with  small  temperature  drops. 
Because  of  unexpected  flow  behavior  in  the  absence  of  gravity, 
however,  the  use  of  two-phase  flow  has  been  intentionally 
avoided  regardless  of  its  favorable  features.  Analysis  of 
two-phase  flow  under  microgravity  conditions  is  formidable 
because  most  of  the  correlations  derived  in  normal-gravity 
environments  are  empirical  in  nature.  These  correlations 
cannot  be  applied  to  microgravity  conditions  without 
considering  the  gravity  effects  on  heat  transfer  mechanisms. 
The  gravity  term  appearing  in  the  ground-based  correlations 
needs  to  be  regarded  as  another  variable  under  microgravity 
conditions. 

The  experimental  efforts  to  investigate  the  effects  of 
gravity  on  boiling  heat  transfer  began  in  the  middle  1950 's 
when  the  need  for  scientific  and  practical  applications  of 
fluid  technology  in  space  became  increasingly  obvious. 
Siegel [4]  published  a comprehensive  review  of  low-gravity 
effects  on  heat  transfer  that  encompasses  the  efforts 
conducted  before  1966.  His  review  covers  the  four  fundamental 
areas  of  heat  transfer;  namely,  free  convection,  pool  boiling, 
forced-convection  boiling,  and  forced-convection  condensation. 

For  pool  boiling  in  saturated  liquids,  quite  a few 
experimental  data  was  collected  under  microgravity 
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conditions [5]  [6]  [7]  . Siegel [4]  reported  that  although  gravity 
had  little  effect  on  nucleate  boiling  on  wires  in  saturated 
liquid,  an  increase  in  heat  transfer  coefficient  was  observed 
for  flat  plate  heaters  upon  the  change  of  gravity.  Keshock 
and  Siegel [8] [9]  also  reported  that  the  bubble  growth  rates  in 
water  and  a 60%  by  weight  aqueous-sucrose  solution  at 
saturation  temperature  were  insensitive  to  gravity  change  in 
microgravity.  However,  Keshock' s experiments  were  conducted 
in  the  drop  tower  facility  providing  a 1.0  second  microgravity 
duration.  It  is  possible  that  the  microgravity  condition  was 
terminated  before  the  wall  temperature  became  steady. 

Compared  to  the  amount  of  pool  boiling  database,  little 
information  is  available  to  date  on  the  forced  convective 
boiling.  Papell[ll]  investigated  the  threshold  velocity 
beyond  which  the  critical  heat  flux  was  insensitive  to  the 
flow  direction.  Evans [12]  observed  bubble  distributions  over 
the  tube  cross  section  of  an  adiabatic  flow  under  microgravity 
conditions  and  found  that  they  were  similar  to  that  of 
vertical  flow  in  normal  gravity. 

More  recently,  a landmark  experiment  for  the  low-gravity 
two-phase  flow  was  conducted  by  Chen,  Downing,  Parish,  and 
Keshock[13].  They  observed  flow  patterns  by  a high  speed 
camera  and  measured  total  pressure  drops  in  adiabatic  straight 
and  curved  test  sections.  Only  slug  and  annular  flow  were 
observed.  The  measured  pressure  drops  in  a 6.7m,  15.8mm  ID 
horizontal  test  section  during  microgravity  operation  exceeded 
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those  measured  during  the  ground  testing  by  20%  at  the  same 
qualities.  This  change  is  caused  by  the  transition  of  flow 
patterns  from  stratified  to  slug  or  annular  flows. 

Determination  of  flow  regime  boundaries  is  another 
controversial  task.  Colin,  Fabre,  and  Dukler[14]  observed  the 
gas-liquid  two-phase  flow  pattern  and  attempted  to  find  the 
regime-transition  criteria.  The  results  were  mapped  on  a flow 
regime  diagram  as  a function  of  liquid  and  gas  superficial 
velocities.  The  predicted  criterion  agreed  with  the  boundary 
between  slug  and  annular  flow. 

More  recently,  Reinarts,  Best,  and  Hill [15]  measured  two- 
phase  condensation  heat  transfer  rates  of  Freon  12  in 
microgravity.  The  results  showed  that  the  condensation  heat 
transfer  was  significantly  reduced  during  the  microgravity 
period.  An  attempt  was  made  to  obtain  the  similar  form  of 
Nusselt  number  with  Dittus-Boelter  equation  assuming 
homogeneous  flow.  However,  gravity  effects  on  heat  transfer 
was  not  clearly  shown  in  their  modeling  and  the  correlation 
derived  was  highly  empirical . 

1.2  Scope  of  the  Study 

Although  numerous  efforts  have  been  devoted  to 
microgravity  pool  boiling  for  the  past  30  years,  the 
fundamental  mechanism  behind  the  observed  phenomena  is  yet 
unclear.  This  is  mainly  due  to  the  highly  empirical  nature  of 
the  existing  correlations  obtained  in  normal  gravity.  The 
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main  objective  of  the  present  work  was  to  identify  the 
influence  of  gravity  on  the  heat  transfer  mechanisms  in  flow 
boiling.  To  accomplish  this  goal,  first,  theoretical  bases  of 
existing  correlations  for  heat  transfer  and  pressure  drop  were 
examined.  Second,  the  following  specific  tasks  were  selected: 

- observation  of  flow  patterns  in  microgravity, 

- measurement  of  total  pressure  drops  of  flow  boiling 
in  a straight  tube, 

- measurement  of  void  fractions  in  saturated  and 
subcooled  boiling,  and 

- measurement  of  two-phase  heat  transfer  coefficients 
in  microgravity. 

Two  sets  of  microgravity  flow  boiling  loops  were  built 
and  tested  under  different  microgravity  environments.  The 
first  series  of  experiments  were  conducted  in  a 15m  drop  tower 
facility  providing  1.25  seconds  of  free  fall.  Flow  pattern 
observation  and  pressure  drop  measurement  were  the  main 
objectives  of  this  series  of  experiments. 

The  second  series  of  experiments  were  conducted  aboard 
the  NASA  Lear jet  allowing  approximately  18  seconds  of 
microgravity  conditions.  The  newly  designed  flow  loop  and 
data  acquisition  system  were  mounted  on  two  racks  matching 
with  the  requirements  of  flight  experiments.  To  maintain  the 
inlet  conditions  constant  and  reduce  the  transient  time,  a 
once-through  flow  line  driven  by  pressurized  pistons  was 
employed.  The  main  objective  of  this  experiment  was  to 
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measure  the  heat  transfer  coefficients  in  subcooled  flow 
boiling  under  microgravity  conditions.  The  investigation  on 
the  subcooled  region  was  particularly  motivated  by  the 
observation  of  significantly  different  flow  patterns  in  the 
subcooled  region  during  the  drop  tower  experiments.  An 
attempt  was  made  to  obtain  velocity  profiles  and  temperature 
distributions  in  the  subcooled  regions  under  microgravity 
conditions,  using  a method  proposed  by  Sato  and  Sadatomi[16] . 
Their  method  was  chosen  in  the  present  analysis  because  the 
contribution  of  bubble  agitation  to  be  affected  by  gravity 
change  was  differentiated  from  forced  convection  heat 
transfer.  Since  the  number  of  data  sets  obtained  from  the 
present  experiments  was  small,  emphasis  was  placed  on 
scrutinizing  particular  data  sets  from  various  aspects,  rather 
than  deriving  new  correlations  for  microgravity  conditions. 

Collected  data  on  flow  patterns,  pressure  drops,  and  two- 
phase  heat  transfer  coefficients  were  compared  with  respective 
data  obtained  in  normal  gravity  using  the  similar  test 
conditions.  Through  comparison,  the  influence  of  gravity 
change  on  forced  convection,  evaporation,  and  bubble  agitation 
was  investigated.  The  validity  of  several  existing 
mechanistic  correlations  for  two-phase  heat  transfer  was 
examined. 


CHAPTER  2 

PREVIOUS  EFFORTS  FOR  MICROGRAVITY  BOILING 

The  influence  of  gravity  on  pool  boiling  and  forced 
convection  boiling  are  reviewed  in  this  section.  Bubble 
growth  in  the  subcooled  and  saturated  liquid,  detachment  from 
the  wall,  and  nucleate  boiling  heat  transfer  for  pool  boiling 
are  essential  issues  of  relevance  to  forced  convective 
boiling.  However,  very  few  existing  correlations  have  ever 
been  applied  to  low-gravity  environments  due  to  lack  of 
understanding  of  two-phase  flow  behaviors  under  variable 
gravities. 

The  initial  approach  is  thus  to  investigate  the  influence 
of  gravity  on  the  parameters  governing  flow  and  heat  transfer. 
First,  previous  studies  on  the  influence  of  gravity  in  pool 
boiling  are  reviewed.  This  section  also  includes  a discussion 
of  the  role  of  the  buoyancy  force  during  the  bubble 
detachment.  Second,  several  experimental  results  for  forced 
convective  boiling  under  microgravity  conditions  are  reviewed. 
The  major  achievements  and  findings  for  both  pool  boiling  and 
forced  convective  boiling  are  listed  at  the  end  of  this 
chapter.  Finally,  the  motivation  of  the  present  study  and 
unique  features  of  the  present  equipment  are  introduced. 
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2.1  Pool  Boiling  in  Microaravity 

Investigation  of  the  gravity  effects  on  pool  boiling  is 
of  critical  importance  to  understand  the  bubble  growth  through 
the  superheated  layer  in  flow  boiling. 

Reduced  gravity  experiments  performed  in  the  early  1960s 
are  noteworthy.  Sherley[5]  obtained  boiling  curves  for  LH2 
under  both  normal  and  microgravity  conditions  using  a drop 
tower  and  a KC-135  aircraft.  Little  difference  was  observed 
between  the  normal  and  microgravity  boiling  curves.  Sherley 
concluded  that  the  effect  of  gravity  on  the  nucleate  boiling 
was  negligible.  However,  it  should  be  noted  that  results  of 
Sherley 's  work  can  only  draw  qualitative  conclusions  since  the 
actual  gravity  level  and  other  flow  conditions  during  the 
microgravity  experiments  are  not  shown  in  the  literature. 

Siegel  and  Keshock[8]  detected  changes  in  superheats  of 
vertical  and  horizontal  wires  in  distilled  water,  ethyl 
alcohol,  and  60%  by  weight  aqueous-sucrose  solution.  As  listed 
in  Table  2-1,  a decrease  in  superheat  by  2-3 °C  in  microgravity 
was  observed  for  the  horizontal  wires  in  water  and  ethyl 
alcohol,  while  an  increase  in  superheat  was  observed  for 
vertical  wires.  No  effect  of  gravity  on  the  superheat  was 
observed  for  sucrose  solution.  This  fact  suggests  that  bubble 
growth  and  detachment  controls  nucleate  boiling  heat  transfer 
without  the  influence  of  surface  tension.  Siegel  et  al.  also 
analyzed  the  effect  of  gravity  on  the  contact  angle  through 
visual  observation  filmed  by  a high  speed  camera  and  found 


Table  2-1  Heat  transfer  data  for  nucleate  boiling  from  an  electrically  heated  wire[8]. 
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that  the  gravity  change  affects  the  contact  angle  only 
slightly  for  the  above  fluids. 

Clodfelter[7]  investigated  the  effect  of  gravity  on  the 
critical  heat  flux  (CHF)  for  0.5mm  platinum  horizontal  wires 
and  6.4mm  wide  ribbons  in  distilled  water  using  a 16.8m  drop 
tower  facility  and  an  airplane.  The  CHF  measured  at  gravity 
fields  less  than  O.Olg  was  about  O.lStimes  the  normal  gravity 
value.  The  comparison  with  the  theoretical  g^''^  dependence 
showed  that  the  existing  CHF  correlations  qualitatively  agree 
with  the  measured  CHF  in  microgravity  as  well.  However, 
little  information  is  available  on  the  CHF  when  zero  gravity 
is  approached.  Since  the  fundamental  wave  length  on  the 
liquid-vapor  interface  is  a function  of  gravity,  the  existing 
CHF  correlations  should  be  reassessed  in  the  low  gravity 
range . 

A comprehensive  boiling  curve  for  microgravity  pool 
boiling  was  first  obtained  by  Merte  and  Clark[6].  Figure  2-1 
shows  the  boiling  curves  for  liquid  nitrogen  obtained  under 
normal  and  microgravity  conditions.  In  the  nucleate  boiling 
region,  the  gradient  of  the  curve,  namely,  the  heat  transfer 
coefficient,  in  microgravity  was  almost  the  same  as  that  under 
normal  gravity.  However,  the  corresponding  heat  fluxes  at  the 
same  superheat  had  slightly  decreased.  The  CHF  in 
microgravity  decreased  by  a factor  of  2.  In  the  film  boiling 
region,  both  the  heat  flux  and  heat  transfer  coefficient 
decreased  probably  due  to  the  increased  irregularity  of  the 
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Figure  2-1  Boiling  curves  obtained  by  Merte  and 

Clark  under  microgravity  conditions [ 6 ] . 
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interface.  These  results  from  Merte's  pool  boiling  experiment 
suggest  that  the  overall  performance  of  heat  exchangers  to  be 
used  under  microgravity  conditions  are  poorer  than  in  normal 
gravity  when  operated  near  the  critical  heat  flux  region. 

The  influence  of  gravity  on  pool  boiling  heat  transfer 
was  examined  using  widely  accepted  correlations  in  normal 
gravity.  Pool  boiling  is  categorized  into  four  regions  in 
terms  of  the  dominant  heat  transfer  mechanism.  The  following 
theoretical  and  semi-empirical  correlations  for  each  region 
were  particularly  selected  because  the  gravity  term  was 
explicitly  included  the  following: 

1)  Natural  Convection  (Re  > 2000)  (Jakob  [17]) 
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2)  Nucleate  Boiling  (q”=constant)  (Rohsenow  [18]) 
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3)  Critical  Heat  Flux  (Zuber  [10]) 
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4)  Film  Boiling  (Frederking/Clark  [19]) 
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Although  all  of  the  above  correlations  include  the 
gravity  term,  the  physical  background  for  each  equation  is 
different.  In  natural  convection,  the  difference  in  specific 
weight  at  different  temperatures  causes  a convective  motion 
within  the  liquid.  In  nucleate  boiling,  the  buoyancy  force 
acting  on  the  vapor  bubble  is  one  of  the  dominant  forces  upon 
bubble  detachment.  Both  the  critical  heat  flux  and  film 
boiling  heat  flux  are  ruled  by  the  Kelvin-Helmholtz 
instability  of  the  liquid-vapor  interface.  However,  the  above 
correlations  derived  under  normal  gravity  conditions  assume 
the  gravitational  acceleration  constant.  All  heat  fluxes 
predicted  by  the  above  correlations  approach  zero  when 
gravity  becomes  extremely  small.  Under  such  conditions,  the 
soundness  of  these  equations  is  doubtful. 

The  above  correlations  were  all  derived  through 
dimensional  analysis.  The  exponents  attached  to  the 
dimensionless  numbers  were  empirically  obtained  through  curve 
fitting  premised  on  a constant  value  for  the  gravitational 
constant,  g=9.8m/s.  Under  variable  conditions,  however,  the 
terms  or  dimensionless  numbers  related  to  the  gravity  term  are 
expected  to  vary  significantly  from  those  in  normal  gravity. 
In  addition,  as  a result  of  the  diminished  gravity  force,  the 
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surface  tension  and  Marangoni  forces  become  dominant  and 
control  momentum  and  heat  transfer.  It  is  likely  that  those 
values  such  as  the  Grashof  number  or  the  Laplace  constant  are 
overestimated  under  microgravity  conditions.  Thus,  the 
exponents  obtained  under  normal  gravity  conditions  should  be 
re-evaluated  under  microgravity  conditions  considering  the 
effect  of  variable  gravity. 

Table  2-2  summarizes  the  gravity  effects  predicted  from 
the  above  correlations.  The  ratio  of  the  pool  boiling  heat 
flux  in  normal  gravity  to  that  in  microgravity  was  calculated 
from  Merte's  experiments  conducted  in  the  drop  tower.  It  is 
noted  that  the  CHF  data  obtained  in  his  experiments  roughly 
follow  as  predicted  in  Eg.  (2-3)  . Latent  heat  transport 
and  bulk  turbulent  process  are  the  two  major  mechanisms  of 
boiling  heat  transfer.  The  dominant  component  has  yet  to  be 
determined  even  under  normal  gravity  conditions.  The 
microgravity  environment,  free  from  buoyancy-driven 
convection,  offers  an  ideal  stage  to  study  the  problem 
proposed  above. 

Schwartz  and  Adelberg[20]  investigated  the  significance 
of  latent  heat  transport  of  saturated  water  at  1 atm  under 
microgravity  conditions  using  an  airplane  providing  a low 
gravity  duration  of  8 - 10  seconds.  Figure  2-2  shows  that  the 
ratios  of  heat  transport  due  to  evaporation  to  the  total  heat 
flux  increased  as  the  acceleration  level  was  reduced.  The 
values  of  latent  heat  were  calculated  from  the  volume  of  vapor 


Table  2-2  Effects  of  gravity  on  heat  transfer  mechanism  in  pool  boiling. 
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Figure  2-2  Pool  boiling  curves  for  water  under  normal 

and  microgravity  conditions,  (a)  Total  heat 
flux,  and  (b)  Latent  heat  transport [20] . 
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bubbles  in  photographs.  The  total  heat  flux  remained 
unchanged  in  normal  and  microgravity.  This  is  presumably  due 
to  the  suppression  of  natural  convection  caused  by  the  loss  of 
buoyancy  force.  It  is  considered  that  under  microgravity 
conditions,  the  heat  transport  due  to  natural  convection  was 
replaced  by  the  enhanced  evaporation  to  maintain  the  same 
total  heat  flux. 

More  recently,  Zell,  Straub,  and  Weinzierl [21]  obtained 
a boiling  curve  near  the  onset  of  nucleate  boiling  for  a gold 
film  in  saturated  and  subcooled  Freon  113  under  microgravity 
conditions.  The  experiment  was  conducted  in  a sounding  rocket 
providing  about  6 minutes  of  microgravity  conditions  on  the 
order  of  10‘®g,  which  is  an  extremely  long  and  small 
acceleration  level  compared  with  previous  microgravity 
environments.  Figure  2-3  displays  the  comparison  of  boiling 
curves  in  normal  and  microgravity  conditions  near  the  onset  of 
nucleate  boiling.  The  boiling  curves  indicated  that  the 
superheat  to  initiate  the  nucleation  in  microgravity  becomes 
larger  than  in  normal  gravity,  and  that  the  slope  of  the 
boiling  curve  after  the  onset  of  nucleate  boiling  is 
significantly  smaller  than  in  normal  gravity.  A prominent 
difference  in  bubble  growth  on  the  heated  surface  in  saturated 
Freon  113  under  microgravity  conditions  is  shown  in  Figure  2- 
4 . The  growing  bubbles  covered  the  total  heating  area  until 
they  reached  the  boundary  of  a container,  and  then  broke  into 
pieces.  It  appeared  that  the  buoyancy  force  was  not  strong 
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enough  to  detach  the  bubble  in  saturated  Freon.  This 
experiment  suggests  that  a steady  state  for  saturated 
pool  boiling  may  not  be  reached  since  bubbles  do  not  detach 
from  the  wall  and  keep  growing.  This  observation  of  pool 
boiling  phenomena  indicate  that  the  heat  transfer  is  expected 
to  deteriorate  under  microgravity  conditions. 

Ervin,  Merte,  Keller,  and  Kirk[22]  conducted  transient 
nucleate  boiling  experiments  in  a drop  tower  facility  and 
observed  boiling  phenomena  of  subcooled  Freon  113  on  a 400A 
thick  gold  film  sputtered  on  a quartz  substrate.  The  boiling 
spread  across  the  heater  surface  at  the  moment  of  boiling 
inception  and  the  propagation  process  was  classified  into  six 
distinct  categories.  Ervin  et  al.  reported  that  the  energetic 
type  of  boiling  propagation  on  the  heating  surface  occurred 
uniquely  in  microgravity.  It  is  assumed  that  the  energetic 
boiling  propagation  was  due  to  build-up  of  enthalpy  in  the 
superheated  layer. 

Bubble  growth  in  pool  boiling  is  mainly  controlled  by  the 
surface  condition  of  the  heater  and  the  thermal  properties  of 
the  fluid.  Zuber[23]  derived  a correlation  for  the  growth  of 
a single  bubble  in  homogeneous  saturated  liquid 


D{t) 


f 1/2 

2 hfgPg{TlCC)^/^ 


2q» 

hfgPg 


t 


(2-5) 


The  rate  of  growth  predicted  by  the  above  equation  is 
proportional  to  t^^^  and  agreed  with  the  experimental  results. 
Eg. (2-5)  also  indicates  that  the  growth  rate  is  independent  of 
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gravity  unless  other  conditions  are  changed.  However, 
according  to  the  results  reported  by  Merte  et  al. [6]  and  Ervin 
et  al.[22],  the  wall  superheats  changed  for  the  same  heat 
fluxes  in  microgravity  as  the  boiling  progressed.  Zell  et 
al.[21]  concluded  that  gravity  dependent  mechanisms  such  as 
bubble  agitation  and  natural  convection  were  completely 
compensated  by  the  phase  change  mechanism  under  microgravity 
conditions,  which  may  influence  the  growth  rate  of  bubbles. 

The  bubble  detachment  diameter  in  pool  boiling  is  another 
non-consensus  issue  under  both  normal  and  microgravity 
conditions.  Fritz [24]  considered  the  balance  of  the  buoyancy 
force  and  the  surface  tension  force  at  the  moment  of 
detachment.  The  model  can  be  applied  to  the  bubble  growth 
with  small  growth  rates  such  as  gas  bubbles.  However,  there 
exist  dynamic  forces  for  a rapidly  growing  bubble  that  affects 
the  bubble  detachment.  An  analysis  of  the  forces  acting  on 
the  bubble  in  normal  and  microgravity  was  conducted  by  Keshock 
and  Siegel[9].  Inertia  of  liquid  initially  assists  the  bubble 
detachment.  However,  it  opposes  the  detachment  after  the 
growth  rate  is  decelerated.  They  concluded  that  the  bubble 
detachment  in  saturated  water  did  not  occur  in  the  absence  of 
gravity  because  the  surface  tension  force  exceeded  the  sum  of 
the  buoyancy  force  and  inertia  force. 

Cochran,  Adylott,  Frysinger[62]  reported  that  the  effect 
of  gravity  became  small  as  subcooling  increased.  A horizontal 
chromel  strip  was  heated  in  distilled  water  at  1 atm  with 
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variable  subcooling  and  subjected  to  2.3  seconds  of 
microgravity  free  fall.  The  vapor  volume  of  a single  bubble 
was  measured  as  a function  of  time  from  the  beginning  of 
growth  until  the  end  of  detachment.  For  subcooling  more  than 
16.7°C,  the  effect  of  gravity  on  the  vapor  volume  was 
negligible.  It  is  inferred  that  this  is  due  to  the  increased 
recondensation  of  vapor  bubbles  and  that  the  detachment 
diameter  is  determined  by  the  thickness  of  the  thermal 
boundary  layer.  The  effect  of  gravity  on  the  number  of 
nucleation  sites  was  also  investigated.  The  visual 

observations  detected  little  influence  of  gravity  on  the 
number  of  nucleation  sites.  However,  the  above  conclusions 
conflict  with  the  results  from  Siegel  and  Keshock's 
experiments [8 ] , reporting  a decrease  in  heater  temperature. 
To  investigate  the  gravity  effect  on  nucleation,  close 
observation  of  activated  nucleation  sites  on  the  heating 
surface  with  visual  aids  is  necessary. 

Recently,  Klausner  and  Zeng[26]  proposed  an  unified 
bubble  detachment  model  to  predict  vapor  bubble  departure 
diameters  in  both  pool  and  flow  boiling  conditions.  A new 
expression  to  describe  the  drag  force  acting  on  the  bubble  was 
obtained  through  numerical  techniques.  The  force  balances 
horizontal  and  perpendicular  to  the  heating  surface  were 
considered  separately  to  explain  the  mechanism  of  sliding 
bubbles  before  they  lift  off  the  wall.  Klausner  et  al. 
concluded  that  the  magnitude  of  the  surface  tension  force 
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accounts  for  only  a fraction  of  other  forces  such  as  unsteady 
drag  and  is  not  significant  at  the  points  of  departure.  For 
the  lift-off  diameter  of  a bubble,  a simple  force  balance 
between  the  buoyancy  and  the  inertia  forces  may  be  used.  They 
also  pointed  out  that  the  departure  process  is  strongly 
dependent  on  the  pressure  through  the  bubble  growth  rate. 

The  following  is  a list  of  major  results  from  pool 
boiling  experiments  under  microgravity  conditions: 

a)  Significance  of  the  surface  tension  in  nucleate  boiling 

was  examined  using  different  heater  geometries  and 

fluids.  When  the  surface  tension  is  dominant  upon  bubble 
detachment,  gravity  effects  on  the  nucleate  boiling  was 
insignificant [8] . 

b)  The  slope  of  the  boiling  curve  obtained  under 

microgravity  in  the  nucleate  boiling  region  remained 

almost  the  same,  indicating  that  the  heat  transfer 
coefficient  was  unchanged.  However,  CHF  and  film  boiling 
roughly  followed  g^^'*  and  g^^^,  respectively,  as  predicted 
by  the  existing  correlations [ 6] . 

c)  The  heat  transfer  coefficient  near  the  ONE  was 

significantly  reduced.  The  wall  superheat  that  initiate 
nucleate  boiling  was  increased [21] . 

e)  The  bubble  detachment  diameter  in  saturated  liquid  became 
significantly  large  under  microgravity [21] . 

d)  For  constant  heat  fluxes,  the  fraction  of  latent  heat 
transfer  tends  to  increase  as  the  acceleration  level 
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decreases [ 20 ] . Bulk  turbulent  convection  appeared  to  be 
suppressed  due  to  the  lack  of  buoyancy-driven  flow. 

f)  In  subcooled  pool  boiling,  the  effect  of  gravity  was 
negligible  for  subcoolings  of  over  16. 7K.  The  number  of 
nucleation  sites  was  unchanged  under  microgravity [8] . 

g)  The  contact  angle  appeared  unchanged  under  microgravity 
conditions [9] . 

2.2  Forced  Convective  Boiling  in  Microqravitv 

Compared  to  abundant  research  work  for  pool  boiling,  only 
a limited  amount  of  study  has  been  performed  on  forced 
convective  boiling  in  microgravity.  This  is  because 
microgravity  environments  are  difficult  to  simulate  and  quite 
often  costly  to  create  in  laboratories.  Although  the 
characteristic  period  of  bubble  formation  is  on  the  order  of 
10'^  seconds,  at  least  several  minutes  are  required  for 
conventional  forced  convective  boiling  loops  to  reach  steady 
states.  The  drop  tower  facilities  provide  2-6  seconds  of 
low-gravity  environment,  while  the  microgravity  aircraft  such 
as  Lear jet  or  KC-135  provide  about  20  seconds.  Regardless  of 
difficulties  described  above,  several  compact  flow  boiling 
loops  have  been  built  and  tested  aboard  the  reduced-gravity 
airplanes.  These  experimental  results  are  reviewed  in  the 
following  paragraph  and  major  contributions  are  summarized  at 


the  end. 
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For  a better  understanding  of  the  changes  in  two-phase 
flow  under  microgravity , visual  observation  of  the  flow 
pattern  change  is  necessary.  Evans [12]  observed  the  bubble 
distribution  of  an  air-water  mixture  in  a 12.7mm  diameter, 
45.7cm  long  adiabatic  plastic  tube  under  microgravity 
conditions  and  found  that  the  flow  patterns  were  similar  to 
that  of  vertical  upward  flow  under  normal  gravity.  Bubbly, 
slug,  and  annular  flow  patterns  were  the  three  major  flow 
regimes  classifying  the  void  distributions  for  two-phase  flow 
under  normal  and  microgravity  conditions.  Current  issues 
under  this  topic  are  to  investigate  the  phase  transition 
mechanism  and  thus  to  determine  the  boundaries  of  these  flow 
patterns.  Although  several  researchers [13] [14]  proposed  flow 
regime  maps  applicable  to  microgravity  conditions,  conclusive 
results  have  not  yet  been  obtained. 

Information  on  the  pressure  drop  is  essential  for  the 
design  of  heat  exchangers  in  space.  Unfortunately,  only  a 
limited  amount  of  research  has  been  done  in  this  area. 

Chen,  Downing,  Parish,  and  Keshock[13]  measured  the 
pressure  drop  in  15.8mm  diameter,  1.83m  long  straight  and 
curved  tubes  using  Freon  114.  A mixture  of  vapor  and  liquid 
with  the  quality  ranging  from  0.05  to  0.9  was  produced  in  a 
swirl  flow  evaporator  and  introduced  into  an  adiabatic 
transparent  test  section  allowing  for  the  observation  of  flow 
patterns.  The  measured  total  pressure  drops  in  microgravity 
were  larger  than  those  obtained  in  the  horizontal  test  section 
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in  normal  gravity  for  the  same  qualities.  However,  the  direct 
comparison  of  pressure  drops  is  questionable  because  the  flow 
patterns  are  different.  Besides,  in  adiabatic  systems,  the 
bubble  formation  on  the  wall  in  the  nucleate  boiling  region  is 
neglected  and  the  frictional  component  of  pressure  drop  is 
underestimated . 

Adiabatic  dispersed  flow  experiments  using  water  and  n- 
butyl  benzoate  of  identical  densities  were  carried  out  under 
normal  gravity  conditions  by  Sridhar,  Chao,  and  Soo[27]. 
Pressure  drops  along  a 2.44m  long,  45mm  square  horizontal  test 
section  were  measured  by  changing  the  ratio  of  n-butyl 
benzoate  in  the  mixture  to  simulate  the  variation  of  void 
fraction.  It  was  observed  that  the  frictional  pressure  drop 
was  unaffected  by  large  changes  in  the  volume  fraction  of 
benzoate  drops  and  their  size  distribution.  However,  the 
adiabatic  systems  fail  to  simulate  the  void  fraction 
distribution  near  the  wall.  In  actual  heat  exchangers,  the 
bubble  formation  on  the  heated  wall  needs  to  be  considered  in 
the  nucleate  boiling  region  since  the  momentum  of  flow  is  lost 
substantially  near  the  wall. 

Klausner[25]  studied  the  influence  of  gravity  on  two- 
phase  flow  by  changing  the  orientation  of  flow.  Heat 
transfer,  pressure  drop,  and  void  fraction  of  Freon  11  in  a 
19.1mm  diameter,  1.48m  long  copper  tube  was  measured  in 
vertical  up-flow,  vertical  down-flow,  and  horizontal  flow 
configurations.  The  observed  flow  pattern  was  mainly  in  the 
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annular  flow  regime.  Gravity  exerts  its  influence  on  vapor 
volume  fraction  through  changes  in  the  relative  velocity 
between  the  liquid  and  vapor  phase.  It  is  also  found  that 
both  the  ebullition  process  and  bulk  turbulent  convection  of 
flow  boiling  are  affected  by  gravity  through  the  thickness  of 
the  liquid  film  in  the  annular  flow  range. 

Under  normal  gravity  conditions,  heat  transfer  in 
vertical  and  horizontal  tubes  are  different.  Under 
microgravity  conditions,  however,  the  absence  of  buoyancy 
force  cancels  the  difference  due  to  flow  orientation  with 
respect  to  the  gravity  field.  Papell,  Simoneau,  and  Brown [11] 
investigated  the  effect  of  flow  orientation  on  the  CHF  for 
liquid  nitrogen  flowing  through  a vertical  heated  tube  1.28cm 
in  diameter  at  0.516MPa  with  a subcooling  of  10. 6K  and  found 
that  the  CHF  was  independent  of  the  flow  orientation  beyond 
flow  velocities  greater  than  1 . 95m/sec (Figure  2-5).  This 
velocity  corresponds  to  a Froude  number  of  31.2,  indicating 
that  the  inertia  force  is  about  30  times  larger  than  the 
gravity  force. 

Papell  [28]  also  observed  the  forced  convective  boiling  of 
distilled  water  in  a 7.9mm  diameter  stainless  tube  in 
microgravity.  A short  flow  interruption  at  the  beginning  of 
zero  gravity  conditions  was  reported,  followed  by  an  increase 
in  the  heat  transfer  coefficient  on  the  order  of  16%.  This 
observation  agrees  with  Siegel's  results  from  the  pool  boiling 
experiments,  in  which  the  heat  transfer  coefficients  had 
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Figure  2-5  Critical  heat  flux  as  a function  of  inlet 

velocity  for  upward  and  downward  flow  of 
liquid  nitrogen  in  a 0.505"  diameter  tube 
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increased  under  microgravity  conditions.  For  pool  boiling,  as 
the  bubble  size  becomes  larger,  the  nuinber  of  bubbles  emitted 
from  the  wall  becomes  smaller  when  a constant  heat  flux  is 
applied.  Zell  et  al.  observed  such  formation  of  large  bubbles 
in  their  microgravity  experiments.  If  evaporation  is  a 
dominant  heat  transfer  mechanism,  the  heat  transfer  rate  of 
flow  boiling  systems  shows  similar  tendency  to  that  of  pool 
boiling  systems.  However,  the  contributions  of  convective  and 
evaporation  heat  transfer  is  still  unclear  under  microgravity 
conditions.  Precise  measurement  of  single-phase  heat  transfer 
rate,  void  fraction,  and  wall  temperatures  along  the  flow  axis 
is  needed  to  resolve  this  issue.  Unfortunately,  these  data 
were  not  available  in  Papell's  experiment. 

Feldmanis[29]  conducted  a flight  experiment  aboard  a KC- 
135  microgravity  aircraft  and  measured  temperatures  and 
pressures  in  a closed,  forced  convective  boiling  and 
condensing  loop.  An  increase  in  the  system  pressure  and  an 
attenuation  of  pressure  oscillations  in  the  9.53mm  diameter, 
0.91mm  long  evaporator  was  observed  during  the  period  of 
microgravity.  No  noticeable  temperature  change  in  the  bulk 
temperature  was  observed  at  the  outlet.  Quantitative  analysis 
had  not  been  made  since  the  data  on  the  flow  rate  and  thermal 
equilibrium  quality  were  not  measured.  The  pressure 
oscillations  are  probably  caused  by  the  formation  of  slug 
bubbles  which  are  usually  observed  in  the  beginning  of 
saturated  boiling.  Feldmanis  agreed  on  the  increase  in  heat 
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transfer  rate  in  microgravity , the  fact  first  pointed  out  by 
Papell,  and  explained  such  an  increase  in  terms  of  the 
transition  of  flow  patterns  from  the  bubbly  to  slug  flow.  As 
a pioneering  work  in  microgravity  flow  boiling,  his  experiment 
introduced  a number  of  questions  as  well.  However,  the 
discussion  was  speculative  due  to  lack  of  experimental 
evidence. 

Crowley  and  Sam [30]  observed  stable  operation  of  the  two- 
phase  test  rig  aboard  the  NASA  KC-135  aircraft  under 
microgravity  conditions.  However,  the  wall  temperatures  on  a 
condenser  during  the  microgravity  trajectories  appear  to  be  in 
transients.  The  system  pressure  increased  about  5%  (ISkPa)  in 
microgravity  and  correspondingly  saturated  fluid  temperatures 
in  the  system  increased  by  IK  to  2K.  Crowley  et  al.  reported 
that  the  heat  transfer  rates  in  the  evaporator  remained 
constant  in  microgravity.  This  result,  however,  contradicts 
the  findings  of  Feldmanis  and  Papell.  Since  wall  temperature 
profiles  of  the  evaporator  are  not  shown  in  the  report,  it  is 
difficult  to  confirm  their  statement.  The  thermal  equilibrium 
quality,  flow  patterns,  and  void  fraction  should  be  measured 
simultaneously  to  evaluate  the  heat  transfer  rate  in 
microgravity . 

Standley  and  Fairchild [31]  examined  an  operation  of  a R- 
11  two-phase  pumped  loop  system  aboard  the  NASA  KC-135 
microgravity  aircraft  as  a preliminary  run  for  subsequent 
suborbital  rocket  tests.  Mass  fluxes  were  lOOkg/m^/s  and 
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144kg/m^/s  and  heater  powers  were  280W  and  420W.  The  use  of 
Rll  allowed  thermal  equilibrium  qualities  large  enough  to 
exceed  85%  of  the  void  fraction  at  the  boiler  exit.  The 
magnitude  of  temperature  increases  in  the  boiler  fluid  and  on 
the  boiler  wall  were  both  1 to  2K,  indicating  a slight 
decrease  in  heat  transfer  rate  under  microgravity  conditions. 
This  fact  opposes  the  findings  of  Papell,  Feldmanis,  and 
Crowley  et  al.  Note  that  the  flow  pattern  at  these  high 
qualities  is  annular,  in  which  thin  liquid  film  remaining  on 
the  heating  wall  surrounds  vapor  flow  at  the  center.  Standley 
et  al . postulated  that  the  above  temperature  increase  was  due 
to  change  in  two-phase  heat  transfer  mechanism.  Due  to  lack 
of  mechanistic  models  and  limited  database  for  microgravity 
flow  boiling,  their  hypothesis  has  not  been  verified  yet. 

The  following  is  a list  of  major  results  and  findings 

made  by  previous  investigators  for  flow  boiling  under 

microgravity  conditions: 

a)  Beyond  a certain  value  of  Froude  number,  the  effect  of 
gravity  on  CHF  can  be  negligible  since  the  inertia  force 
becomes  dominant  over  the  gravity  force [11]. 

b)  The  three  similar  flow  patterns  were  observed  under 

microgravity  although  the  structure  of  each  flow  pattern 
was  different [ 12 ] . 

c)  Several  flow  regime  maps  have  been  proposed  to  determine 

the  boundaries  of  each  flow  regime  for  microgravity 
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use[13].  However,  none  of  them  successfully  classified 
the  flow  patterns  observed  in  microgravity. 

d)  The  total  pressure  drop  in  the  horizontal  straight  and 
curved  test  section  for  two-phase  flow  increased  probably 
due  to  flow  pattern  change [13]. 

e)  The  heat  transfer  coefficient  under  microgravity  was 
increased  by  16% [11].  However,  other  investigators 
contended  otherwise [ 30] . 

The  review  of  previous  work  has  revealed  a substantial 
lack  of  database  on  flow  boiling  in  microgravity.  Especially, 
the  change  of  two-phase  flow  heat  transfer  rate  in 
microgravity  is  unclear.  Regardless  of  strenuous  experimental 
work  done  in  the  past,  no  mechanistic  models  capable  of 
accounting  for  gravity  change  has  been  proposed  yet.  This  is 
because  the  drop  towers  or  microgravity  aircraft  cannot 
provide  enough  period  of  microgravity  for  each  component  of 
flow  boiling  loops.  The  effects  of  gravity  on  boiling, 
condensing,  and  adiabatic  two-phase  mixtures  are  still 
unidentified.  All  the  previous  experimental  facilities 
reviewed  here  consist  of  evaporators,  condensers,  pumps,  and 
flow  lines  operated  simultaneously  under  unknown  influence  of 
gravity.  As  Feldmanis  pointed  out,  19  seconds  or  less  is  too 
short  for  reaching  equilibrium  for  the  entire  flow  loop 
because  of  long  transit  time  in  the  flow  loop.  Within  a 
limited  amount  of  time  allowed  for  drop  tower  and  flight 
experiments,  focused  items  need  to  reach  a steady  state  to 
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measure  temperatures,  pressures,  and  flow  rates  under  thermal 
equilibrium  conditions. 

The  present  study  focuses  exclusively  on  convective  flow 
boiling  under  microgravity  conditions.  To  eliminate  the 
feedback  effect  on  the  boiling  section,  pumps  were  replaced 
with  a once-through  flow  line  and  two  separate  reservoirs. 
The  pistons  equipped  with  the  reservoirs  are  driven  by 
pressurized  air  provided  from  a high  pressure  tank  through  a 
pressure  regulator.  Under  the  test  conditions  considered 
here,  the  transit  time  in  a 0.5m  boiling  test  section  is 
approximately  2 to  5 seconds.  Test  conditions  were 
deliberately  confined  to  low  flow  rates  and  low  heat  fluxes, 
since  the  influence  of  gravity  may  be  less  significant  in  high 
flow  rates  where  the  inertia  force  is  overwhelmingly  larger 
than  the  buoyancy  force. 


CHAPTER  3 

EXPERIMENTAL  DESIGN  AND  PROCEDURE 
3 . 1 Drop  Tower  Experiment 

A simple  flow  boiling  loop  using  Freon  113  was  built  for 
the  drop  tower  experiment.  Schematic  diagram  of  the  boiling 
loop  is  shown  in  Figure  3-1.  The  loop  mainly  consists  of  the 
boiling  test  section,  vapor-liquid  separator,  cooler, 
positive-displacement  gear  pump,  and  pre-heater.  The  vapor 
generated  in  the  test  section  was  released  into  the  atmosphere 
and  only  liquid  was  recovered  and  recirculated.  The 
preheater,  consisting  of  a coiled  copper  tube  wrapped  with  an 
adjustable  lOOW  heater,  allowed  fine  adjustment  of  subcooling 
at  the  inlet  of  the  test  section  by  controlling  the  applied 
voltage  to  the  heater.  The  liquid  recovered  in  the  separator 
was  cooled  while  passing  through  the  cooler  immersed  in  iced 
water.  The  selection  of  a positive-displacement  gear  pump, 
capable  of  providing  a constant  liquid  flow  rate  regardless  of 
the  pressure  in  the  downstream,  was  made  to  deal  with  the 
unexpected  pressure  change  in  the  boiling  test  section.  Each 
component  was  connected  with  4mm  I.D.  Teflon  and  copper  tubes. 
Flow  rate  was  controlled  by  throttling  the  flow  with  a valve 
on  the  line  bypassing  the  pump. 
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Figure  3-1  Schematic  diagram  of  a drop  tower  flow 

boiling  loop. 
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Mass  flow  rate  was  measured  in  the  upstream  of  the 
boiling  test  section  with  a turbine  flowmeter  which  is 
insensitive  to  flow  orientation  and  gravity  change.  The 
measurable  range  of  the  flow  meter  was  0 - 300cc/min.  Since 
the  loop  is  open  to  ambient  at  the  separator,  pressure  at  the 
exit  was  assumed  latm.  Bulk  temperatures  at  the  inlet  and 
outlet  temperatures  were  measured  by  two  C-C  type 

thermocouples  immersed  in  the  flow.  The  junctions  of  the 
thermocouples  were  placed  at  the  center  of  the  flow. 

The  volume  flow  rate  and  bulk  temperatures  at  the  inlet 
and  exit  of  the  test  section  were  measured  on  real  time  and 
the  output  signals  were  transmitted  to  a high  speed  sampling 
recorder. 

3.1.1  Flow  Pattern  Observation 

The  first  series  of  the  drop  tower  experiment  was 
intended  to  observe  the  flow  patterns.  A rectangular 

transparent  boiling  test  section  allowed  visual  observation  of 
the  boiling  process  on  the  heating  surface  and  flow  pattern 
development  along  the  flow  direction.  Figure  3-2  shows 
details  of  the  boiling  test  section  and  its  cross  section 
view.  Four  pyrex  plates  pressed  by  stainless  steel  plates 
were  fastened  together  to  form  a 5mmx5mmx500mm  square  channel, 
inside  which  transparent  electric  heaters  were  coated.  The 
heaters  consist  of  a thin  film  made  of  Cl2Sn  and  its  total 
resistivity  of  four  surfaces  was  308Q.  A llOV  of  dc  power  was 
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Figure  3-2  Test  section  used  in  drop  towerexperiments . 
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applied  to  the  heater  to  produce  55. 4W  of  total  electrical 
input  into  the  test  section.  It  was  found  later  that  an 
average  of  54%  of  electrical  power  was  used  to  heat  the  fluid 
passing  through  the  channel.  The  resulting  heat  flux  in  the 
test  section  was  2.9kW/m^.  The  exit  qualities  ranged  from 
-4.5%  to  17.5%  in  the  present  study.  The  electrical  power  for 
the  test  section  and  preheater  was  continuously  supplied  by 
extension  cables  during  the  free  fall. 

Outside  the  test  section,  a 30mmxl20mm  surface  mirror  was 
attached  beside  the  flow  channel.  The  reflected  image  of  flow 
patterns  through  20mmxl00mm  windows  were  photographed  with  a 
35mm  still  camera  with  back-lit  configuration.  The  camera, 
triggered  with  a remote  controller  at  the  top  of  the  drop 
tower,  was  motor-driven  and  capable  of  taking  5 frames  every 
0.25  seconds  during  1.25  seconds  of  free  fall. 

3.1.2  Pressure  Drop  Measurement 

The  main  objective  of  the  second  series  of  the  drop  tower 
experiment  was  to  measure  pressure  drops  under  microgravity 
conditions.  The  same  flow  loop  was  used  except  for  the  test 
section.  The  test  section  used  in  the  pressure  drop 
measurement  consists  of  a 4mm  I.D.  SUS  tube  wrapped  uniformly 
with  a nichrome  wire  whose  total  resistance  was  13.0  ohms.  A 
26V  dc  power  was  applied  to  this  heater  to  generate  52W  of 
electrical  power  input.  Heat  loss  to  the  ambient  was  measured 
by  reading  the  temperature  drop  across  the  insulator  around 
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the  test  section.  Approximately  12%  of  heat  flux  was  lost  to 
the  ambient  and  the  resulting  net  heat  flux  transferred  to  the 
flow  was  6.63kW/m^.  The  test  section  was  arranged  vertically 
to  reduce  the  transient  time  from  l.Og  to  0.02g. 

A diaphragm-type  differential  pressure  transducer  was 
used  to  measure  pressure  drops  along  the  test  section. 
Pressure  taps  were  made  at  the  inlet  and  exit  of  the  channel 
and  connected  to  the  high  pressure  and  low  pressure  ports  on 
the  transducer  by  4mm  I.D.  flexible  tubes  filled  with  liquid 
Freon  R113.  Under  variable-gravity  conditions,  the  static 
pressure  inside  the  leading  tubes  needs  to  be  considered  since 
it  also  changes  with  gravity.  The  differential  pressure  was 
sampled  every  2 msec,  during  the  free  fall  and  recorded  in  the 
same  high  speed  sampling  recorder  used  in  the  flow  pattern 
experiment . 

The  above  two-phase  loop  was  mounted  on  a two-staged, 
500mm  diameter  cylindrical  drop  module  with  four  5cm  diameter, 
Im  long  spikes  underneath  to  mitigate  the  impact  of  landing 
which  was  expected  to  be  on  the  order  of  50g.  The  velocity  of 
the  module  at  the  end  of  the  free  fall  was  calculated  to  be 
19.4m/s.  Upon  landing,  the  four  spikes  were  almost  completely 
buried  into  the  sand.  The  drag  shield  hanging  from  the  module 
was  used  to  reduce  the  air  drag  induced  during  the  free  fall. 
The  test  module  with  the  fluid  loop  mounted  is  shown  in  Figure 
3-3 (a).  The  total  weight  of  the  test  module  was  87kg.  The 
magnitude  of  gravity  during  the  free  fall  was  previously 
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A 15m  drop  tower  facility  and  flow  boiling 
loop  mounted  on  a drop  module. 
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measured  with  an  accelerometer  and  found  to  be  0.02g.  The 
experiment  was  conducted  in  a 15m  drop  tower  facility  owned  by 
Ishikawajima  Harima  Heavy  Industry  Inc.,  Yokohama,  Japan,  in 
June  and  July,  1989.  The  free  fall  distance  of  the  drop  tower 
facility  was  12m  (Figure  3-3 (b) ) and  allowed  a microgravity 
duration  of  1.25  seconds. 

3.1.2  Pressure  Drop  Measurement 

First,  the  pump  and  measurement  systems  were  turned  on 
and  the  flow  was  initiated.  After  leaks  on  the  flow  lines 
were  checked,  the  voltage  on  the  preheater  was  gradually 
increased  until  a desired  inlet  temperature  was  attained. 
Then,  the  main  heater  was  turned  on  and  boiling  was  started. 
It  took  approximately  3-4  minutes  before  the  boiling  was 
stabilized  and  a steady  state  was  reached.  Meanwhile,  the 
batteries  of  the  camera  motor  and  remote  controller  were 
checked.  Using  the  built-in  crane,  the  operator  lifted  up  the 
drop  module  to  the  top  of  the  tower  slowly  with  maintaining 
vertical  attitude.  About  30  seconds  were  waited  until  the 
swinging  motion  of  the  module  completely  damped  out.  Then, 
the  keblar  wire  hanging  the  module  was  cut  either  manually  or 
automatically  with  a video-monitored  scissors. 
Simultaneously,  the  motor-driven  camera  was  triggered  and 
photographed  5 to  6 frames  of  flow  patterns  during  the  free 
fall  consecutively.  After  landing,  the  heater  power  was 
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turned  off  and  the  data  stored  in  the  recorder  were  printed 
out. 


3.2  Leariet  Flight  Experiment 

The  experimental  package,  as  shown  in  Figure  3-4 (a)  and 
(b) , consists  of  two  racks  that  meet  dimensional  requirements 
of  the  NASA  Learjet  and  a boiling  test  section.  Rack  No.l 
contained  the  entire  flow  boiling  apparatus,  while  Rack  No. 2 
contained  a computer  and  data  acquisition  system.  The  test 
section  stood  vertically  between  the  two  racks.  The  flow 
boiling  apparatus  mainly  consisted  of  a boiling  test  section, 
two  0.178mO.D.  acrylic  plastic  pistons,  coiled  preheater, 
condenser,  high  pressure  air  tank,  and  flow  meter  system. 
Schematic  diagram  of  the  experimental  apparatus  is  shown  in 
Figure  3-5.  Specifications  of  each  component  are  described  in 
the  ensuing  text. 

Test  section.  Figure  3-6  shows  a detailed  schematic  of 
the  test  section.  The  test  section  consisted  of  a 1.28cm 
diameter,  0.56m  long  (0.48m  for  heated  length)  copper  tube 
with  0.711mm  wall  thickness,  a WATLOW  500W  coiled  cable 
heater,  and  one  flow  pattern  window  at  the  end  of  the  copper 
tube.  The  mineral-insulated,  stainless-sheath  heating  cable 
with  3 . 2mm  diameter  was  coiled  around  the  copper  tube  and 
covered  approximately  24%  of  the  surface.  The  cable  was 
clamped  every  5cm  to  secure  solid  contact  with  the  tube.  The 
test  section  was  covered  with  12.3mm  thickness  ceramic  fiber 
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Figure  3-4  Schematics  of  Lear jet  flight  research  instrument. 
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insulating  strips  to  minimize  the  heat  loss  to  ambient.  Two 
sets  of  thermocouples  were  installed  inside  and  outside  the 
insulator  to  measure  the  heat  loss.  The  measurement 
showed  that  about  15%  of  heater  power  was  dissipated  through 
the  insulator. 

Two  12.8mm  three-way  valves  were  connected  at  the  ends  of 
the  test  section  to  change  the  orientation  of  the  flow.  The 
entire  test  section  was  fastened  on  a 4 . 8mm  thickness  aluminum 
plate  arranged  vertically.  An  array  of  17  E-type 

thermocouples  equally  spaced  on  the  copper  tube  measured 
outside  wall  temperature  change  during  the  experiment.  The 
bulk  fluid  temperatures  were  measured  by  1.58mm  diameter 
stainless  sheathed  thermocouples  mounted  on  both  ends  of  the 
test  section.  The  thermocouple  junctions  at  the  tip  of  the 
sheath  extended  6.4mm  into  the  flow.  The  accuracy  of  the 
temperature  measurement  determined  by  calibration  was  on  the 
order  of  ±0.5°C. 

Two  OMEGA  gauge  pressure  transducers  (PX800-030GV  and 
PX800-050GV)  were  mounted  at  the  inlet  and  outlet  of  the  test 
section  perpendicular  to  flow  to  measure  the  system  static 
pressure  as  well  as  the  pressure  drop  across  the  test  section. 
Nominal  accuracy  of  the  transducers  is  ±0.1%,  or  on  the  order 
of  ±27.6Pa  for  both  transducers.  The  voltage  output  signal 
was  digitized  and  converted  into  psig  unit  through  an  ACCES 
12-bit  data  acquisition  system. 
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A mixture  of  vapor  and  liquid  passed  through  the  flow 
pattern  window  at  the  exit  of  the  test  section.  The  flow 
channel  was  formed  between  two  transparent  acrylic  blocks 
2.2mm  apart  from  each  other.  The  blocks  were  fastened  with 
screws  and  sealed  with  Nitrile  seals.  This  thin  flow  area  was 
designed  to  prevent  overlapping  of  vapor  phase  distributed  in 
the  liquid  and  provides  a 4x5cm  view  for  flow  pattern 
observation.  Flow  patterns  were  recorded  by  a video  camera  on 
high  speed  mode. 

Plastic  pistons  and  pressurized  air  tank.  Before  the 
test,  about  9.8  liters  of  liquid  Freon  113  was  contained  in 
Piston  1.  The  pistons  were  made  of  3.2mm  acrylic  plastic 
cylinder  and  19.1mm  end  plates  hermetically  sealed  with 
nitrile  0-rings.  The  interior  of  the  cylinder  was  divided  by 
a free  piston  which  also  served  as  a boundary  between  liquid 
and  pressurized  air (Figure  3-7) . Three  0-rings  on  the  piston 
were  lubricated  with  NOSOL  Uta  350,  a product  of  Kluber 
Lubrication  which  is  compatible  with  Freon  113 . 

Piston  1 was  pressurized  up  to  68.9kPa  with  compressed 
air  stored  in  a high  pressure  tank.  A regulator  controls  air 
flew  into  the  cylinders  and  maintained  pressure  constant  even 
in  the  expanded  volume.  Figure  3-8  shows  the  calibration 
curve  for  compressed  air  pressure  vs.  flow  rate.  Since  the 
flow  rate  increased  linearly  with  the  air  pressure,  the 
required  flow  rate  was  obtained  by  controlling  the  air 


pressure. 


49 


05 

O 

(/) 

O) 


05 


O CQ 


: ::;f.;:;X.;:;X;X-Xv  1 

E 


3^ 


(O 

Q. 

O 


(O 

05  (1) 

^ c 

- 

CO  .9 


T3 

05 

0 


O 

to 

• MM 

Q. 

'E' 

O 

< 


t 


Figure  3-7  A 0.178mO.D.  air-driven  plastic  piston. 
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Figure  3-8 


Air  driving  pressure  vs.  liquid  flow  rate 
for  the  piston  system. 


51 


Preheater  and  condenser.  The  preheater,  located  in  the 
upstream  of  the  test  section,  consisted  of  9.5mm  diameter, 
3.4m  long  coiled  copper  tube  and  a 600W  heavily  insulated 
heating  tape.  Although  the  preheater  was  originally  installed 
to  reduce  inlet  subcooling,  we  found  during  the  ground 
operation  that  start-up  period  of  the  preheater  was  about  3 
minutes,  which  would  sacrifice  the  number  of  trajectories 
available  during  the  flight.  Therefore,  the  preheater  was  not 
activated  during  the  present  experiments  and  the  test 
conditions  were  confined  to  high  subcoolings. 

After  leaving  the  test  section,  a mixture  of  vapor  and 
liquid  was  completely  condensed  inside  a 9.5mm  diameter,  2.3m 
copper  tube  coiled  and  immersed  in  ice  and  water  contained  in 
a well-insulated  therm  flask.  Inlet  and  outlet  temperatures 
were  constantly  monitored  and  recorded  throughout  the 
experiment.  Although  the  temperature  of  the  iced  water  was 
not  monitored,  previous  ground  tests  showed  that  it  was 
approximately  1“C. 

Flow  meter  and  signal  conditioner.  An  OMEGA  turbine  flow 
meter  (FTB-201) , which  was  insensitive  to  the  gravity  change, 
measured  the  liquid  flow  rate  and  sent  voltage  signal  to  a 
data  acquisition  system  through  a signal  conditioner  (FLSC- 
28)  . The  flow  meter  and  signal  conditioner  were  calibrated 
separately  before  the  experiments  and  showed  an  accuracy  on 
the  order  of  ±0. 126cc/sec. 
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Electrical  controls.  Rack  No. 2 contained  a power  control 
unit  and  a data  acquisition  system.  The  power  control  unit 
for  the  main  heater  was  housed  in  a 50.8cm  long  by  12.7cm  high 
by  7.6cm  deep  aluminum  enclosure.  The  control  of  power  was 
made  by  switching  the  circuit  with  different  resisters.  An 
ammeter  indicated  the  current  provided  to  the  main  heater 
(Figure  3-9) . The  controller  draw  the  nominal  28V  dc  power 
from  the  aircraft  power  system  and  provided  regulated  current 
to  the  main  heater.  The  controller  was  protected  by  a 20  amp 
fuse  on  the  power  input  line.  A relay  switch  activated  by  12V 
dc  was  used  to  prevent  electrical  shock. 

The  personal  computer  for  data  acquisition  and  a light 
for  flow  pattern  recording  were  connected  directly  to  the 
nominal  60Hz  115VAC  power  system  on  board. 

Data  accmisition  system.  The  data  acquisition  system 
included  an  ACCES  12bit  A/D  converter  installed  inside  an  IBM 
compatible  286  personal  computer  and  two  16  channel 
multiplexers.  All  32  channels  were  assigned  to  properties 
such  as  temperatures,  pressures,  and  flow  rate,  and  scanned 
these  data  at  interval  of  about  0.7  seconds  (Table  3-1).  The 
minimum  time  interval  was  determined  by  the  performance  of  the 
host  computer.  On  the  other  hand,  the  channel -to-channel 
scanning  interval  was  on  the  order  of  a micro  second.  Since 
this  interval  was  much  smaller  than  the  transient  time  of 
thermal  properties,  it  was  regarded  that  the  properties  were 
practically  sampled  at  the  same  time.  The  outputs  from 
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Figure  3-9  Electrical  power  and  interconnect  diagram. 
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Channel  assignments  on  the  multiplexer  card. 


Gain 

Properties 

Unit 

— 

Thermocouple  cold  junction 

(-) 

xlOO 

TWl  (Inlet  Wall  temperature) 

(C) 

xlOO 

TW2 

(C) 

XlOO 

TW3 

(C) 

XlOO 

TW4 

(C) 

XlOO 

TW5 

(C) 

XlOO 

TW6 

(C) 

XlOO 

TW7 

(C) 

XlOO 

TW8 

(C) 

XlOO 

TW9 

(C) 

XlOO 

TWIO 

(C) 

XlOO 

TWll 

(C) 

XlOO 

TW12 

(C) 

XlOO 

TWl  3 

(C) 

XlOO 

TWl  4 

(C) 

XlOO 

TWl  5 

(C) 

XlOO 

TWl  6 

(C) 

XlOO 

TW17  (Exit  Wall  temperature) 

(C) 

XlOO 

Condenser  inlet  temperature 

(C) 

XlOO 

Condenser  exit  temperature 

(C) 

XlOO 

Bulk  temperature  at  the  inlet 

(C) 

XlOO 

Bulk  temperature  at  the  exit 

(C) 

XlOO 

Insulator  inside  temperature 

(C) 

XlOO 

Insulator  outside  temperature 

(C) 

XlOO 

Preheater  inlet  temperature 

(C) 

XlOO 

Preheater  exit  temperature 

(C) 

XlOO 

Cylinder  1 exit  temperature 

(C) 

XlOO 

Preheater  inlet  temperature 

(C) 

x50 

Air  pressure  in  Cylinder  1 

(psig) 

XlOO 

Test  section  inlet  pressure 

(psig) 

XlOO 

Test  section  exit  pressure 

(psig) 

xl 

Volume  flow  rate 

(GPM) 
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thermocouples  and  pressure  transducers  on  the  order  of  1/1000 
volt  was  amplified  100  times  in  the  multiplexers.  The 
calibrated  data  was  stored  in  a 3 2 OK  virtual  memory. 
Resolution  of  the  present  data  acquisition  system  was  2.44mV. 
The  corresponding  resolutions  for  the  temperature, 
pressure, and  flow  rate  were  0.41‘’C,  48.1Pa,  and  0.132cc/sec, 
respectively . 


CHAPTER  4 

FLOW  PATTERN,  VOID  FRACTION,  AND  PRESSURE  DROP  UNDER 

MICROGRAVITY  CONDITIONS 

Selection  of  appropriate  models  for  hydrodynamic  and  heat 
transfer  analyses  in  two-phase  flow  depends  on  the  prediction 
of  topology  of  the  flow  or  commonly  known  as  flow  patterns. 
The  particular  flow  pattern  is  determined  by  the  system 
pressure,  flow  rate,  heat  flux,  channel  geometry,  and  gravity 
level.  The  influence  of  the  gravity  on  flow  patterns  had  not 
been  intensively  discussed  until  the  application  of  two-phase 
flow  in  space  became  a practical  issue.  However,  it  is  easily 
imagined  that  the  buoyancy- induced  motion  of  the  vapor  phase 
would  become  weak  or  disappear  in  microgravity.  The  first 
task  of  the  present  study  is,  therefore,  to  identify  the  flow 
patterns  under  microgravity  conditions  and  obtain  an  insight 
into  the  influence  of  gravity  on  two-phase  flow  structure. 

In  this  chapter,  the  results  of  flow  pattern  observations 
in  drop  tower  experiments  and  Lear jet  flight  experiments  are 
presented.  Flow  patterns  in  microgravity  were  compared  to 
those  obtained  in  normal  gravity  under  similar  flow  conditions 
and  their  differences  were  discussed.  Next,  results  from  the 
pressure  drop  measurement  in  drop  tower  experiments  were 
reported.  The  influence  of  gravity  on  each  component  of 
pressure  drop  was  examined. 
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4 . 1 Observation  of  Flow  Patterns 

Four  test  conditions  were  chosen  to  study  flow  patterns 
and  boiling  evolution  in  the  channel.  The  experimental 
conditions  in  these  tests  are  summarized  in  Table  4-1.  Flow 
rates  and  subcoolings  in  the  tests  were  37 . 7kg/m^/sec  and 
114 . 4kg/m^/sec,  and  4.6K  and  17.  OK,  respectively.  Two 
relatively  high  and  low  subcoolings  were  selected  to 
investigate  the  influence  of  gravity  on  bubble  growth  in  the 
superheated  layer.  The  test  section  was  arranged  vertically 
in  the  test  module  to  minimize  the  transient  effects  due  to 
gravity  change. 

The  experimental  domain  of  these  tests  may  be 
characterized  by  Froude  number,  Weber  number,  and  Bond  number. 
These  dimensionless  groups  are  defined  as 


Fr  = 

(4-1) 

gD 

We  = 

(4-2) 

o/D 

(p,-p,)gD 

(4-3) 

a/D 

where  Ufo  denotes  the  liquid  velocity  at  the  inlet  of  the  test 
section,  g the  gravity  constant  on  the  ground  (9.8m/s^),  D the 
diameter  of  the  tube,  and  o the  surface  tension.  Fr  is  the 
ratio  of  the  inertia  to  the  gravity  force,  while  We  is  the 
ratio  of  the  inertia  to  the  surface  tension  force.  Bo  is 


Table  4-1  Test  conditions  for  flow  pattern  observation  in 

drop  tower  experiments. 
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defined  as  the  ratio  of  gravity  to  surface  tension  forces. 
Gravitational  acceleration  reduces  to  about  1/100  of  normal 
gravity  for  airplane  experiments,  and  about  1/50  for  drop 
tower  experiments.  Figure  4-1  shows  the  variation  of  Fr,  We, 
and  Bo  for  the  drop  tower  experiments.  Under  normal  gravity 
conditions,  the  gravity  force  is  significantly  larger  than  the 
other  forces.  Under  microgravity  conditions,  the  magnitudes 
of  surface  tension  and  liquid  inertia  become  comparable  to 
that  of  gravity.  Bo  for  normal  and  microgravity  conditions 
were  20.4  and  0.41,  respectively,  indicating  that  the  surface 
tension  became  as  important  as  the  gravity  force.  The  above 
result  shows  that  0.02g  of  microgravity  level  in  the  drop 
tower  experiment  was  capable  of  successfully  simulating  the 
hydrodynamic  aspect  of  microgravity  two-phase  flow  phenomena. 

During  the  drop  tower  experiment,  growth  of  flow  patterns 
in  a uniformly  heated  transparent  channel  was  photographed 
every  0.25  seconds  with  a motor-driven  35mm  still  camera 
during  the  1.25  seconds  of  free  fall.  Schematic  diagrams  and 
photographs  of  typical  flow  patterns  observed  under  micro  and 
normal  gravity  conditions  are  shown  in  Figure  4-2 (a)  and  (b) . 
There  were  three  distinguishable  flow  patterns  in  vertical 
upward  flow  in  normal  gravity;  bubbly,  slug,  and  annular  flow. 
Under  microgravity  conditions,  three  flow  patterns  were  also 
identified;  although  their  details  were  quite  different  from 
the  corresponding  patterns  in  normal  gravity. 


0.013<Fr<5.89 
0.260  < We  < 2.34 
0.408  < Bo  < 20.4 


Gravity  dominant 


Bo 


We 


Figure  4-1 


Dominant  forces  under  the  conditions  of 
the  drop  tower  experiment. 
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Figure  4-2  Characteristic  flow  patterns  in  a vertical 

tube  under  various  acceleration  levels. 
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In  the  bubbly  flow  region,  two  types  of  bubbles  exist  in 
terms  of  the  way  they  are  formed  and  the  process  they 
contribute  to  the  heat  transfer.  Those  observed  on  the  wall 
are  nucleation  bubbles  as  a result  of  direct  evaporation  from 
a thin  microlayer.  They  transfer  significant  amount  of  heat 
from  the  wall.  The  size  of  these  bubbles  are  determined  by 
hydrodynamic  force  balance  and  local  thermal  conditions  on  the 
nucleation  site.  Relatively  large  bubbles  observed  in  the 
liquid  core  region  are  formed  after  several  detached 
nucleation  bubbles  coalesce  into  a single  larger  bubble  in  the 
course  of  their  migration  into  the  flow.  In  vertical  flow 
systems  where  the  buoyancy  force  acts  in  the  direction  of  the 
flow  axis,  bubble  motion  in  the  liquid  is  accompanied  by 
momentum  and  heat  exchange  on  a local  level . This  process 
causes  further  turbulence  even  in  laminar  flow  and  promotes 
heat  transfer. 

The  photographs  in  Figure  4-2  indicates  that  the  size  of 
nucleation  bubbles  was  similar  under  normal  and  microgravity 
conditions,  whereas  the  size  of  coalesced  bubbles  in 
microgravity  was  approximately  3-4  times  larger  than  under 
normal  gravity  conditions.  It  is  also  shown  that  the 
nucleation  bubbles  were  constantly  removed  by  the  liquid  flow 
even  under  microgravity  conditions  without  assistance  of  the 
buoyancy  force.  Bubble  growth  is  mainly  controlled  by  the 
thermal  condition  of  a thin  superheated  layer  adjacent  to  the 
heating  wall.  On  the  other  hand,  the  formation  of  coalesced 
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bubbles  is  governed  by  relative  motion  of  vapor  phase.  The 
almost  spherical  shape  of  each  bubble  indicates  that  the 
relative  velocity  was  extremely  small  or  zero. 

It  is  considered  that  the  size  of  a coalesced  bubble  is 
related  to  the  Helmholtz  instability  of  the  liquid-vapor 
interface  having  relative  motion.  The  wavelength  of  an 
interface  that  leads  to  instability  is  a function  of  the 
Laplace  constant  defined  as 


2 TZ 

N 


o 

(Pf-Pg)  g 


(4-4) 


where  Xq  is  the  wave  length  to  provide  a minimum  slip  velocity 
maintaining  interface  stability.  As  the  acceleration  level  g 
becomes  smaller,  A,q  becomes  longer.  Under  the  gravity  level 
of  0.02g,  Iq  is  7.1  times  longer  than  that  of  normal  gravity. 
There  seems  to  be  a certain  relation  between  the  enlarged 
bubble  size  and  the  extended  stabilizing  wavelength  under 
microgravity  conditions. 

Small  bubbles  were  observed  between  the  heating  surface 
and  the  vapor  core  in  both  slug  and  annular  flow  regimes. 
This  phenomenon  was  more  pronounced  in  the  annular  flow  regime 
in  which  small  bubbles  were  trapped  in  a liquid  film.  The 
existence  of  these  small  bubbles  significantly  perturbs  the 
surface  of  the  liquid  film  formed  on  the  heating  surface.  In 
the  annular  flow,  this  highly  disturbed,  irregular  liquid 
interface  causes  an  increase  in  pressure  drop. 
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The  results  of  flow  pattern  observation  along  with  the 
exit  quality  and  superficial  velocities  are  summarized  in 
Table  4-2.  The  term  "Frothy"  was  used  to  describe  a liquid 
film  containing  small  bubbles.  Flow  patterns  of  upper  30cm 
part  of  the  test  section  were  photographed  by  a motor-driven 
still  camera.  The  transition  of  flow  patterns  from  1-g  to 
0.02-g  under  the  four  experimental  conditions  are  shown  in 
Figure  4-3  (a)  through  (h) . Figure  4-3  (a)  and  (c)  for  the  exit 
qualities  of  17.1%  and  9.5%,  respectively,  indicate  that  once 
the  annular  flow  regime  is  established  in  normal  gravity,  the 
pattern  sustains  under  microgravity  conditions.  Figure  4-3  (e) 
and  (f)  shows  a transition  of  flow  pattern  from  bubbly  to  slug 
in  the  lower  part  and  from  slug  to  annular  flow  in  the  upper 
part  of  the  channel.  The  flow  pattern  transition  under  Case 
D0FP4  shown  in  Figure  4-3  (g)  and  (h)  suggest  that  the  shear 
stress  and  velocity  profiles  under  microgravity  conditions 
were  quite  different  from  those  in  normal  gravity.  Since  the 
thermal  equilibrium  quality  at  the  exit  for  this  case  was 
-1.7%,  the  bulk  temperature  of  the  liquid  flow  was  still  less 
than  the  saturation  temperature  in  normal  gravity.  The 
bubbles  tend  to  recondense  upon  touching  the  cooler  liquid 
flow.  As  a result,  most  bubbles  were  observed  near  the  wall 
region  rather  than  in  the  liquid  core. 

The  second  series  of  experiments  was  conducted  using  a 
NASA  Lear jet  which  provides  approximately  18  seconds  of 
microgravity  conditions (0. Olg) . The  flow  patterns  were 


Table  4-2  Vapor  volume  fraction  and  flow  patterns  in  drop  tower  experiments. 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP1  (z=30cm  - 40cm) 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP2  (2=40cm  - 50cm^ 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP2  (z=30cm  - 40cm) 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP3  (z=40cm  - 50cm) 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP3  (z=30cm  - 40cm) 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP4  (z=40cm  - 50cm^ 


Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 


Case:  D0FP4  (z=30cm  - 40cm) 
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Figure  4-3  Flow  pattern  transitions  during  the  free  fall. 
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recorded  with  a video  camera  at  the  exit  of  the  test  section 
through  two  transparent  plastic  blocks  between  which  was  a 
2 . 6mmx38 . 4mm  rectangular  channel.  This  rectangular  channel 
has  the  same  cross  sectional  area  as  the  circular  boiling  test 
section  attached  below.  This  prevents  a pressure  drop  at  the 
junction  and  is  intended  to  provide  two-dimensional  void 
distribution  across  the  cross  section.  A high-speed  mode, 
video  camera  was  used  to  investigate  continuous  motion  of 
phase  change  and  to  measure  bubble  velocity  in  the  dispersed 
flow.  After  several  trials,  a high-speed  mode  of  1/2000 
second  with  back-lit  setting  was  found  to  be  appropriate  to 
record  the  moving  boundaries  of  bubbles. 

During  the  experiment,  the  high-speed  images  were  blurred 
due  to  power  failure  of  a video  light.  The  video  camera  was 
switched  to  an  automatic  mode  and  recording  was  continued. 
The  recorded  images  qualitatively  showed  the  fundamental 
difference  of  vapor  phase  distribution  during  the  gravity 
transition.  The  observation  of  flow  patterns  in  subcooled 
flow  identified  a similar  change  of  vapor  phase  distribution 
to  that  observed  in  the  drop  tower  experiment.  For  Case  F1015 
with  the  exit  quality  of  -0.41%,  the  slug  flow  observed  in 
normal  gravity  changed  into  the  annular  flow  with  smooth 
surface  in  microgravity.  This  observation  agrees  with  the 
previous  result  from  D0FP3.  The  conditions  at  the  exit  of  the 
test  section  for  F1020  and  F1035  were  highly  subcooled. 
Formation  of  large  bubbles,  similar  to  those  observed  in  D0FP4 
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in  the  drop  tower  experiment,  was  also  observed  during  the 
microgravity  period. 

Through  the  two  series  of  microgravity  experiments, 
three  types  of  typical  flow  patterns  were  also  identified 
under  0.02g  and  O.Olg  microgravity  conditions.  The  flow 
pattern  in  the  subcooled  boiling  region  in  both  normal  and 
microgravity  was  classified  as  bubbly  flow,  although  the 
diameter  of  coalesced  bubbles  in  microgravity  was 
significantly  larger  than  that  in  normal  gravity.  The 
concentrated  bubbly  flow  and  slug  flow  in  normal  gravity  tends 
to  shift  into  the  annular  flow  in  microgravity.  The  influence 
of  gravity  appeared  less  significant  on  the  hydrodynamic 
aspect  of  annular  flow. 

Flow  pattern  mapping.  Observed  flow  patterns  in  the  drop 
tower  experiments  are  plotted  on  the  Dukler  reduced-gravity 
flow  pattern  map  as  a function  of  liquid  and  vapor  volumetric 
fluxes,  Ui3  and  U^3[53].  Among  various  flow  regime  maps [13], 
Dukler 's  map  was  chosen  because  the  volumetric  flux  was  a 
measurable  quantity  independent  of  gravity.  Figure  4-4 (a)  and 
(b)  show  the  distributions  of  flow  patterns  under  normal  and 
microgravity  conditions,  respectively.  Although  the  number  of 
data  sets  is  small,  each  flow  pattern  is  well  classified  on 
these  maps.  Large  bubbles  observed  in  the  low-quality  range 
were  also  classified  as  bubbly  flow. 

Transition  boundaries.  The  relationship  between  the  two 


velocities  is  given  by 


Uls  (m/s)  Uls  (m/s) 


(a)  g/gn=1.0 


(b)  g/gn=0.02 


Figure  4-4 


Flow  pattern  transitions  on  the  Dukler 
flow  regime  map. 
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U,s  = 


\ 


-1 


u - 

gs  r 


(4-5) 


where  Uj,  is  the  bubble  rising  velocity  and  expressed  as 


11=0.35 


(Pf-Pg)  gD 
pf 


\l/2 


(4-6) 


for  slug  flow,  while  Ub  for  bubbly  flow  is  given  by [51] 


\0.52 


K = 0.33 g 


0.76 


1.28 


(4-7) 


for  bubble  diameters  less  than  1mm.  For  low  pressures.  Eg.  (4- 
6)  is  reduced  to  0.35(g*  D)^^^.  Cq  in  Eg.  (4-5)  is  the 
distribution  parameter  defined  in  the  drift  flux  model.  A 
value  of  1.2  was  assumed  throughout  the  analysis  for 
simplicity.  The  transition  boundaries  from  bubbly  to  slug  and 
slug  to  annular  flow  were  experimentally  obtained  and  shown  in 
Figure  4-4.  The  critical  void  fractions  for  each  boundary  was 
calculated  from  Eg. (4-5) . The  value  of  Ub  was  obtained  by 
assuming  g=0.196m/s^  in  microgravity. 

The  boundaries  between  bubbly  and  slug  flow  are  given  by 
aBs=0. 538±0. 025  (l.Og),  and 

Obs=0.378±0.130  (0.02g), 

while  those  between  slug  and  annular  flow  are 

asA=0.757±0.085  (l.Og),  and 

asA=0.794±0.135  (0.02g). 

It  is  interesting  to  note  that  these  boundaries  are  close  even 
though  the  acceleration  is  changed.  The  void  fraction  is  a 
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function  of  slip  velocity  and  significantly  varies  when  the 
gravity  level  is  reduced.  However,  the  critical  void 
fractions  are  mainly  determined  by  the  geometry  of  flow 
channel  and  the  velocity  distribution  of  flow.  It  appears 
that  the  critical  void  fractions  are  independent  of  gravity 
change.  In  other  words,  the  same  flow  pattern  map  can  be  used 
for  variable  gravity  conditions  when  void  fraction  is  used  as 
a transition  criterion.  To  verify  the  above  hypothesis,  more 
data  must  be  collected  under  microgravity  conditions. 

Colin  et  al.  [14]  and  Dukler  et  al.  [53]  proposed 
different  values  of  0.2  and  0.45,  respectively,  as  a boundary 
between  bubbly  and  slug  flow  based  on  their  adiabatic, 
horizontal  two-phase  flow  experiments  carried  out  aboard  a 
reduced-gravity  aircraft.  The  present  analysis  qualitatively 
agreed  with  their  results.  This  implies  that  the  hydrodynamic 
mechanism  leading  to  bubble  coalescence  is  independent  of 
gravity  change.  The  boundary  for  the  transition  from  slug  to 
annular  flow  was  not  reported  since  their  system  was 
horizontal . 

Another  issue  neglected  in  the  above  analysis  is  the 
estimation  of  the  distribution  parameter,  Cq.  As  described  in 
the  previous  section,  the  radial  distribution  of  void  fraction 
is  quite  different  especially  in  the  subcooled  region. 
Selection  of  different  values  for  Cq  leads  to  different 
criteria  for  boundaries.  An  estimation  of  the  value  of  Cq 
under  microgravity  conditions  was  attempted  using  the  radial 
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distribution  of  the  void  fraction  and  discussed  in  the 
following  section. 


4.2  Void  Fraction 

The  void  fraction  was  measured  by  image  analysis  of  the 
same  photographs  used  for  identifying  flow  patterns.  The 
Macintosh  program  IMAGE  1.28P  converts  the  grayscaled  images 
into  binary  files.  By  adjusting  the  threshold,  the  area  of 
vapor  and  liquid  phases  were  converted  into  black  and  white 
pixels,  respectively.  Since  the  visual  observation  of  flow 
patterns  described  above  identified  a significant  increase  in 
void  fraction,  the  present  analysis  on  void  fraction  was 
focused  only  on  the  low  quality  range.  Cases  D0FP3  and  D0FP4 
with  exit  qualities  of  3.7%  and  -4.5%,  respectively,  were 
selected  for  this  purpose.  In  the  high-quality  region  whose 
flow  patterns  are  slug  and  annular,  the  irregular  boundaries 
of  vapor  phase  was  difficult  to  define  with  the  present 
method . 

The  binary  images  of  flow  patterns  are  shown  in  Figure  4- 
5.  IMAGE  1.28p  then  counted  the  number  of  bubbles  in  the 
designated  area  to  yield  the  average  diameter  of  each  bubble. 
The  accumulated  area  of  vapor  phase  was  divided  by  the  total 
area  of  the  photographed  region  to  obtain  the  average  vapor 
volume  fraction.  Since  the  flow  patterns  of  interest  are 
discrete  bubbly  and  slug  flow  in  the  low  quality  region,  the 
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Figure  4-5  Binary  images  of  vapor  phase  distributions  in 

the  radial  directions,  (a)  D0FP4  (z=40-50cm)  , 
(b)  D0FP4  (z=30-40cm) , (c)  D0FP3  (z=40-50cm) , 
and  (d)  D0FP3  (z=30-40cm) . 
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error  caused  by  the  overwrapped  vapor  phase  was  assumed 
negligible. 

The  void  fraction,  defined  as  the  area  ratio  of  vapor 
phase  across  the  cross  section,  is  calculated  from  the  volume 
fraction  using  the  following  relationship: 


(l-«) 

a 


where  s is  the  slip  ratio  given  by 


s 


\0.205 


Re 


-0.016 


(4-8) 


(4-9) 


for  normal  gravity  cases[25].  The  above  correlation  for  the 
slip  ratio  does  not  account  for  gravity  change.  The  average 
void  fraction  in  microgravity  was  significantly  larger  than 
that  with  normal  gravity  in  the  entire  low  quality  range. 
However,  the  void  fraction  profiles  for  normal  and 
microgravity  are  likely  to  merge  in  the  high-quality  region. 
This  is  confirmed  by  the  photographs  for  Cases  DOFPl  and 
D0FP2 , in  which  the  once-established  annular  or  concentrated 
slug  flow  patterns  in  normal  gravity  remained  as  they  were  in 
microgravity.  The  experimental  results  showed  that  the 
hydrodynamic  aspect  of  two-phase  flow  was  significantly 
influenced  by  gravity  change  in  the  low-quality  region.  In 
the  present  study,  the  drift  flux  model,  Sekoguchi's  model, 
and  the  slip  correlation  were  employed  to  predict  void 
fraction  profiles. 
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Drift  flux  model.  The  drift  flux  model  proposed  by 
Kroger  and  Zuber  includes  a gravity  term  explicitly  upon 
evaluating  the  weighted  mean  drift  velocity  <Ugj>.  The 
expression  for  <Ugj>  is  flow-pattern  dependent  and  different 
forms  are  proposed  for  particular  flow  patterns.  For  the  case 
of  one-dimensional,  vertical  up-flow  of  isolated  bubbles  with 
little  coalescence [ 32 ] : 


<u^.>  = 1.53  (l-a)2 


g] 


qg(pf-pg) 

pf' 


ll/4 


C(,=1.0  (4-10) 


while,  for  vertical-up  flow  in  the  bubbly-slug  region. 


<Ugj>  = 1.41 


og(pf-pg) 

pf 


1/4 


Co  =1.2 


(4-11) 


where  Cq  is  the  distribution  parameter  defined  as 


r = <«J> 

V-'O  I 

“ <axj> 


(4-12) 


At  this  moment,  the  distribution  parameter  characterizing  the 
radial  distribution  of  vapor  phase  in  microgravity  is  unknown. 
For  simplicity,  the  same  values  derived  for  normal  gravity 
conditions  were  used  throughout  the  present  analysis.  Note 
that  <Ugg>  is  proportional  to  g^^^  for  Eq.  (4-11)  , and  decreases 
by  62%  under  0.02g.  The  void  fraction  for  the  drift  flux 
model  is  given  by 
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(4-13) 


It  should  be  noted  that  in  an  extreme  case  of  zero  gravity 
(g=O.Om/s^),  the  above  equation  reduces  to 


If  Co  is  equal  to  1.0,  Eq.  (4.14)  yields  a homogeneous  model  in 
which  the  slip  ratio  is  equal  to  1.0.  The  slip  ratio  in  the 
drift  flux  model  is  given  by: 


In  the  above  equation,  Cq  is  the  distribution  factor  and 
indicates  the  degree  of  vapor  phase  distribution  across  the 
cross  section[47]. 

For  example,  a hypothetical  case  of  adiabatic,  laminar 
two-phase  flow  is  considered.  As  shown  in  Figure  4-6 (a),  the 
second  order  polynomials  are  assumed  for  the  volumetric  flux 
(j/jc)  and  the  void  fraction  profile  (a-a^/a^-K) 


(4-14) 


(4-15) 


2 


(4-16) 


2 


(4-17) 


CO  >1.0  CO  <1.0 
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Figure  4-20  Local  flow  conditions  near  the  bubble  detachment 

under  microgravity  conditions. 
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Here,  subscripts  c and  w denote  the  values  at  the  center  line 
and  the  wall,  respectively.  The  profile  of  the  total 
volumetric  flux  j approaches  that  of  single-phase  liquid  flow 
when  the  void  fraction  is  small.  For  simplicity,  a„=0  is 
assumed  for  the  adiabatic  flow.  Profile  (A)  in  Figure  4-6 (a) 
corresponds  to  the  case  where  the  vapor  phase  concentrates 
near  the  wall  and  simulates  the  flow  pattern  to  be  observed  in 
the  nucleate  boiling  region  in  a vertical  pipe.  The 
calculated  value  of  Cq  from  Eq. (4-12)  yields  0.667.  Profile 
(B)  is  the  flow  pattern  to  be  observed  in  slug  flow  and  the 
vapor  phase  concentrates  near  the  center  line  where  the  liquid 
velocity  is  fastest.  The  value  of  Cq  in  such  a case  is  1.60. 
When  the  vapor  phase  is  uniformly  distributed  across  the  tube, 
as  Profile  (C) , Cq  becomes  1.0. 

The  bubbly  flow  pattern  under  normal  gravity  conditions 
observed  in  the  drop  tower  experiments  showed  a concentration 
of  bubbles  near  the  wall,  while  the  concentration  of  vapor 
phase  shifted  toward  the  center  line  under  microgravity 
conditions.  Therefore,  it  is  inferred  that  the  value  of  Cq  in 
microgravity  may  be  larger  than  that  of  normal  gravity.  To 
estimate  the  value  of  Cq,  the  radial  void  fraction 
distributions  obtained  in  the  drop  tower  experiments  were 
approximated  with  the  6th-order  polynomials  and  calculated 
from  Eq.  (4-12)  . The  quadratic  form  of  j in  Eq.  (4-16)  was 
maintained  for  the  present  analysis.  The  distributions  of 
void  fraction  in  Case  D0FP4  were  selected  as  an  example 
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because  the  region  photographed  was  the  nucleate  boiling 
region  where  the  largest  difference  in  void  fraction  was 
observed  in  microgravity . The  constants  for  normal  and 
microgravity  void  fraction  profiles  are  as  follows 


i.oq 

0.02a 

ao 

+0.00474 

+0.43019 

ai 

-0.06263 

-1.34286 

a2 

+0.51416 

+19.9703 

^3 

-2.07937 

-105.916 

a4 

+7.78731 

+234.904 

as 

-10.7704 

-229.471 

ae 

+4.60682 

+81.4233 

The  values  of  Cq  for  normal  and  microgravity  conditions  are 
0.997  and  1.272,  respectively.  As  expected,  the  value  of  Cq 
in  normal  gravity  is  slightly  smaller  than  1.0  due  to  the  void 
concentration  near  the  wall.  However,  Cq  became  larger  than 
1.0  in  microgravity  mainly  due  to  coalescence  of  large  bubbles 
near  the  center  line.  According  to  Eg. (4-15) , the  average 
slip  ratio  becomes  larger  than  1.0  even  under  zero  gravity 
condition.  However,  this  prediction  contradicts  the  physical 
interpretation  of  bubble  movement  described  in  Figure  4-6 (b). 
Under  zero  gravity  conditions  in  the  nucleate  boiling  region, 
the  maximum  velocity  of  detached  bubbles  cannot  exceed  the 
liquid  velocity.  As  a result,  the  overall  slip  ratio  becomes 
less  than  1.0.  The  drift  flux  model  fails  to  predict  the  slip 
ratio  in  the  nucleate  boiling  region  when  the  gravitational 
acceleration  approaches  zero.  However,  for  simplicity,  the 
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normal-gravity  value  of  Cq  was  assumed  throughout  the  present 
calculation. 

The  void  fraction  profiles  for  Cases  D0FP3  and  D0FP4 
under : 

2 . 4g maximum  gravity  during  the  pull-up  maneuver  in 

the  Lear jet  flight  experiment, 

l.Og normal  gravity,  of  course, 

0.02g microgravity  level  in  the  drop  tower 

experiment , 

O.Olg microgravity  level  in  the  Learjet  flight 

experiment,  and 

1.0E-6g — residual  microgravity  level  experienced  in 

spacecraft  in  orbits. 

are  plotted  in  Figure  4-7,  4-8,  and  4-9,  respectively.  From 
Figure  4-7  (a),  one  can  notice  that  the  deviation  of  void 
fraction  from  the  l.Og  prediction  is  more  significant  in 


microgravity  than  at  higher  gravities.  Although  all  curves 
appear  to  merge  at  higher  qualities,  the  difference  in  void 
fraction  becomes  larger  as  the  quality  becomes  smaller.  Under 
normal  gravity  conditions,  the  measured  void  fractions  were 
significantly  lower  than  those  predicted  by  the  drift  flux 
model.  This  is  mainly  due  to  the  small  heat  flux  used  in  the 
present  experiment.  The  amount  of  heat  fluxes  used  in  Kroger 
and  Zuber's  experiment  was  about  5-15  times  larger  than  that 
of  the  present  experiment.  Besides,  the  material  consisting 
of  the  heating  surface  was  coated  glass  which  has  a fairly 
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Figure  4-9 


Void  fraction  profiles  predicted  by 
the  slip  model. 
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Void  fraction  profiles  predicted  by 
Sekoguchi's  model. 
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Figure  4-9  Void  fraction  profiles  predicted  by 

the  slip  model. 
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smooth  surface.  It  is  considered,  therefore,  that  the  number 
of  activated  nucleation  sites  in  the  present  study  was  much 
smaller  than  that  in  Kroger  and  Zuber's  experiment. 
Recondensation  of  vapor  phase  was  estimated  by  Ahmad  method 
and  included  in  the  calculation  of  exit  qualities.  The  void 
fraction  predicted  for  the  imaginary  case  of  10"®g  is 
significantly  larger  than  that  of  l.Og.  This  indicates  that 
the  slip  velocity  still  exists  under  0.02g  conditions  and  that 
assuming  s=l . 0 overestimates  the  void  fraction  in 
microgravity . 

Sekoauchi's  method.  For  Case  D0FP4,  the  exit  condition 
is  still  subcooled.  In  the  subcooled  region,  the  thermal 
equilibrium  quality  becomes  negative  and  correlations  used  in 
saturated  boiling  are  no  longer  valid  for  the  negative 
quality.  However,  vapor  phase  exists  even  in  this  region  and 
plays  a significant  role  for  heat  transfer.  Because  of  the 
existence  of  vapor  phase,  the  net  or  actual  quality  has  a 
positive  value. 

Sekoguchi,  Tanaka,  Esaki,  and  Imasaka  developed  a semi- 
theoretical  model  to  predict  the  actual  quality  in  the 
subcooled  boiling  region,  assuming  a temperature  and  velocity 
profile  in  the  radial  direction[33] . The  methodology  is 
summarized  in  Chapter  5.3.2  since  it  is  closely  related  to  the 
heat  transfer  problem.  His  model  was  particularly  selected 
here  because  its  theoretical  nature  provides  an  insight  into 
the  gravity  dependance  of  the  flow  and  heat  transfer  more 


93 


explicitly  than  other  entirely  empirical  correlations.  His 
model  solves  the  void  fraction  profile  in  the  axial  direction 
simultaneously  through  iteration  process.  For  the  void 
fraction,  however,  he  employed  Smith's  correlation  including 
a factor  derived  empirically  and  set  to  0.4 [34].  The  gravity 
effects  on  void  fraction  are  not  explicitly  shown  in  Smith's 
correlation. 

The  single-phase  heat  transfer  coefficient  consists  of 
laminar  flow  and  natural  convection  components.  The  latter  is 
a function  of  the  Grashof  number  and  proportional  to 
Through  the  gravity-dependent  single-phase  heat  transfer 
coefficient,  the  true  quality  as  well  as  void  fraction 
incorporates  the  gravity  effects  in  Sekoguchi's  model. 

In  the  saturated  region,  the  thermal  equilibrium  quality 
alone  was  used  for  the  above  calculation. 

The  resulting  void  fraction  profiles  under  variable 
gravity  conditions  are  shown  in  Figure  4-8 (a)  and  (b)  . For 
Case  D0FP3,  the  void  fraction  profiles  for  various  gravities 
merge  in  the  lower  quality  region  than  that  of  the  drift  flux 
model.  The  difference  of  the  predicted  void  fraction  in  0.02g 
and  that  of  l.Og  is  less  pronounced  compared  to  the  drift  flux 
model.  The  void  fraction  in  the  saturated  boiling  region  is 
still  overestimated  and  its  difference  is  larger  than  that  of 
the  drift  flux  model.  In  the  highly  subcooled  region,  the 
predictions  qualitatively  agreed  with  the  experimental  data  of 
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normal  gravity.  This  is  because  Sekoguchi's  method  was 
exclusively  developed  for  the  subcooled  boiling  analysis. 

Slip  correlations.  The  slip  ratio  under  normal  gravity 
conditions  is  estimated  by  Eg. (4-9) . For  the  present  system 
in  l.Og,  the  value  of  s was  2.63.  However,  Eg. (4-9)  is  unable 
to  predict  the  slip  ratio  under  microgravity  conditions. 
Three  values  of  slip  ratio (1.0,  3.0,  and  5.0)  were  arbitrarily 
chosen  and  substituted  into  the  following  eguation  to  obtain 
the  void  fraction: 


<a>  = 


1 


1+sil  -l) 

fp.l 

\X  ) 

IpJ 

(4-18) 


The  curves  are  plotted  in  Figure  4-9 (a)  and  (b) . The  result 
shows  that  the  slip  ratio  must  be  much  larger  to  match  the 
present  experimental  data.  Besides,  these  curves  tend  to  vary 
linearly  and  do  not  reflect  the  appropriate  void  fraction 
behavior  around  x=0.0.  In  the  subcooled  region,  the 
prediction  grossly  underestimates  the  experimental  data. 

Other  well-known  models  for  subcooled  boiling  consists  of 
Levy's  and  Larsen-Tong ' s models [35] [36] . However,  Levy's 
model  reguires  subcooling  at  the  bubble  detachment  point, 
which  was  unavailable  in  the  present  experiment.  Levy 
considered  the  hydrodynamic  force  balance  of  a detaching 
bubble  to  estimate  the  liguid  temperature  at  the  top  of  the 
bubble.  Since  the  magnitude  of  the  buoyancy  force  is  expected 
to  vary  in  microgravity,  the  liguid  temperature  at  the  top  of 


95 


the  detaching  bubble  is  unknown.  He  concluded  that  the 
buoyancy  force  on  bubble  detachment  was  negligible  through 
correlating  his  experimental  data  with  his  model.  Therefore, 
gravity  effects  on  boiling  heat  transfer  are  not  clearly 
represented  with  his  model  unless  the  force  balance  is 
refined.  Similarly,  Larsen-Tong ' s model  requires  information 
about  the  thickness  of  a bubble  boundary  layer  growing  on  the 
heating  surface  in  order  to  calculate  the  heat  flux 
transferred  to  the  liquid  core.  This  quantity  was  not 
measured  in  the  present  experiment. 

The  flow  boiling  experiment  in  microgravity  showed 
significant  differences  in  the  average  void  fraction  and  the 
void  distribution  in  the  radial  direction.  Among  the  models 
described  above,  the  drift  flux  model  is  capable  of 
incorporating  these  changes  in  microgravity  with  less  effort. 
This  model  includes  the  gravity  term  explicitly  as  a form  of 
bubble  rising  velocity  and  is  able  to  express  the  change  of 
void  fraction  distribution  in  the  radial  direction  by  choosing 
an  appropriate  distribution  parameter  Cq  for  variable  gravity 
conditions. 

Under  a gravity  level  of  0.02g,  all  of  the  above  models 
and  experimental  results  showed  a large  increase  in  volume 
fraction  in  the  subcooled  region.  As  will  be  stated  in 
Chapter  5,  the  onset  of  nucleate  boiling  seems  insensitive  to 
gravity  change.  These  facts  lead  to  a conclusion  that  the 
void  fraction  rapidly  increases  in  the  subcooled  region  under 
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microgravity  conditions,  rather  than  gradually  increasing  as 
under  normal  gravity  conditions.  Such  a sudden  increase  in 
void  fraction  at  the  beginning  of  the  two-phase  region  may 
cause  an  instability. 

The  size  of  a coalesced  bubble  flowing  with  liquid  is 
normally  determined  by  hydrodynamic  balance.  However,  growth 
of  such  large  bubbles  is  restricted  by  the  liquid  temperature 
profile  and  geometry  of  the  heating  surface.  The  above 
analysis  also  suggests  that  under  low  gravity  conditions  there 
is  a possibility  for  the  slip  ratio  to  become  less  than  1.0  in 
microgravity . 

4.3  Velocity  and  Temperature  Profiles 

As  a result  of  different  void  distributions  under 
microgravity  conditions,  the  velocity  profile  and  heat 
transfer  mechanism  are  expected  to  change,  especially  in  the 
low-quality  region.  The  histograms  in  Figure  4-5  show  the 
average  radial  distributions  of  vapor  phase  scanned  every  0.25 
seconds  and  integrated  in  the  axial  direction  for  Cases  D0FP3 
and  D0FP4.  The  integer  on  the  y-axis  is  the  count  of  black 
pixels  on  a single  axial  scanning  line  normalized  for  2®=256. 
In  the  bubbly  flow  regime  in  normal  gravity,  the  bubbles 
generated  on  the  wall  form  two  peaks  adjacent  to  the  wall, 
while  they  are  merged  into  one  broad  peak  at  the  center  in 
microgravity.  For  the  slug  flow  in  D0FP3,  the  radial 
distribution  had  not  changed  much,  although  the  average  volume 
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fraction  increased  and  the  flow  pattern  approached  the  annular 
flow  in  microgravity . 

An  attempt  was  made  to  reflect  the  above  change  in  flow 
patterns  on  the  analysis  of  two-phase  flow  under  microgravity 
conditions.  To  investigate  the  effect  of  phase  change  in 
microgravity  on  the  velocity  and  temperature  profiles,  the 
contribution  of  vapor  phase  to  the  momentum  and  heat  transfer 
must  be  specified.  Sato,  Sadatomi,  and  Sekoguchi  developed  a 
theoretical  model  to  include  the  turbulent  effect  of  bubbles 
for  dispersed  bubbly  flow  and  verified  their  theory  with  an 
adiabatic  vertical  flow  experiment  using  an  air-water 
mixture [16]  [37]  [38] . Their  analysis  is  useful  for  the  present 
study  since  the  contribution  of  the  turbulence  caused  by 
bubbles  is  considered  separately.  In  the  following  paragraph, 
their  theory  is  briefly  outlined. 

They  introduced  an  additional  turbulent  stress  due  to 
bubble  agitation  into  the  two-dimensional  shear  stress 
equation  as  follows 

T =pf(l-a)  (Vf +€'+€'') -^  (4-19) 

where  a is  the  local  void  fraction,  Ul  is  the  average  liquid 
velocity  in  the  flow  direction,  Vf  is  the  kinetic  viscosity, 
e'  is  the  inherent  eddy  diffusivity  of  the  turbulence,  and  €” 
is  the  eddy  diffusivity  due  to  bubble  agitation.  Using  non- 
dimensional  parameters.  Eg. (4-19)  is  given  by 
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du^*  _ 


X* 


dy*  (i-a)  (V, +€'+€")/ 


(4-20) 


where  each  parameter  in  the  above  eqpiation  is  expressed  as 


X 


♦ 


(4-21) 


Here,  x„  and  Ul*  denote  the  wall  shear  stress  and  the  friction 
velocity,  respectively.  The  shear  stress  distribution  x*  in  a 
vertical  tube  was  obtained  from  the  force  balance  between 
pressure  and  shear  friction.  The  resulting  form  of  the  shear 
stress  yields 


X 


/ 


1 


(4-22) 


where  B=gR/UL*^.  Note  that  x*  is  reduced  to  unity  when: 

a)  the  void  fraction  becomes  small  and  the  integrations  in 
Eq. (4-22)  approach  zero,  or 

b)  the  term  B becomes  zero  regardless  of  void  distribution 
as  the  gravity  approaches  zero. 

In  the  present  experiments,  the  residual  gravities  of  0.01  and 
0.02g  are  still  large  enough  to  modify  the  shear  stress 
distribution  and  cannot  be  neglected. 
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The  eddy  dif fusivities  for  fully  developed  turbulent 
bubbly  flow  are  expressed  as  follows 


€'  = 


l-exp(-^) 

I\ 


\2 


6 R*  3 R*  3 R* 


(4-23) 


where  the  damping  factor  ( 1-exp (yVR^)  was  added  to  account 
for  the  attenuation  of  turbulence  and  the  eddy  diffusivity 
near  the  wall.  Based  on  the  experimental  data,  Laufer[63] 
chose  0.4  and  16  for  k and  A^,  respectively. 

For  the  eddy  diffusivity  due  to  bubble  agitation,  they 
proposed  the  following  expression 


(4-24) 


where  is  the  empirical  constant  and  assumed  to  be  1.2,  dg 
the  bubble  diameter,  and  Ug  the  relative  velocity  of  bubbles 
to  the  liquid  velocity.  Eg. (4-20) , Eg. (4-22) , Eg. (4-23)  and 
Eg. (4-24)  were  solved  simultaneously  to  obtain  the  velocity 
profile  with  the  following  boundary  condition 

u^*  = 0 for  y*=0  (4-25) 

Correction  of  was  iterated  until  the  mass  flux  G matched 
the  given  value. 

A similar  approach  was  employed  to  obtain  the  temperature 
distribution  in  the  radial  direction  using  the  velocity 
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distribution  obtained  above.  The  dimensionless  form  of  the 
temperature  profile  is  given  by 


dy* 


-Inu 

2 


(1-a) 


.//« 


1 +Pr. 


e'  +e" 


f )\ 


(4-26) 


where  Tl*,  q”*,  and  Nu  are  the  dimensionless  temperature,  heat 
flux,  and  Nusselt  number,  respectively,  defined  as 


m • _ ( ^CL ) 


(4-27) 


\ (1-a)  u.rdr 


*=  g = -R  0 


I (1-a)  u,rdr 


(4-28) 


(4-29) 


The  two-phase  flow  heat  transfer  coefficient  hip  is  given  by 
q„'V  (T„-Tcl)  . Under  the  following  boundary  condition  at  the 
wall, 

Tl*=1,  for  y*  = 0 (4-30) 

Eg. (4-26)  and  Eq. (4-28)  were  solved  simultaneously  by  the  4th 
order  Runge-Kutta  method.  The  value  of  T„  was  initially 
guessed  and  iteration  was  repeated  until  Td,  converges  into  the 
measured  value.  To  deal  with  the  extremely  large  temperature 
gradient  near  the  wall  region,  the  radial  distance  was  divided 
into  165  meshes  using  a non-uniform  grid  system. 
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In  the  above  analysis,  the  void  fraction  profile  must  be 
given  as  input  data.  The  above  method  was  applied  to  the 
results  of  Case  D0FP4  obtained  in  the  first  series  of  drop 
tower  experiment.  The  flow  rate  was  114kg/mVsec  and  the 
channel  width  of  the  boiling  test  section  was  5mm.  The  heat 
flux  applied  to  the  0.5m  test  section  was  2.93kW/m^.  The 
measured  center  line  temperature  was  46.0°C  for  both  normal 
and  microgravity  conditions.  Laminar  flow  analysis  was  used 
since  the  liquid  Reynolds  number  was  1540  for  Case  D0FP4.  To 
evaluate  e”,  dg  and  Ug  need  to  be  specified.  The  average 
values  of  bubble  diameter  under  normal  and  microgravity 
conditions  were  read  from  the  photograph,  while  the  slip 
velocity  Ug  was  replaced  with  the  rising  velocity  of  a single 
bubble  to  reflect  the  gravity  effects.  The  values  of  dg  were 
0.0011m  and  0.0032m  for  normal (l.Og)  and  microgravity  (0.02g) 
conditions,  respectively. 

The  results  from  the  above  calculation  in  normal  and 
microgravity  for  Case  D0FP4  in  the  drop  tower  experiments  are 
shown  in  Figure  4-10  (a)  through  (e)  . Because  of  the  slip 
velocity  of  vapor  phase,  bubbles  with  small  diameters  tend  to 
concentrate  near  the  wall  under  normal  gravity.  On  the 
contrary,  coalescence  of  bubbles  seem  to  be  promoted,  forming 
larger  bubbles  under  microgravity  conditions.  The  average 
diameter  of  coalesced  bubbles  in  microgravity  at  the  same 
position  was  about  2.9  times  larger  than  those  of  normal 
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(a)  Void  fraction  distributions 
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(b)  Liquid  velocity  distributions 

Figure  4-10  Two-phase  flow  properties  under  1 . Og  and 

0.02g  for  Case  D0FP4  calculated  by 
Sato's  method. 


Temperature  (C)  Tau 
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(c)  Normalized  shear  stresses 


(d)  Liquid  temperature  distributions 


Figure  4-10 


continued 


(s/3vUi)  nu 
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(e)  Total  eddy  dif fusivities 


Figure  4-10  


continued. 
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gravity.  This  is  related  to  the  Helmholtz  instability  of  a 
free  surface  moving  with  a certain  velocity. 

Using  the  void  fraction  profile  in  Figure  4-10 (a),  the 
velocity  distributions  were  calculated  for  both  normal  and 
microgravity  conditions  (Figure  4-10 (b)).  Compared  to  the 
single-phase  laminar  flow  profile,  liquid  velocities  of  two- 
phase  flow  in  both  normal  and  microgravity  became  larger  owing 
to  concentration  of  voids  near  the  wall.  When  a slip  velocity 
existed,  the  velocity  profile  of  microgravity  two-phase  flow 
approached  that  of  turbulent  flow  even  though  the  flow  was 
laminar.  However,  the  velocity  profile  for  two-phase  bubbly 
flow  in  microgravity  without  the  slip  velocity  showed  that  the 
profile  was  similar  to  that  of  single-phase  flow.  Under  such 
conditions,  the  existence  of  bubbles  decreases  the  flow  area 
and  increases  the  liquid  velocity  if  the  mass  flux  is  kept 
constant.  Although  the  single-phase  heat  transfer  increases, 
no  contribution  is  made  to  the  generation  of  turbulence  by  the 
bubbles  under  zero-gravity  conditions.  Shear  stress 
distribution  in  microgravity  were  similar  to  those  of  single- 
phase laminar  flow  in  normal  gravity  (Figure  4-10 (c)). 
However,  the  profile  for  two-phase  flow  in  normal  gravity  was 
significantly  different  from  others.  This  is  because  the 
parameter  B in  Eq. (4-22)  includes  the  gravity  constant  g and 
the  profile  for  x*  approaches  that  of  single-phase  flow 
regardless  of  the  void  fraction  profile. 


Table  4-3  Computed  wall  shear  stresses  and  wall  temperatures  of  laminar  two-phase  flow 

for  Case  D0FP4 . 
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The  calculated  wall  shear  stresses  for  normal  and 
microgravity  two-phase  flow  are  shown  in  Table  4-3.  The  ratio 
of  the  wall  shear  stress  in  two-phase  flow  to  that  of  single- 
phase flow  in  normal  gravity  was  2.07.  Under  microgravity 
conditions,  this  ratio  became  less  than  2.07  when  Ug  was  less 
than  0.003m/s. 

Temperature  distribution  in  the  radial  direction  was 
obtained  using  the  velocity  profile  previously  computed. 
Figure  4-10 (d)  shows  the  temperature  profiles  for  both  normal 
and  microgravity  two-phase  flows.  The  temperature  profile  for 
single-phase  laminar  flow  is  also  shown  for  reference.  The 
calculated  wall  temperatures  for  microgravity  conditions  were 
higher  than  that  under  normal  gravity,  indicating  that  the 
heat  transfer  rate  was  deteriorated  in  microgravity.  The  wall 
temperatures  and  Nusselt  numbers  are  computed  and  listed  in 
Table  4-3.  The  ratio  of  Nu  for  two-phase  flow  to  that  of 
single-phase  flow  Nu/Nu^p  for  l.Og  and  0.02g  are: 
(hag+hfc)/hfo=20 . 3 for  l.Og,  and 

(hag+hfc)/hfo=10.4  for  0.02g  (Ub=0. 005m/s)  . 

where  hag  denotes  the  heat  transfer  coefficient  accounting  for 
bubble  agitation,  hf^  the  heat  transfer  coefficient  accounting 
for  single-phase  forced  convection,  and  hfo  the  heat  transfer 
coefficient  for  liguid  alone.  This  result  indicates  that  heat 
transfer  in  the  two-phase  flow  for  the  same  amount  of  heat 
flux  was  deteriorated  under  microgravity  conditions. 
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It  is  considered  that  the  turbulence  generated  by  the 
bubbles  were  weakened  in  microgravity  in  spite  of  the  increase 
in  bubble  size.  Figure  4-10 (e)  shows  the  intensity  of  total 
eddy  diffusivity  for  bubbly  flow.  Since  the  eddy  caused  by 
bubble  agitation  is  much  larger  than  the  kinetic  viscosity  Vl, 
the  total  eddy  diffusivity  follows  the  void  fraction  profile. 
However,  e"+Vf  in  microgravity  becomes  smaller  as  the  gravity 
constant  becomes  smaller.  Below  a certain  value  of  gravity 
constant,  e"+Vf  in  microgravity  becomes  smaller  than  that  of 
normal  gravity  although  the  average  void  fraction  in  O.Olg  is 
much  larger  than  that  of  normal  gravity.  Under  zero-gravity 
conditions,  the  slip  velocity  disappears  and  e"  becomes 
completely  zero. 

The  same  analysis  was  applied  to  Case  F1035  in  the 
Lear jet  flight  experiment.  The  liquid  Reynolds  number  for 
this  case  was  9616  and  falls  in  the  turbulent  flow  range.  The 
heat  flux  used  in  this  case  is  about  10  times  larger  than  that 
of  Case  D0FP4.  The  eddy  diffusivity  inherent  to  the  turbulent 
flow  e'  is  introduced  in  the  calculation.  Selection  of  Case 
F1035  is  based  on  the  fact  that  the  pressure  change  during  the 
gravity  transition  was  smallest.  Since  the  void  fraction 
profile  was  unavailable  for  this  series  of  experiments,  the 
same  void  fraction  profiles  of  Case  D0FP4  was  employed. 

The  resulting  velocity  profiles  are  shown  in  Figure  4- 
11(a).  Because  of  the  void  fraction  profile,  the  velocity 
near  the  wall  is  slightly  increased  in  normal  gravity.  Under 


Liquid  veiocity  (m/s) 
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(a)  Liquid  velocity  distributions 


(b)  Total  eddy  diffusivity 


Two-phase  flow  properties  under  1 . Og  and 
0.02g  for  Case  F1035  calculated  by 
Sato's  method. 


Figure  4-11 
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(c)  Eddy  dif fusivities  (l.Og) 


(d)  Eddy  dif fusivities  (O.Olg) 


Figure  4-11 


continued 
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(e)  Normalized  shear  stresses 


(f)  Liquid  temperature  profiles 


Figure  4-11 


continued 


Table  4-4  Computed  wall  shear  stresses  and  wall  temperatures  of  turbulent  two-phase 

flow  for  Case  F1035. 
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microgravity  conditions,  the  profile  is  similar  to  that  of 
single-phase  flow.  This  is  because  e"  becomes  smaller  and  e' 
became  larger  than  e"  in  microgravity  (Figure  4-11  (b)  and  (d) ) . 
However,  in  normal  gravity,  e"  and  €'  were  comparable  to  each 
other (Figure  4-11 (c)).  The  shear  stress  distributions  shown 
in  Figure  4-11 (e)  were  similar  to  those  of  Case  D0FP4,  except 
that  the  magnitude  of  x*  in  normal  gravity  was  smaller.  The 
calculated  wall  shear  stresses  for  normal  and  microgravity 
conditions  are  shown  in  Table  4-4.  The  wall  shear  stress  in 
zero-gravity  was  1.18  times  larger  than  that  of  normal 
gravity.  The  above  computation  indicates  that  larger 
frictional  pressure  drops  are  expected  in  microgravity  since 
the  relationship  between  the  wall  shear  stress  and  the 
frictional  pressure  drop  gradient  is  given  by 


rdPi 


dz 

2*,f 

' dP' 
dz 

(4-31) 


where  x„^2<|(  and  x„  are  the  wall  shear  stresses  for  two-phase 
and  single-phase  flow,  respectively.  The  measurement  of 
pressure  drops  in  microgravity  will  be  discussed  later  in 
Section  4.4. 

The  temperature  profiles  in  Figure  4-11 (f)  showed  a 
decrease  in  heat  transfer  rate  under  microgravity  conditions. 
The  calculated  wall  temperatures  under  normal  and  microgravity 
conditions  (Ub=0. 005m/s)  were  55.9°C  and  63.2"C,  respectively, 
for  the  measured  bulk  temperatures.  The  wall  temperature  for 
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the  single-phase  turbulent  flow  was  76.3°C.  The  ratios  of  the 
two-phase  bulk  convection  heat  transfer  to  the  single-phase 
convection  heat  transfer  for  l.Og  and  O.Olg  were: 
(hag+hfc)/hfo=l. 68  for  l.Og,  and 

(hag+hfc)/hfo=l. 37  for  O.Olg  (Ub=0. 005m/s)  . 

The  bulk  convection  heat  transfer  rate  of  two-phase  flow  in 
O.Olg  was  about  18.4%  smaller  than  that  of  l.Og.  The  above 
result  indicates  that  the  bulk  convection  of  two-phase  flow 
consisting  of  single-phase  flow  plus  bubble  agitation  was 
deteriorated  under  microgravity  conditions.  It  also  shows 
that  the  contribution  of  bubble  agitation  to  heat  transfer 
becomes  less  significant  in  microgravity. 

The  total  heat  transfer  rate  including  evaporation  was 
measured  in  the  Lear jet  flight  experiment  and  their  results 
are  discussed  in  Chapter  5. 

4 . 4 Pressure  Drop 

The  pressure  drop  for  two-phase  flow  consists  of 
gravitational,  accelerational,  and  frictional  components. 
Among  the  three  components,  frictional  pressure  drop  is  of 
paramount  importance  since  its  gradient  is  directly  related  to 
the  wall  shear  stress,  hence  heat  transfer.  In  this  chapter, 
the  results  from  pressure  drop  measurements  conducted  in  the 
drop  tower  facility  is  reported  and  compared  with  the 
prediction  by  a homogeneous  model.  Then,  the  influence  of 
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(P1  -P2)  = APread  + pf  (h1  +h2)  g 
= AP1  + AP2  + pf  g H 


Figure  4-12 


Pressure  drop  measurement  along  a vertical 
test  section  in  the  drop  tower  experiment. 
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gravity  on  the  accelerational  and  frictional  pressure  drops  is 
discussed. 

The  pressure  drop  along  the  axis  of  the  4ininI.D.x500inm  SUS 
tube  was  measured  by  a diaphragm-type  differential  pressure 
transducer.  Pressure  taps  were  attached  perpendicular  to  the 
flow  at  the  exit  and  inlet  of  the  test  section.  Two  lead 
lines  filled  with  liquid  Freon  113  were  connected  to  the  high 
and  low  pressure  ports  of  the  transducer  elevated  by  hi  from 
the  inlet  of  the  test  section  (Figure  4-12)  . The  output  of  the 
transducer  includes  the  static  pressure  of  the  liquid  inside 
the  vertical  lead  lines.  The  total  pressure  drop  in  the  test 
section  (P1-P2)  is  given  by 

(P,-P,)  (4-32) 

As  long  as  the  lead  lines  are  filled  with  liquid  at  the  same 
temperature,  the  elevation  of  the  pressure  transducer  does  not 
affect  the  measurement.  Under  variable  gravity  conditions  for 
vertical  flow  systems,  however,  the  static  pressure  changes 
significantly.  Due  to  the  fact  that  the  static  pressure  in 
l.Og  for  the  drop  tower  experiment  is  approximately  50  times 
larger  than  that  of  0.02g,  continuous  use  of  the  same  range  in 
the  pressure  transducer  for  normal  and  microgravity  conditions 
was  impossible.  A finer  range  was  selected  for  measurement  of 
the  total  pressure  drop  in  microgravity.  The  output  signal 
ranging  from  0.96V  to  5.0V  was  recorded  in  a high  speed 
sampling  recorder  every  2 msec,  in  real  time.  The  voltage 
output  became  0.96V  when  the  differential  pressure  was  zero. 


Table  4-5  Test  conditions  for  pressure  drop  measurement  in  drop  tower  experiments. 
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Calibration  of  the  pressure  transducer  was  conducted  under 
normal  gravity  conditions  using  a column  of  Freon  113  at  room 
temperature . 

The  test  conditions  of  the  pressure  drop  measurements  are 
summarized  in  Table  4-5.  With  heat  flux  constant,  pressure 
drops  for  five  different  mass  flow  rates  were  measured.  The 
exit  guality  ranged  from  0.54%  to  20.7%.  Inlet  temperatures 
were  maintained  near  43 “C  to  provide  about  5K  of  subcooling. 
A dc  power  of  26V  was  applied  to  the  nichrome  heater  of  13.0 
ohms  resistance,  supplying  52. OW  of  electrical  power  to  the 
test  section.  Heat  loss  was  estimated  by  reading  thermocouple 
outputs  inside  and  outside  of  the  insulator  wrapped  around  the 
tube.  Of  total  heat  input,  12%  of  heat  flux  was  lost  to 
ambient.  The  resulting  net  heat  flux  was  6.63kW/m^. 

Pressure  drop  recordings  for  the  five  cases  are  shown  in 
Figure  4-13.  The  flow  rate  change  during  the  free  fall  was 
within  2%.  On  the  recording  chart,  three  types  of  total 
pressure  drop  were  read  from  the  curve;  a)  the  maximum  value 
appearing  in  the  middle  of  the  free  fall,  b)  the  minimum  value 
appearing  at  the  end  of  the  free  fall,  and  c)  the  average 
value  during  1.25  seconds  of  the  free  fall.  The  above  three 
types  of  pressure  drops  measured  during  the  free  fall 
experiment  are  plotted  in  Figure  4-14  against  the 
thermodynamic  quality  for  the  five  cases.  As  the  quality 
increases,  the  total  pressure  drop  in  microgravity  increases 
in  the  low  quality  region  less  than  1%,  while  they  tend  to 
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Figure  4-13  Total  pressure  drops  measured  in  drop 

tower  experiments.  L=0.5m,  D=5mmx5mm, 
q"w=2 . 93kW/m'^2 . 
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Figure  4-14  Total  pressure  drops  in  a vertical  tube 

under  1 . Og  and  0.02g  conditions. 
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decrease  in  the  higher  quality  region.  This  is  due  to  two 
competing  factors:  one  is  the  decreasing  mass  flux  and  the 
other  is  the  increasing  quality  and  void  fraction. 

Two  measurements  of  pressure  drops  in  normal  gravity  are 
also  shown  in  Figure  4-14.  The  pressure  drops  measured  in 
normal  gravity  indicated  much  larger  total  pressure  drop  due 
to  elevated  static  pressure. 

A homogeneous  model,  in  which  the  velocities  of  the  vapor 
and  liquid  is  considered  equal,  was  used  to  calculate  each 
component  of  the  total  pressure  drop.  The  test  section  is 
simply  divided  into  single  and  two-phase  regions.  For  systems 
with  an  uniform  heat  flux,  the  length  of  the  single-phase 
region  is  given  by 

Lip  ^ (hf-hf  J (4-:?:?) 

L M 

where  Ah  is  the  increase  in  enthalpy  provided  by  the  heater. 
Ah  is  computed  from  the  total  heat  input  Q and  the  mass  flow 
rate  m as  follows 

Ah  = -?  (4-34) 

m 

In  the  present  system  with  a constant  heat  flux,  the  single- 
phase region  becomes  larger  as  the  mass  flow  rate  increases. 
At  the  boundary  of  the  two  regions,  the  thermodynamic  quality 
is  zero.  In  the  two-phase  region,  the  following  are  the 
expressions  for  each  component  of  the  pressure  drop  based  on 
the  homogeneous  model 
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where  fT?  is  the  two-phase  friction  factor  defined  as 


for  turbulent  flow  and 


0.079 


(4-38) 


16 

i?6pp 


(4-39) 


for  laminar  flow.  An  equation  proposed  by  McAdams [39]  was 
adopted  here  for  the  mean  two-phase  viscosity  /iip 


Mtp=x  Mg  + (1-^)  Mf  (4-40) 

In  the  single-phase  region,  the  accelerational  component  is 
negligible.  The  gravitational  and  frictional  components  are 
given  by 


(4-41) 


G^VpLpp 


(4-42) 


In  Eg. (4-36) , the  single-phase  friction  factor  ffo  defined  as 
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f ^ 0.079 


for  turbulent  flow  and 


(4-43) 


16 

Re 


(4-44) 


for  laminar  flow.  In  the  above  calculation,  the  value  of  g 
was  set  to  be  0.196m/s^  for  the  case  of  microgravity. 

Due  to  lack  of  buoyancy  force,  the  slip  ratio  is 
considered  to  be  close  to  unity.  Under  such  circumstances, 
calculation  by  the  homogeneous  two-phase  flow  model  is 
expected  to  provide  a reasonable  estimate  for  pressure  drops 
in  microgravity.  The  predicted  total  pressure  drops  for 
corresponding  cases  are  plotted  in  Figure  4-14  together  with 
experimental  data.  The  solid  line  indicates  the  predicted 
total  pressure  drops.  The  predicted  values  agrees  with 
experimental  results  in  the  high  quality  region  larger  than 
2%,  while  the  prediction  tends  to  underestimate  in  the  low 
quality  region.  In  the  low  quality  region  where  the  flow 
pattern  is  mainly  bubbly  flow,  the  increase  in  void  fraction 
is  significantly  large,  as  shown  in  Figure  4-3 (g) . This 
result  indicates  that  the  void  fraction  and  pressure  drop  in 
the  subcooled  region  are  no  longer  negligible. 

To  investigate  further  the  influence  of  gravity  on 
pressure  drop,  each  component  in  the  two-phase  region  for 
normal  gravity  and  microgravity  conditions  was  separated  and 
compared  to  each  other  in  Figure  4-15 (a)  and  (b) . In 
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(a)  l.Og 


Quality,  X 


(b)  0.02g 


Figure  4-15  Influence  of  gravity  on  pressure  drop 

components  in  a vertical  tube. 
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microgravity,  the  frictional  component  is  dominant  in  both  the 
two-phase  and  single-phase  regions.  Since  AP2f  and  AP2a 
include  the  square  root  of  the  mass  flux,  both  components  tend 
to  decrease  in  the  high  quality  region  where  G is  reduced  to 
attain  higher  qualities  for  the  present  study.  Single-phase 
pressure  drop  linearly  decreases  as  L^p  becomes  smaller.  On 
the  other  hand,  in  normal  gravity,  the  static  pressure  in 
the  single-phase  region  is  dominant  for  the  present  range  of 
quality,  and  the  increase  in  frictional  and  accelerational 
components  in  the  two-phase  region  hardly  contributes  to  the 
total  pressure  drop. 

As  indicated  previously,  the  frictional  component  of  the 
pressure  drop  is  of  paramount  importance  to  the  heat  transfer 
problem  in  microgravity.  The  wall  shear  stress  for  a tube  is 
correlated  with  the  gradient  of  frictional  pressure  drop  by 


(4-45) 


With  the  measurement  procedure  used  in  the  present  experiment, 
(AP/Az)2f  cannot  be  measured  directly.  Assuming  that  the  void 
fraction  profile  is  known,  the  calculated  values  of  APp,  AP2g, 
and  DP2a  were  subtracted  from  the  measured  total  pressure  drop 
to  obtain  APpf.  Then,  AP2f  was  divided  by  L2P  to  yield  the 
average  frictional  pressure  drop  gradient.  The  same  procedure 
was  applied  to  the  data  obtained  under  normal  gravity 
conditions.  The  result  is  shown  in  Figure  4-16,  with  the 
predictions.  The  maximum,  average,  and  minimum  values  of  total 
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Figure  4-16 


Frictional  pressure  drop  gradient 
in  a vertical  tube  under  l.Og  and 
0.02g  conditions. 
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pressure  drop  measured  during  the  microgravity  period  are 
plotted  separately.  (AP/Az)2f  in  the  low  quality  region  under 
microgravity  conditions  became  significantly  larger  than  those 
of  normal  gravity  conditions.  A conceivable  reason  for  the 
increased  (AP/Az)2f  in  the  low  quality  region  is,  again,  the 
formation  of  significantly  large  bubbles  in  the  liquid  flow. 
Taking  into  account  the  correlation  provided  by  Eg.  (4-45)  , the 
wall  shear  stress  in  microgravity  in  the  low-quality  region 
seems  to  become  larger  than  that  in  normal  gravity  under  the 
conditions  of  present  experiment.  This  result  indicates  that 
the  pressure  drop  in  the  subcooled  region  is  no  longer 
negligible  under  microgravity  conditions. 

Separated  flow  models  can  provide  more  accurate 
predictions  if  the  void  fraction  profile  along  the  flow  is 
given.  The  axial  evolution  of  void  fraction  was  not  available 
in  the  present  experiment,  although  the  volume  fraction  near 
the  test  section  exit  was  determined  from  the  photographs. 
The  existing  void  fraction  correlations  are  flow-pattern 
sensitive  and  tend  to  deviate  from  each  other  depending  on  the 
heat  flux,  mass  flow  rate,  and  geometry  of  heaters,  especially 
in  the  low-quality  region.  The  test  conditions  of  the  present 
study,  characterized  by  low  flow  rate  and  low  heat  flux,  are 
very  unique  and  difficult  to  find  existing  correlations 
applicable  to  the  present  experiments.  Because  of  the  above 
reasons,  the  separated  flow  model  was  not  used  in  the  present 
study . 


CHAPTER  5 

GRAVITY  EFFECTS  ON  SUBCOOLED  BOILING  HEAT  TRANSFER 

5.1  Flight  Experiment  Results 

On  September  16,  1992,  two  racks  containing  the  flow 

boiling  equipment  and  data  acquisition  system  were  aboard 
NASA's  LeRC  microgravity  Lear jet  and  subjected  to  a series  of 
lO'^g  microgravity  tests.  Each  microgravity  trajectory  lasted 
for  approximately  18  seconds. 

Test  conditions  of  the  Lear  jet  experiments  are  summarized 
in  Table  5-1.  The  ground  test  results  with  the  preheater 
turned  on  showed  that  approximately  three  minutes  were 
required  for  start-up.  The  thermal  equilibrium  qualities  at 
the  exit  of  the  test  section  for  these  tests  ranged  from 
-0.04%  to  -12.1%  and  maintained  within  the  subcooled  region. 
The  mass  fluxes  were  changed  from  165  to  384kg/mVsec* 
Without  using  the  preheater,  the  inlet  temperatures  were  kept 
between  15.2°C  to  20.3°C.  A constant  dc  current  of  14. 7A  was 
supplied  to  the  main  heater.  The  electrical  power  input  was 
411W.  The  heat  flux  calculated  from  the  current  and  the 
resistance  measured  was  24.3kW/m^  at  the  operating 
temperatures.  Heat  loss  to  ambient  was  measured  by  reading 
temperatures  inside  and  outside  the  insulator  around  the  test 
section  and  found  to  be  0.4kW/m^. 
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Table  5-1  Learjet  flight  experiment  conditions 


129 


(0 

CO 

r- 

MO 

CN 

00 

MO 

r- 

Oi  — 

OO  o 

o 

tH  00 

00  MO 

rH  O 

rH  CN 

cn  CO 

tH  rH 

iH  CN 

S-H 

O CN 

CN  OO 

o o 

O rH 

CN  O 

O O 

O CN 

CN  rH 

o o 

cn 

o • 

o • 

o • 

O • 

O • 

O • 

o • 

O • 

o • 

cu  a 

• rH 

• CN 

• o 

• tH 

• CO 

• O 

• rH 

• CO 

• o 

<a 

O 

O 

O 

O 

O 

O 

O 

O 

O 

(0 

— ^ 

^ 

^ 

V 

^ 

CM 

CO 

o 

o 

00 

r- 

tH 

o 

s ^ 

00  VD 

lO  00 

^ cr» 

cn 

CO  cn 

00  CN 

VX)  MO 

00  rH 

o cn 

Ct.H 

1 — 1 vx 

<N  MO 

cn  o 

rH  CO 

^ rH 

cn  r- 

iH  CO 

CN  O 

rH  ^ 

e CO 

rH  • 

o • 

tH  • 

tH  • 

o • 

rH  • 

rH  • 

rH  • 

o a 

• ^ 

• MO 

• t-H 

• ^ 

• 00 

• iH 

• ^ 

• KD 

• CO 

CM 

O 

O 

O 

O 

O 

O 

O 

O 

O 

(rt 

Cm 

m 

VD 

r- 

MO 

CN 

cn 

CO 

CO 

S ^ 

O lO 

MO  CO 

CO  o 

cn  CO 

CN  cn 

OO  cn 

OO  CN 

^ cn 

O VX> 

tH  'sr 

o r- 

cn  o 

O CN 

CN  iH 

cn  v£) 

O rH 

O OO 

rH 

a CO 

rH  • 

tH  • 

o • 

t-H  • 

rH  • 

o • 

rH  • 

iH  • 

tH  • 

° CM 

• CO 

• CN 

• t-H 

• CO 

• MO 

• iH 

• CO 

• CN 

• CO 

CM 

O 

O 

O 

o 

O 

O 

O 

O — 

O 

CO 

\ 

cvj 

00 

(Tk 

cn 

VD 

00 

CN 

MO 

yv 

CO 

O 

cn 

CN 

r- 

<D 

OO 

tr>S 

rH 

o 

CN 

f-H 

rH 

CN 

CO 

CO 

CO 

CM 

iT)  . 

• 

O • 

00  • 

'X)  • 

• 

CO  • 

O • 

^ • 

o 

VD  O 

r-  o 

rH  O 

cn  o 

o 

CN  o 

r-  o 

'X)  O 

00  o 

O 

I— 1 

CO 

(N 

1 — 1 

1 — 1 

CN 

CO 

CO 

CO 

O 

<D 

cn 

VD 

CN 

MO 

o 

CN 

lT) 

oO 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

KO 

r* 

tH 

tH 

VD 

VD 

<X) 

H 

CO 

CO 

CO 

CO 

CN 

CN 

CN 

c 

•H 

cr> 

'X) 

<D 

cn 

rH 

lO 

00 

«o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

H o 

rH 

CO 

MO 

'X) 

CN 

rH 

r- 

CN 

<1  — ' 

CO 

CO 

CN 

CO 

CO 

CO 

CO 

CO 

n 

CO 

CO 

00 

CN 

CO 

CN 

CN 

CO 

■hO 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

o 

cn 

m 

IT) 

MO 

r- 

r- 

r- 

Eh 

CN 

CN 

rH 

tH 

rH 

t 

tH 

iH 

tH 

•H 

T3 

O — 

1 

1 

1 

1 

1 

1 

1 

1 

1 

o r- 

1 — 1 kO 

CN 

rH  cn 

vx)  r- 

VD  ^ 

MO  r- 

CO  CO 
CN  iH 

•H  O 

M 0 

O CJ^ 

O CO 

MO  00 

00 

MO  cn 

O CO 

o cn 

MO  r- 

0 CO 

o o 

CN  CN 

00  cn 

o o 

CO  CO 

MO  MO 

r-  r- 

cn  o 

04 

CN  CN 

CN  CN 

CN  CN 

tH  tH 

CN  CN 

CN  CN 

rH  tH 

tH  tH 

rH  CN 

No. 

rH 

CN 

O 

tH 

CN 

O 

iH 

CN 

o 

LT) 

lD 

MO 

o 

o 

O 

MO 

MO 

MO 

0 

rH 

tH 

tH 

CN 

CN 

CN 

CO 

CO 

CO 

CO 

o 

o 

o 

O 

o 

O 

O 

o 

o 

(0 

tH 

I— 1 

rH 

1—1 

rH 

rH 

tH 

rH 

t — 1 

U 

Cp 

pp 

Cp 

Pp 

Pp 

Pp 

Pp 

Pp 

Pp 

• 

TJ 

Q) 

s: 

(0 

•H 

I— H 

fO 

-P 

(0 

(1) 

(0 

ft 

0) 


•iH 

4-1 

• o 

O M 
pc:  a 
(D 

PI  0) 
M 

< 

CO  (0 
< (0 
Z 0) 

p 

-M  a 

ft 

ft 

ro 

(T>  4-> 
\ C 
CO  ft 
«H  4-J 
\ (0 
C 
O 
C U 
O 

C 

TD  (D 
0 x: 

-p  > 

o 

D TD 
T5  O 
C -H 
O P 
CJ  0 

a 

0 

p 0 

0 x: 
> p> 

CO  0 
-p  p) 
C 0 
0 O 
S *H 
•P  XJ 

P c 

0 *rH 

a 

X CO 
0 0 
B 

-p  *p 
x:  4-> 

cn 

■H  0 
rH  X 

Eh 


+ -)< 


130 


Although  two  flights  were  scheduled,  only  the  first 
flight  was  completed  because  of  a malfunction  of  the  airplane 
control  system.  Five  parabolas  were  flown  during  the  single 
flight.  The  qualities  and  flow  patterns  at  the  exit  of  the 
test  section  are  summarized  in  Table  5-2.  During  the  five 
parabolas,  two  data  sets  in  virtual  memory  obtained  at  the 
first  and  fourth  trajectories  were  erased  because  the  computer 
power  was  lost  just  before  entering  the  low-gravity  maneuver. 
The  probable  cause  of  the  power  loss  is  an  unsecured 
electrical  contact  of  the  power  supply  inside  the  computer, 
which  failed  during  the  transition  from  2.4g  to  Og. 

Figure  5-1  (a)  through  (c)  show  the  acceleration  histories 
in  x-(fore-aft) , y- (right-left)  , and  z-(up-down)  directions 
recorded  in  a built-in  accelerometer.  Here,  the  z-axis  is 
parallel  to  the  flow  direction  in  the  vertical  test  section. 
It  was  attempted  to  null  acceleration  in  the  z-direction.  The 
first  and  second  curves  in  Figure  5-1  indicate  the 
accelerations  in  the  z-directions  recorded  in  coarse (l/20g  per 
scale)  and  fine(l/200g  per  scale)  modes,  respectively. 
Although  the  z-axis  acceleration  fluctuated  during  the  initial 
period  of  each  trajectory,  it  damped  out  within  ±l/100g  as 
time  progressed.  Acceleration  along  the  x-axis  was  almost 
nullified,  while  that  of  the  y-axis  was  maintained  slightly 
positive  to  follow  the  parabolic  trajectory. 

Bulk  temperature  profiles,  flow  rates,  pressure  drops, 
and  wall  temperature  profiles  for  Cases  F1015,  F1020,  and 


Table  5-2  Summary  of  experimental  conditions  for  Lear jet  flight  experiments  . 
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Figure  5-1  Acceleration  histories  recorded  by  an  on-board 

accelometer  during  the  Learjet  flight  experime 
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Figure  5-1  Acceleration  histories  recorded  by  an  on-board 

accelometer  during  the  Learjet  flight  experime 
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Figure  5-1  Acceleration  histories  recorded  by  an  on-board 

accelometer  during  the  Learjet  flight  experime 
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F1035  are  shown  in  Figure  5-2  through  5-4,  respectively.  Flow 
rates  of  165,  198,  and  373kg/m^/sec  were  tested  while 
maintaining  a constant  heat  flux  of  23.9kW/m^. 

Bulk  temperatures  at  the  exit  of  the  test  section  for 
Case  F1015  are  shown  in  Figure  5-2 (a).  The  profile  indicates 
that  the  system  requires  about  200  seconds  for  start-up  for 
the  same  inlet  temperature.  The  exit  bulk  temperature  rose 
rapidly  at  t=212sec  because  the  flow  stopped  during  the 
pushover  before  entering  the  microgravity  trajectory.  Upon 
entering  the  microgravity  trajectory,  the  flow  started  again 
and  the  bulk  temperature  reached  a steady  state  after  a few 
fluctuations.  The  acceleration  level  estimated  from  the 
static  pressure  was  2.4g  during  this  period.  The  observed 
flow  pattern  transition  was  slug  to  annular,  indicating  that 
the  liquid  was  deficient  in  the  upper  portion  of  the  test 
section.  Expansion  of  the  vapor  phase  allowed  the  liquid  to 
stay  in  the  lower  portion  of  the  test  section  during  the  high 
acceleration  period.  The  decrease  in  flow  rate  became  less 
significant  during  the  succeeding  subcooled  boiling  tests 
where  liquid  phase  was  continuous  along  the  test  section.  The 
above  incident  implies  that  experiments  for  higher  thermal 
equilibrium  qualities  may  encounter  flow-rate  instabilities  in 
vertical  flow  systems  during  the  high  acceleration  period 
unless  sufficient  driving  pressure  is  supplied.  To  avoid  the 
flow  interruption,  the  boiling  condition  at  the  test  section 
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(a)  Bulk  temperatures 

Figure  5-2  Variation  of  flow  properties  for  Case  F1015  during 

the  flight  experiment. 
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(b)  Flow  rate 

Figure  5-2  Variation  of  flow  properties  for  F1015  during 

the  flight  experiment. 
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(c)  Pressure  drop 

Figure  5-2  Variation  of  flov/  properties  for  Case  F1015  during 

the  flight  experirr:ent , 


Wall  temperature,  (K) 
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(d)  Wall  temperatures 


Figure  5-2 


Variation  of  flow  properties  for  Case  F1015 
the  flight  experiment. 
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Figure  5-2  Variation  of  flow  properties  for  Case  F1015  during 

the  Lear jet  flight  experiment. 
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(f)  Wall  superheats 

Figure  5-2  Variation  of  flow  properties  for  Case  F1015 

during  the  flight  experiment. 
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(g)  Subccolings 

Figure  5~2  Variation  of  flow  properties  for  Case  F1015 

during  the  flight  experiment. 
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(a)  Bulk  temperatures 

Figure  5-3  Variation  of  flow  properties  for  Case  F1020  during 

the  flight  experiment. 
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(b)  B’iow  rates 

Figure  5-3  Variation  of  flow  properties  for  Case  F1020  during 

the  flight  experiment. 
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(c)  Pressure  drops 

Figure  5-3  Variation  of  flow  properties  for  Case  F1020  during 

the  flight  experiment. 


Wall  temperature,  (K) 


146 


195.5sec(1.0g) 

210.55ec(2.4c) 

■ ■■A.' 

233.3sec(0.01g) 

— — 

Tf(lg) 

•A  • 

Tf(0.01g) 


(d)  Wall  temperatures 


Figure  5-3  Variation  of  flow  properties  for  Case  F1020  during 

the  flight  experiment. 
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Figure  5-3  Variation  of  flow  properties  for  Case  F1020  during 

the  Lear jet  flight  experiments. 
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(f)  Wall  superheats 


Figure  5-3  Variation  of  flow  properties 

during  the  flight  experiment. 
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(g)  Subcoolings 


Figure  5-3  Variation  of  flow  properties  for  Case  F1020 

during  the  flight  experiment. 
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Figure  5-4  Variation  of  flow  properties  for  Case  F1035  during 

the  flight  experiment. 
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(b)  Flow  rates 

Figure  5-4  Variation  of  flow  properties  for  Case  F1035  during 

the  flight  experiment. 
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(c)  Pressure  drops 


Figure  5-4  Variation  of  flow  properties  for  Case  F1035  during 

the  flight  experiment. 
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Figure  5-4  Variation  of  flow  properties  for  Case  F1035  during 

the  Lear jet  flight  experiments. 
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(f)  Wall  superheats 

Figure  5-4  Variation  of  flow  properties  for  Case  F1035 

during  the  flight  experiment. 
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Variation  of  flow  properties  for  Case  F1035 
during  the  flight  experiment. 


157 


exit  for  the  present  system  was  confined  to  subcooled 
qualities. 

The  variation  of  gauge  pressures  located  at  the  inlet  and 
outlet  of  the  test  section  during  the  experiments  are  shown  in 
Figure  5-2 (b) . The  ambient  pressure  in  the  cabin  on  that  day 
was  0.85atm.  Since  all  the  test  conditions  at  the  exit  fall 
in  the  subcooled  region,  the  acceleration  pressure  drop  was 
small  compared  with  the  static  pressure  drop.  In  normal 
gravity,  the  gravitational  pressure  drop  was  dominant  compared 
to  other  components.  According  to  error  propagation  formula, 
the  uncertainty  range  of  the  differential  pressure  was  ±38.5Pa 
based  on  the  accuracy  of  ±27.5Pa  for  a single  pressure 
transducer.  Since  the  built-in  accelerometer  was  independent 
of  our  data  acquisition  system,  matching  between  the  events 
and  elapsed  time  was  made  by  reading  the  pressure  data. 

As  shown  in  Figure  5-2 (c),  the  flow  rates  increased 
during  the  microgravity  trajectories  with  the  decrease  in 
static  pressures.  Since  the  acceleration  in  microgravity 
became  about  1/100  of  normal  gravity,  the  static  pressure  of 
Freon  113  in  a 0.5m  vertical  tube  changed  from  7585Pa  to 
75.9Pa  at  room  temperature.  Correspondingly  the  flow  rates 
increased  by  27.1%,  14.6%,  and  2.9%  for  Cases  F1015,  F1020, 
and  F1035,  respectively.  This  change  was  undesirable  because 
the  objective  of  the  present  study  was  to  investigate  the 
gravity  effects  on  heat  transfer  with  other  conditions 
unchanged.  Especially  for  Case  F1015,  the  flow  completely 
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stopped  at  t=222.0sec.  As  can  be  seen  from  the  pressure  data, 
this  coincides  with  the  time  when  the  highest  acceleration  was 
achieved.  For  Case  F1035,  however,  the  flow  rate  change  as 
well  as  the  pressure  change  was  small  enough  to  allow  for 
comparison  of  microgravity  data  with  normal  gravity  data.  On 
the  analysis  of  Cases  F1015  and  F1020,  changes  in  the  flow 
rate  and  pressure  were  considered. 

To  successfully  operate  the  present  system,  the  liquid 
phase  must  be  continuous  to  avoid  the  flow  interruption,  or 
the  driving  pressure  must  be  large  enough  to  exceed  the 
pressure  losses  in  the  loop.  The  calculation  of  pressure 
losses  in  the  loop  is  shown  in  Appendix  A. 

Wall  temperature  and  estimated  bulk  temperature 
distributions  under  normal  and  micro  gravity  conditions  are 
shown  in  Figure  5-2 (d).  The  bulk  temperatures  at  the  inlet 
and  exit  of  the  test  section  were  measured  by  l.SmmO.D.  E-type 
thermocouples  immersed  in  the  flow.  In  the  two-phase  region, 
the  axial  bulk  temperatures  were  estimated  by  Ahmad's  equation 
without  using  the  condensation  parameter [40] . In  the  single- 
phase region,  bulk  temperatures  were  estimated  simply  by  the 
following  heat  balance 


T,(z)  = 


\ 


A I r. 


(5-1) 


In  the  subcooled  boiling  region,  the  dimensionless  subcooling 
is  expressed  as 
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(T,-r,) 


(5-2) 


where  Cj  is  the  heat  transfer  parameter  and  z*  is  the 
dimensionless  distance  (=z/L2p)  . Ci  is  defined  as 

(5-3) 

where  hfo  is  the  single-phase  forced  convective  heat  transfer 
coefficient,  Lpp  is  the  significant  boiling  length  given  by 

(5-4) 

(J.  D 

ATpn  and  AT^  denote  the  subcoolings  at  the  inlet  and  the  bubble 
detachment  point,  respectively.  An  array  of  17  fine-gauge 
thermocouples  were  mounted  on  the  0.5m  long,  0.71mm  wall 
thickness  copper  tube  at  equal  intervals.  The  outputs  from 
the  thermocouples  were  calibrated  by  9th  order  polynomials. 
Accuracy  of  the  temperature  reading  was  ±0.5°C.  Wall 
temperatures  at  both  ends  were  lower  than  others  because  of 
pressure  taps  and  tube  fittings  diffusing  heat  from  the  copper 
tube.  Therefore,  these  2 data  points  at  each  end  were 
discarded. 

For  Case  F1015,  the  wall  temperatures  decreased 
significantly  along  the  test  section  in  microgravity, 
indicating  a large  increase  in  heat  transfer  rate.  Although 
the  liguid  flow  once  stopped  and  re-started,  it  was  stabilized 
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at  the  end  of  the  microgravity  trajectory  of  18  seconds 
(Figure  5-2 (b)).  The  flow  pattern  changed  from  slug  flow  to 
annular  flow  in  microgravity.  Our  previous  results  in  the 
drop  tower  experiments  also  showed  a similar  transition  from 
slug  to  annular  flow  in  microgravity. 

As  shown  in  Figure  5-2  (d)  , a peak  of  the  wall  temperature 
profile  near  the  test  section  inlet  disappeared  under 
microgravity  conditions.  From  the  data  set  obtained  in  this 
experiment,  it  is  difficult  to  correlate  the  wall  temperature 
profile  to  the  void  fraction  profile  in  the  flow  direction. 
However,  a similar  subcooled  experiment  using  Freon  113  was 
conducted  by  Hino  and  Ueda[41].  They  measured  wall 
temperature  profiles  and  photographed  the  void  fraction 
distributions  simultaneously.  It  was  found  that  the  peak  in 
the  wall  temperature  curve  corresponds  to  the  onset  of 
nucleate  boiling  point  (ONE) . The  profiles  in  microgravity 
indicate  that  the  wall  temperatures  at  the  ONE  decreases 
significantly. 

The  wall  temperatures  in  the  subcooled  region  near  the 
test  section  outlet  under  microgravity  conditions  decreased  by 
1-2 °C.  Since  the  system  pressures  have  decreased  during  the 
transition  from  2.4g  to  O.Olg,  the  change  of  saturation 
temperature  was  taken  into  account  to  investigate  the  gravity 
effect  on  heat  transfer.  The  overall  wall  temperature 
transition  along  the  test  section  is  shown  in  Figure  5-2 (e). 
The  wall  superheat  and  liquid  subcooling  are  separately  shown 
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in  Figure  5-2 (f)  and  5-2 (g),  respectively.  The  increase  in 
wall  superheats  near  the  test  section  outlet  indicates  that 
the  bulk  convection  heat  transfer  was  decreased  in 
microgravity . 

Similar  tendencies  were  observed  in  Figure  5-3  and  5-4 
for  Cases  F1020  and  F1035,  respectively.  Since  the  exit 
condition  was  highly  subcooled  for  Case  F1035,  the  static 
pressure  change  was  minimal  among  the  three  data  sets.  Case 
F1035,  therefore,  exhibits  the  gravity  effects  with  a minimum 
interference  of  pressure  change. 

Characteristic  boiling  regions.  For  highly  subcooled 
liquid  flowing  in  a heated  channel,  the  heat  transfer  process 
evolves  from  forced  convective  single-phase  flow  to  subcooled 
flow  boiling,  then  to  saturated  flow  boiling.  The  subcooled 
boiling  region  is  further  divided  into  2 regions  at  the  bubble 
detachment  point:  the  highly  subcooled  region  with  low  void 
fraction  where  bubbles  are  attached  or  slide  along  the  wall, 
and  the  slightly  subcooled  region  with  large  void  fraction 
where  the  detached  bubbles  flow  with  the  liquid. 

As  boundaries  of  the  above  regions,  two  characteristic 
points  are  defined.  One  is  the  onset  of  nucleate  boiling 
(ONB)  point  at  which  the  first  bubble  formation  is  observed, 
and  the  other  is  the  bubble  detachment  point  at  which  bubbles 
break  off  the  wall  and  enhance  the  stirring  motion  in  the 
liquid.  Figure  5-5  shows  the  typical  flow  pattern  and  the 
variation  of  surface  and  liquid  temperatures.  The  range  of 
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the  present  study  is  confined  to  the  single-phase  region, 
highly  subcooled  region,  and  a part  of  slightly  subcooled 
region.  In  microgravity , the  heat  transfer  mechanisms  are 
considered  to  be  different  from  those  in  normal -gravity.  The 
existing  models  to  predict  the  boundaries  of  the  above  regions 
cannot  be  applied  to  microgravity  heat  transfer  without 
verification.  In  the  following  section,  the  gravity  effects 
on  the  single-phase  and  subcooled  heat  transfer  mechanisms 
will  be  discussed,  followed  by  discussions  of  gravity  effects 
on  the  onset  of  nucleate  boiling. 

5.2  Forced  Convective  Sinale-Phase  Flow 


A simple  heat  balance  in  this  region  yields  the  liquid 
temperature  profile.  Eg. (5-1)  was  used  to  estimate  the  liquid 
temperature  within  the  single-phase  region.  The  wall 
temperature  is  given  by 


r„(z)  =Tf(z) 


(5-5) 


where  hfo  is  the  single-phase  heat  transfer  coefficient.  The 

1 

minimum  and  maximum  Reynolds  numbers  for  the  present  tests 
were  4145  and  9656,  respectively,  which  fall  in  the 
intermediate  range  between  laminar  and  turbulent  flow.  The 
decrease  in  system  pressure  during  the  microgravity  maneuver 
induced  an  increase  in  the  mass  flow  rate.  Although  the  inlet 
temperature  was  maintained  constant  throughout  the  experiment, 
this  increase  in  mass  flow  rate  altered  bulk  temperature 
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profiles.  Among  three  cases,  F1035  showed  a minimum  decrease 
in  bulk  temperature  of  0.6K  at  the  test  section  exit,  while  a 
maximum  decrease  of  I.IK  was  observed  for  F1020. 

5.2.2  Contribution  of  Natural  Convection 

It  is  assumed  that  the  effects  of  gravity  change  is 
negligible  on  single-phase  turbulent  flow  since  no  phase 
change  is  involved.  However,  natural  convection  heat 
transfer,  which  is  driven  by  the  buoyancy  force,  is 
significantly  affected  by  the  gravity  change.  In  general,  the 
Nusselt  number  is  expressed  as  a function  of  Reynolds  number 
Re,  Prandtl  number  Pr,  and  Grashof  number  Gr  defined  as 

Gr-.£l££!£i^  (5-6) 

During  the  microgravity  trajectories,  the  acceleration  level 
reached  ±0.01.  The  corresponding  Gr  is  approximately  1/100  of 
Gr  under  normal  gravity.  In  laminar  flow  correlations,  Nu  is 
proportional  to  g°’^  and  a change  in  g by  1/100  reduces  Nu  by 
37%.  In  a vertical  flow  system,  if  the  flow  rate  is  low  and 
the  liquid  Reynolds  number  is  less  than  2000,  the  natural 
convection  is  significantly  affected  by  the  gravity  level. 

For  the  relatively  low  Reynolds  numbers  in  this 
experiment,  the  contribution  of  natural  convection  heat 
transfer  to  that  of  laminar  and  turbulent  heat  transfer  was 
estimated  by  Jackson's  correlation [55]  derived  for 
intermediate  Reynolds  number  ranges  near  the  entrance  region. 
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This  correlation  was  chosen  particularly  because  the 
contributions  of  laminar,  turbulent,  and  natural  convection 
for  forced  convective  flow  in  a pipe  was  distinguishable 

Nu  = 1.74  [ (0.785RePrR +Re^Pr^'^xio-^)^ 

Z ( V ) 

+(GrPrf^*  xiO'^]^''^ 

The  first,  second,  and  third  term  on  the  RHS  of  Eq.  (5-7) 
indicate  the  contributions  of  laminar,  turbulent,  and  natural 
convection,  respectively.  The  effects  of  gravity  change  in 
the  single-phase  flow  on  Nu  were  negligible  since  the  relative 
magnitude  of  the  natural  convection  accounts  for  less  than  1%. 
In  fact,  the  wall  temperature  profiles  in  this  region  showed 
no  change  during  the  gravity  transition,  as  shown  in  Figure  5- 
3(d)  and  5-4 (d).  In  Case  F1015,  the  present  analysis  was  not 
applicable  since  the  flow  pattern  changed  from  slug  to  annular 
flow. 

5.2.3  Single-Phase  Heat  Transfer  Coefficient 

Due  to  the  unique  test  conditions  and  configuration  of 
the  test  section,  none  of  the  existing  correlations 
successfully  predicts  the  single-phase  heat  transfer  rate  for 
the  present  study.  However,  the  single-phase  heat  transfer 
still  accounts  for  substantial  fractions  in  the  subcooled 
region.  Successful  prediction  of  the  gravity  effects  on 
boiling  heat  transfer  depends  on  the  accurate  prediction  of 
single-phase  flow  heat  transfer. 
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A new  correlation  specifically  applicable  to  our  study 
was  developed  based  on  the  turbulent  flow  heat  transfer  model 
proposed  by  Sieder-Tate.  Nuq  for  fully  developed  turbulent 
flow  is  given  by 


Nu=0.023 

' GD 

0.8 

^pl 

1/3 

Mf  J 

J 

(5-8) 


Since  the  inlet  of  the  vertical  tube  is  located  0.14m  above 
the  bend,  the  flow  was  not  completely  developed  in  the  inlet 
region.  Normally,  turbulence  caused  by  the  edge  of  the  tube 
increases  the  heat  transfer  rate.  Eg. (5-8)  alone 
underestimates  the  heat  transfer  rate  in  this  region. 
According  to  [55],  the  local  Nu(z)  in  the  entrance  region 
should  be  multiplied  by  the  following  correction  factor 


Nu(z)  =1. 57  Nu.  I — 


0.054 


(5-9) 


The  validity  of  the  above  equation  is  Re  > 8,000  and  (L/D)  < 
60.  The  factor  1.57  was  modified  to  fit  the  range  of 
appropriate  Reynolds  numbers.  Substitution  of  Eg (5-8)  into 
Eg. (5-9)  yields 


Nu(z) 


/ \0.14 


0.054 


(5-10) 


where  Ci  and  C2  are  constants  to  be  determined  experimentally 
through  least  square  fitting.  Three  data  sets  obtained  in  the 
flight  experiment  were  used  to  correlate  Re  with  Nu  in  the 
transitional  range  between  laminar  and  turbulent  flow.  The 


Nusselt  Number,  Nu 
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Reynolds  Number 


Single-phase  heat  transfer  in  a 
vertical  tube  in  Learjet  flight 
experiments . 


Figure  5-6 
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result  of  curve  fitting  is  shown  in  Figure  5-6.  The  pressure 
in  the  Learjet  cabin  was  0.0862MPa,  while  that  of  ground  tests 
was  O.lOlMPa.  The  coefficients  for  the  new  correlation  are 
Ci=2.192,  and 

€2=0.330. 

The  validity  of  the  above  correlation  is  confined  to  4145  < Re 
< 9656  and  q„”=23 . 9kW/m^ . Nu  for  the  experiments  were 

calculated  as  follows 


^1^(2)  exp 


n 

D 

kAT-T,) 


(5-11) 


The  wall  temperature  T„  was  measured  on  the  outer  surface  of 
the  tube,  while  the  bulk  temperatures  between  the  inlet  and 
outlet  of  the  test  section  were  estimated  by  Eq. (5-1) . 

Two  other  correlations  for  forced  convective  heat 
transfer  were  also  plotted  against  Re.  One  is  Jackson's 
correlation  (Eq.(5-8)),  and  the  other  is  that  of  Warmeatlas 
[42]  derived  for  the  transition  between  laminar  and  turbulent 
flow 


Nu  = 


(a/8) (Re -1000)  Pr 
1+12.7  (a/8)^^^ 


2/3 


r,  \0.11 

Pr 


Pr 


xO.86 


D 
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-0.16 


wy 


/ 


a=  (1.821ogj^gRe-1.64) 


(5-12) 


Eq. (5-8)  tends  to  overestimate  the  heat  transfer  rate  since 
this  relationship  was  developed  for  smaller  Pr.  The  value  of 
Pr  for  Freon  113  was  on  the  order  of  7.0.  Eq.(5-12),  on  the 
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other  hand,  underestimates  the  heat  transfer  rate  for  small 
Re.  This  is  probably  due  to  the  relatively  small  flow  rate 
used  in  the  present  experiment.  Eq. (5-10)  follows  the 
variation  of  Nu  over  the  entire  range  of  Re  in  the  present 
experiments . 

It  should  be  noted  that  this  correlation  was  derived  only 
for  the  present  study  from  a limited  number  of  data.  Thus, 
this  equation  may  not  be  applicable  elsewhere. 

5.2.4  Effects  of  Pressure  on  Nusselt  Number 

Vertical  flow  lines  including  the  boiling  test  section 
were  subject  to  the  transition  of  gravity  varying  from  l.Og 
through  2.4g  and  to  O.Olg  during  the  trajectory.  As  shown  in 
Table  5-1,  the  largest  pressure  decrease  due  to  gravity  change 
was  3.5psi  experienced  in  Case  F1015.  This  pressure  change 
induced  an  increase  in  flow  rate  by  25.8%.  As  a result,  Nu 
was  increased  by  7.7%.  Using  Eq. (5-10) , the  ratio  of  Nu  under 
Ig  and  O.Olg  conditions  is  given  by 

The  increase  in  Nu  in  the  single-phase  region  for  F1015, 
F1020,  and  F1035  were  +4.6%,  +5.0%,  and  +1.0%,  respectively. 
This  result  indicates  that  the  change  of  heat  transfer  rate 
observed  for  Case  F1035  under  microgravity  conditions  reflects 
the  gravity  effects  with  least  effect  of  pressure  change.  The 


Pr 


Og 


Pr 


1 a 


(5-13) 
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increase  in  flow  rate  for  Case  F1035  in  microgravity  was  3.0%, 
the  smallest  of  the  three  cases. 

If  one  hypothetically  maintains  the  same  flow  rate  with 
the  above  pressure  change,  Nu  would  decrease  by  1.4%.  In 
general,  a decrease  in  pressure  reduces  the  heat  transfer 
rate.  However,  this  decrease  in  Nu  was  exceeded  by  the 
effects  of  the  increase  in  the  flow  rate. 

5.3  Onset  of  Nucleate  Boiling 

In  the  beginning  of  the  subcooled  region,  the  wall 
temperatures  in  microgravity  significantly  decreased.  In  this 
section,  theories  of  the  onset  of  heterogeneous  nucleation  are 
introduced.  Then  the  effect  of  pressure  change  on  the  ONE  is 
examined  to  identify  the  gravity  effects. 

5.3.1  Superheat  at  ONE 

Existing  analytical  solutions  to  predict  the  heat  flux 
and  the  wall  superheat  at  the  onset  of  nucleate  boiling  were 
reviewed  by  Dougall  and  Lippert[56].  Among  several  theories, 
those  proposed  by  Eergles  and  Rohsenow,  Hsu,  and  Frost  and 
Dzakowic  were  applied  to  the  present  study. 

Hsu  postulated  that  the  criterion  for  nucleation  is  that 
the  temperature  of  the  liquid  at  the  top  of  the  bubble  should 
exceed  the  equilibrium  temperature  of  the  vapor  in  the 
bubble[57].  If  a full  spectrum  of  cavity  size  is  assumed  to 
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exist  on  the  heating  surface,  the  criterion  is  given  by  the 
tangent  of  liquid  temperature  profile 


T,(z) 


\ 


^ y 


(5-14) 


and  the  vapor  equilibrium  temperature  profile 

^g  = r3(l  + -pr^)  (5-15) 

■^crit 

where  A is  a combination  of  fluid  properties.  As  values  of  A, 
the  following  values  were  employed 
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{Frost-Dzakowic [61] ) 


At  the  tangent,  the  following  relationships  are  satisfied 


Tf  = T and  , ^ = — i-^f  at  y = nr  ,,  (5-17) 

I g (jy 

where  n*  r^^it  indicates  the  distance  from  the  wall  to  the 
tangent  of  two  profiles.  Frost  and  Dzakowic  and  Hsu  assumed 
that  the  temperature  isotherm  above  the  bubble  is  distorted 
due  to  the  presence  of  the  bubble  and  proposed  Pr^  and  2.0, 
respectively,  for  the  value  of  n.  Bergles  and  Rohsenow 
suggested  n=1.0  without  considering  the  isotherm  distortion. 
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Figure  5-7  Liquid  and  vapor  equilibrium 

temperature  profiles  at  the 
onset  of  nucleate  boiling. 
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Figure  5-7  illustrates  the  above  conditions  schematically. 
Eq. (5-14)  and  Eq. (5-15)  were  solved  under  the  conditions  of 


Eg. (5-17) . 


The  predicted  wall  svmerheat 
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The  wall  superheats  for  the  present  study  calculated  by  Eq. (5- 
18)  , however,  significantly  underestimated  the  experimental 
values  (Figure  5-8)  except  those  of  Frost  and  Dzakowic. 
However,  considering  that  Pr=7.2  for  Freon  113,  it  is  unlikely 
that  the  temperature  isotherm  is  distorted  sharply  above  the 
bubble  within  a thin  thermal  layer.  According  to  Collier[58], 
the  maximum  size  of  cavities  on  the  smooth  copper  tube  is  on 
the  order  of  0.1/im.  The  critical  cavity  sizes  predicted  by 
the  above  models  are  much  larger  than  the  cavities  on  the 
heating  surface.  In  such  a case,  the  condition  given  by 
Eq. (5-17)  is  not  satisfied.  To  initiate  nucleation,  the 
liquid  and  wall  temperatures  should  increase  further  until  r„in 
in  Figure  5-7  matches  the  existing  maximum  cavity  size  r„ax,s  0*1 
the  heating  surface.  The  wall  superheat  under  this  condition 
is  given  by 


Incipient  heat  flux,  q*'inc  (W/m'^2) 
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Figure  5-8  Onset  of  nucleate  boiling  measured  in 

the  Lear jet  flight  experiment. 
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where  r„ax,s  is  calculated  from 


(5-20) 
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For  simplicity,  Hsu's  model,  in  which  the  temperature  isotherm 
without  distortion,  was  employed  in  the  following  analysis. 

The  wall  superheats  for  the  present  study  calculated  by 
Eq. (5-20)  are  tabulated  in  Table  5-3.  The  decrease  in  wall 
superheat  observed  during  the  microgravity  period  is  predicted 
reasonably  by  the  above  model.  It  is  found  that  r„ax,s  and  T„- 
Tg  are  sensitive  to  the  single-phase  heat  transfer  coefficient 
and  that  the  decrease  in  wall  superheat  is  due  to  the  increase 
in  flow  rate  during  the  microgravity  period,  not  due  to  the 
gravity  change.  In  addition,  the  nucleation  initiates  within 
the  laminar  sublayer  having  a thickness  on  the  order  of  Ijum. 
Formation  of  bubbles  in  a superheated  layer  is  a metastable 
phenomenon  and  considered  hardly  affected  by  the  change  of 
buoyancy  force. 

The  large  discrepancy  for  Case  F1015  is  caused  by  the 
flow  pattern  change  from  slug  to  annular  flow  and  the  same 
model  cannot  be  applied  to  this  case. 


Table  5-3  Critical  radii  for  nucleation  and  predicted  wall  superheats  at  the  ONE. 
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5.3.2  Pressure  Effects  on  Fluid  Properties 

The  change  of  thermal  properties  of  Freon  113  as  a 
function  of  pressure  is  illustrated  in  Appendix  C.  For 
example,  a decrease  in  the  system  pressure  for  Case  F1015 
during  the  gravity  transition  from  l.Og  to  O.Olg  was  25.7%. 
The  corresponding  changes  in  the  vapor  density  and  the  liquid 
viscosity  are  -7.1%  and  +6.4%,  respectively.  If  the  heat  flux 
is  unchanged,  the  superheat  for  nucleation  decreases  by  2.1% 
for  the  above  pressure  decrease.  This  change  is  too  small 
compared  to  the  wall  superheat  observed  in  the  experiment. 
The  same  argument  is  applicable  to  other  cases.  Therefore, 
the  change  of  properties  due  to  the  decrease  in  static 
pressure  cannot  cause  the  wall  temperature  decrease  observed 
in  the  present  study.  It  is  considered  that  the  decrease  in 
wall  superheats  near  the  ONB  was  caused  by  the  enhanced  forced 
convection  due  to  an  increase  in  mass  flux  induced  by  reduced 
system  pressure. 

The  effect  of  saturation  temperature  change  on  the 
nucleation  criterion  for  Cases  F1015,  F1020,  and  F1035  is 

illustrated  in  Figure  5-9 (a) -(c),  respectively.  The  vapor 
equilibrium  temperature  asymptotically  approaches  the 
saturation  temperature.  When  the  system  pressure  decreases, 
this  curve  shifts  downward  since  the  asymptote  Tg  is  reduced 
as  the  pressure  is  reduced.  As  the  change  of  thermal 

conductivity  was  small  and  the  heat  flux  was  constant,  the 
slope  of  the  liquid  temperature  profile  remained  unchanged. 


O.OE-KX)  1.0E-05  2.0E-05  3.0E-05  4.0E-05 

Distance  from  surface,  (m) 

(a)  Case  F1015 
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Distance  from  surface,  (m) 


(b)  Case  F1020 


0.0E*KX)  1.0E-05  2.0E-05  3.0E-05  4.0E-05 

Distance  from  surface,  (m) 


(c)  Case  F1035 

Figure  5-9  Liquid  and  vapor  equilibrium  temperatures 

near  the  heating  surface. 
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The  wall  temperatures  decreased  due  to  enhanced  bulk 
convection.  As  a result,  the  range  of  cavity  size  to  be 
activated  became  smaller.  The  change  of  r„ax  is  insignificant 
for  the  wall  superheat  because  such  a cavity  size  is  rare  on 
the  heating  surface. 

5.3.3  Superheated  Laver  Thickness  Near  the  ONB 

In  pool  boiling,  suppression  of  the  natural  convection 
induces  an  increase  in  the  thickness  of  the  thermal  layer  near 
the  heating  surface.  The  difference  in  the  thickness  of  the 
superheated  layer  in  pool  boiling  and  flow  boiling  systems 
under  microgravity  conditions  is  discussed  in  the  following 
section. 

Assuming  that  6^*  is  the  thickness  of  superheated  layer, 
one  can  obtain 

(5-22) 

<3« 

from  Eg. (5-14) . The  above  relationship  is  valid  when  a linear 
temperature  profile  is  assumed.  The  thickness  of  the 
superheated  layer  is  thus  easily  estimated  from  Eg. (5-22) , 
ranging  approximately  30  - 70/xm  for  the  present  study.  When 
the  heat  fluxes  are  constant,  6,.*  is  only  a function  of  wall 
superheat.  The  decrease  in  the  wall  superheat,  caused  by 
enhanced  bulk  convection,  reduced  the  thickness  of  the 
superheated  layer  during  the  microgravity  period. 


180 


The  above  facts  characterize  the  difference  of  heat 
transfer  at  the  ONB  between  the  pool  boiling  and  forced 
convection.  Ervin  et  al.[22]  conducted  microgravity  pool 
boiling  experiments  in  a drop  tower  facility  and  observed 
transient  nucleate  boiling  on  a transparent  gold  film.  Their 
measurement  also  included  the  wall  temperature  variation  at 
the  moment  of  boiling  inception.  They  observed  a formation  of 
large  superheated  thermal  boundary  layers  and  an  energetic 
propagation  of  bubble  growth  associated  with  the  ONB,  which 
were  unique  to  the  microgravity  conditions.  The  result  of 
wall  temperature  measurements  showed  that  during  such 
energetic  bubble  growth,  wall  superheat  was  larger  than  other 
modes  of  bubble  formation.  They  suggested  that  this  is  due  to 
larger  internal  energy  stored  inside  the  superheated  layer. 
Similarly,  Zell  et  al.[21]  reported  that  the  wall  superheat  at 
the  ONB  in  microgravity  was  higher  than  that  in  normal  gravity 
and  that  the  heat  transfer  near  the  ONB  in  microgravity  was 
deteriorated . 

Under  microgravity  conditions,  natural  convection  in  pool 
boiling  may  disappear.  This  will  cause  thickening  of 
superheated  layer  and  the  increase  in  wall  temperature  since 
only  heat  conduction  is  available  in  such  a case.  In  normal 
gravity,  the  natural  convection  assists  to  diffuse  the 
internal  energy  stored  in  the  superheated  layer  otherwise  to 
be  transferred  to  the  vapor  phase. 
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In  flow  boiling  systems  with  high  subcoolings,  forced 
convection,  instead  of  natural  convection,  carries  away  a 
portion  of  internal  energy  in  the  superheated  layer. 
Therefore,  the  energetic  bubble  formation  described  in  Ervin's 
work  do  not  occur  even  under  microgravity  conditions. 

5.4  Subcooled  Boiling  Heat  Transfer 


5.4.1  Heat  Balance 

Bowring [46]  assumed  that  the  total  heat  flux  in  the 
subcooled  region  is  described  as  a sum  of  forced  convective 
heat  flux  q"fc  and  boiling  heat  flux  q"bo 


<7tP  ■'■Qbo 

_ ^ ^ „n 

Qfc  ^ Qag  9ev 


(5-23) 


where  the  total  two-phase  heat  transfer  coefficient  is  defined 
by 


grp 


(5-24) 


The  forced  convective  heat  transfer  coefficient  is  given  by 


h 


f c 


(r„-Tr) 


(5-25) 


and  the  boiling  heat  transfer  coefficient  is  given  by 
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(5-26) 


In  Eq.  (5-23),  the  boiling  heat  flux  q''t,o  for  vertical  flow 
systems  is  further  divided  into  the  heat  flux  transferred  by 
evaporation  q”ev  and  the  heat  flux  transferred  by  bubble 
agitation  q"ag. 

In  the  existing  correlations,  the  bubble  agitation  heat 
flux  q”ag  is  implicitly  included  in  the  boiling  heat  flux. 
However,  q^g''  induced  by  bubble  motion  is  considered  to  vary 
significantly  under  variable  gravity  conditions.  Therefore, 
contributions  of  qev”  and  qag”  are  separately  discussed  in  the 
present  study.  Note  that  the  heat  transfer  coefficients  for 
two-phase  and  forced  convective  flow  are  defined  for 

(T„-Tf)  =(T^-T^)  =AT3,,+AT3,,  (5-27) 

while  that  of  boiling  contribution  is  defined  only  for  the 
superheat  (T„-Ta) . In  the  two-phase  region  where  liquid  and 
vapor  coexist,  the  forced  convective  heat  transfer  coefficient 
hfc  is  larger  than  that  in  the  single-phase  region.  Instead 
of  using  the  single-phase  Reynolds  number,  the  two-phase 
Reynolds  number  including  the  vapor  voidage  was  used 


J 

U-«i 

(5-28) 


Substituting  the  above  equation  into  the  Nusselt  number 
correlation  for  single-phase  flow,  one  obtains 
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(5-29) 


The  values  of  2.192  and  0.33,  obtained  for  the  present  study, 
were  used  as  Ci  and  C2,  respectively.  Eq. (5-29)  indicates  that 
hfc  increases  as  the  void  fraction  increases  in  the  low-quality 
region,  while  hf^  ceases  to  increase  as  the  quality  approaches 
unity.  The  physical  significance  of  the  above  equation  is 
that  the  local  heat  transfer  is  increased  by  the  reduced  flow 
area  in  the  existence  of  growing  void  fraction.  Substitution 
of  Eq. (5-24)  and  Eq. (5-29)  into  Eq. (5-23)  and  rearranging. 


Sekoguchi[33]  evaluated  the  term  (l-q”bo/<3"Tp)’^  experimentally. 
Therefore,  his  method  cannot  be  applied  instantly  to 
microgravity  heat  transfer  problems  without  considering  the 
gravity  effects  on  q”bo*  One  of  the  objectives  of  the  present 
study  is  to  investigate  the  gravity  effects  on  the  ratio 
(hTp/hfo)  • This  ratio  was  selected  as  an  index  of  comparison 
for  microgravity  boiling  heat  transfer  since  the  single-phase 
heat  transfer  coefficient  hfo  also  varies  depending  on  the  flow 
conditions.  In  the  following  section,  the  terms  [(l-x)/(l- 
a)]°'^^  and  (l~g"bo/<3"Tp) in  Eq.  (5-30)  will  be  evaluated 


gives 


(5-30) 


separately  to  examine  the  gravity  effects. 
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5.4.2  Forced  Convective  Heat  Transfer  in  the  Subcooled  Region 


The  exit  qualities  for  Cases  F1015,  F1020,  and  F1035 
under  l.Og  conditions  are  -0.4%,  -6.7%,  and  -12.1%, 
respectively.  The  negative  values  indicate  that  the  flow 
conditions  were  subcooled.  To  account  for  local  boiling  in 
the  subcooled  region,  the  actual  quality  Xg  was  introduced. 
Among  several  theories  for  Xa,  we  turn  to  Sekoguchi's  work 
since  his  analysis  is  based  on  a theoretical  model  and 
specifically  developed  for  subcooled  flow  boiling[33]. 
Assuming  the  l/7th  law  for  the  temperature  and  velocity 
profiles,  the  actual  quality  Xg  is  given  by 


x,=Bo* 


(5-31) 


where 


1 -X. 


\0.33 


{1-a) 


(5-32) 


/ 


\ 


^00 . 67  p^3/A 


2.1917 


Mf 


\0.14 


(5-33) 


Bo  = 


\ 


// 


(5-34) 


Again,  instead  of  0.8,  the  value  of  0.33  was  used  in  the 
definition  of  Bo*.  Although  the  original  expression  for  Bo* 
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uses  the  Dittus-Boelter  equation,  Eq. (5-10)  was  employed  in 
this  analysis. 

Void  fraction  and  Slip  ratio.  For  the  void  fraction, 
Sekoguchi  suggested  the  use  of  Smith's [34]  correlation 


Smith  recommended  0.4  as  the  value  for  (J).  However,  this 
factor  <j)  implicitly  includes  the  effects  of  the  slip  between 
the  liquid  and  vapor  phases  under  normal  gravity  conditions. 
Under  microgravity  conditions,  the  value  of  (J)  will  be 
different  because  the  slip  ratio  is  smaller  in  microgravity. 
Instead  of  Eq. (5-35) , Eq. (4-16)  was  employed  in  the  present 
study.  To  incorporate  the  gravity  effects  on  s,  Eq. (4-11) 
derived  by  Zuber  and  Findlay  was  used [32].  Although  the  flow 
pattern  of  interest  is  bubbly  flow,  the  above  equation  is 
widely  applicable  without  reference  to  flow  pattern. 

For  the  change  in  gravity  from  l.Og  to  2.4g  to  O.Olg,  the 
ratio  <Ugj>/<Ugj>ig  varies  from  1 to  1.257  to  .316, 
respectively.  Since  the  gravity  level  during  the  microgravity 
experiment  was  O.Olg,  the  slip  ratio  approaches  unity  but 
still  remains  larger  than  1.0,  according  to  the  drift  flux 
model . 

Eq. (5-31) , Eq. (5-32) , Eq. (5-33) , and  Eq. (5-34)  were 


Pf/Pg+<l>(V^a-l) 


(5-35) 


+ (1  -<|))  — -1 


1/2 


Pf 


1+<|)(1/X3-1) 


solved  simultaneously  to  obtain  Xg  and  a for  variable 
gravities.  Since  the  equations  were  recursive,  17  to  18 
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Quality,  Xth 


(a)  Case  F1015 


(b)  Case  F1020 


(c)  Case  F1035 


Figure  5-10  Thermal  equilibrium  qualities  for  the 

Lear jet  flight  experiments. 
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iterations  were  needed  before  convergence  using  the  bi-section 
method.  The  thermal  equilibrium  qualities,  void  fractions, 
and  actual  qualities  evaluated  at  the  exit  of  the  test  section 
for  each  case  were  calculated  simultaneously.  The  actual 
qualities  ranged  from  0.0068% (Case  F1035)  to  0.970%  (Case 
F1015) . 

Figure  5-10 (a)  through  (c)  indicate  the  evolution  of 
thermal  equilibrium  qualities  along  the  test  section  for  Cases 
F1015,  F1020,  and  F1035,  respectively.  The  qualities 
increased  linearly  since  the  heat  input  through  the  wall  was 
uniform.  The  inlet  qualities  under  microgravity  conditions 
decreased  since  the  flow  rates  increased  due  to  the  decrease 
in  static  pressure.  The  change  in  gradient  of  each  curve  is 
caused  by  the  liquid  enthalpy  change  due  to  varied  pressure 
and  liquid  temperature. 

The  profiles  for  the  actual  qualities  are  different  from 
thermal  equilibrium  qualities (Figure  5-11 (a)  through  (c) ) . 
The  actual  quality,  defined  as  the  ratio  of  the  enthalpy 
stored  in  the  superheated  layer  to  the  total  enthalpy  carried 
by  flow,  is  given  by 


yi 

CpJ  (r-TJ  Ul  (l-y*)dy* 

0 

I 

h,Ju^(l-y*)dy* 

0 


(5-36) 


where  ys*=(ys/R)  is  the  dimensionless  superheated  layer 
thickness.  Eq. (5-36)  implicitly  includes  the  effects  of 
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Actual  quality,  Xa 


Figure  5-11 


Actual  quality  profiles  calculated  by 
the  modified  Sekoguchi's  method. 
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Void  fraction 


(a)  Case  F1015 


(b)  Case  F1020 


(c)  Case  F1035 

Void  fraction  profiles  predicted  by 
Sekoguchi ' s and  Ahmad's  methods. 


Figure  5-12 
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gravity  through  the  void  fraction  and  liquid  velocity.  Under 
microgravity  conditions,  the  actual  qualities  for  each  case 
were  lower  than  those  of  normal  gravity  since  the  increased 
bulk  convection  removes  some  part  of  the  enthalpy  otherwise 
built  up  in  the  superheated  layer. 

The  void  fraction  profiles  along  the  test  section  are 
shown  in  Figure  5-12 (a)  through  (c)  . The  void  fraction 
becomes  largest  when  the  slip  ratio  is  1.0.  However,  the 
prediction  by  Eg. (4-16)  indicates  that  the  above  assumption 
overestimates  the  void  fraction.  This  results  suggests  that 
there  still  remains  a relative  velocity  between  the  phases 
even  in  O.Olg.  The  void  fraction  calculated  from  Ahmad's 
correlation,  which  is  only  valid  under  normal  gravity 
conditions,  is  also  shown  for  reference. 

Effect  of  Superheated  Laver  Thickness.  It  is 
particularly  interesting  to  examine  gravity  effects  on  the 
thickness  of  the  superheated  layer,  in  the  subcooled 
region.  Assuming  the  l/7th  law  for  the  liquid  temperature 
profile,  y^*  is  given  as 


where  is  a constant (=.83 3)  and  Bo*  is  defined  in  Eq.  (5-33)  . 
Bo*  includes  the  void  fraction  inside  the  term  [(1-Xa)/(1- 
a)]°'^^  and  is  affected  by  gravity  through  the  change  in  slip 
ratio.  The  calculated  superheated  layer  thicknesses  are  shown 
in  Figure  5-13 (a)  through  (c)  for  Cases  F1015,  F1020,  and 


(5-37) 


ys,  Rmax  (m) 


(a)  Case  F1015 


(b)  Case  F1020 
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Figure  5-13  Growth  of  superheated  layer  and  bubble 

diameter  in  the  axial  direction. 
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F1035,  respectively.  Although  y^*  grows  larger  along  the  test 
section,  y^*  decreased  only  slightly  in  microgravity . This 
indicates  that  the  build-up  of  enthalpy  in  the  superheated 
layer,  which  was  observed  in  pool  boiling,  would  not  occur  in 
flow  boiling  since  the  enthalpy  near  the  wall  is  carried  away 
by  forced  convection.  This  is  also  supported  by  the  following 
facts: 

a)  The  contribution  of  natural  convection  to  the  total  heat 
flux  was  less  than  1%  for  the  present  study. 

b)  The  radial  velocity  profiles  calculated  in  Chapter  4 show 
that  the  change  of  the  liquid  temperature  profile  near  the 
wall  in  microgravity  was  small  for  turbulent  flow. 

Using  the  actual  qualities  and  void  fractions  obtained 
through  the  above  method,  the  factor  [ (l-x^)/ (1-a)  was 
evaluated  along  the  test  section.  Figure  5-14  (a) -(c)  shows 
axial  evolutions  of  this  factor  for  each  case.  The  value  of 
[ (I-X3)/ (1-a)  increases  rapidly  near  the  exit  because  the 
gradient  of  quality  dx^/dz  is  nearly  constant  throughout  the 
test  section,  while  da/dz  increases  rapidly  in  the  latter  half 
of  the  test  section.  Due  to  significant  increase  in  void 
fraction,  the  values  of  [ (l-x^)/ (1-a)  in  microgravity 
exceeded  those  in  normal  gravity  for  Cases  F1015  and  F1020. 
However,  for  Case  F1035  in  which  the  void  fraction  as  well  as 
the  true  quality  was  smaller  than  other  cases,  [(1-Xa)/(1- 
(X)]°-33  close  to  1.0.  This  implies  that  single-phase  heat 
transfer  was  still  dominant  for  Case  F1035. 


[(1-Xa)/(1-a)  ]^.33 


(b)  Case  F1020 


Variation  of  forced  convective  factor 
under  microgravity  conditions. 


193 


Figure  5-14 
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5«4.3  Boiling  Heat  Flux 

Estimation  of  the  ratio  of  boiling  heat  flux  to  the  total 
heat  flux  q"bo/q''w  from  the  present  data  sets  is  difficult 
since  the  boiling  conditions  such  as  bubble  diameters  or  the 
number  of  nucleation  sites  are  unknown.  However,  the  flow 
conditions  such  as  flow  rates,  system  pressures,  and  qualities 
were  measurable  with  certain  accuracy.  Therefore,  q”bo  can 
be  obtained  more  accurately  by  subtracting  q"fc  from  the  total 
heat  flux,  rather  than  calculating  q''bo  independently.  Figure 
5-15  (a)  through  (c)  show  the  distributions  of  q"bo/q"w  along 
the  flow  axis  for  Cases  F1015,  F1020,  and  F1035,  respectively. 
For  all  cases,  the  forced  convection  ratio  qfc'Vqw"  increased 
in  the  subcooled  region  near  the  outlet  of  the  test  section 
under  microgravity  conditions.  Conversely,  the  boiling  heat 
flux  ratio  qbo'Vqw"  decreased  because  the  bubble  agitation 
component  q^g”  became  smaller.  The  evaporation  heat  flux  qg^" 
calculated  by  the  mass  transport  model  was  found  to  be  much 
smaller  than  any  other  component  in  the  present  study, 
although  qgv”  was  slightly  increased  in  microgravity. 

Figures  5-16 (a)  through  (c)  show  the  boiling  curves  for 
the  subcooled  boiling  under  microgravity  conditions  using 
Rohsenow's  correlation.  The  value  of  g was  set  to  0.098m/s^ 
for  microgravity  conditions.  Rohsenow's  correlation  is  given 
by 


(q’*bo/q'*w)  and  (q**fc/q**w) 
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(a)  Case  F1015 


Figure  5-15 


Development  of  boiling  heat  flux  in  the 
axial  direction. 
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(b)  Case  F1020 


Figure  5-15 


continued . 
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(c)  Case  F1035 


Figure  5-15 


continued . 


(Zviu/AA) 
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(a)  Case  F1015 

Figure  5-16  Boiling  curves  under  variable  gravity  conditions 
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(b)  Case  F1020 


Figure  5-16 


continued 
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(c)  Case  F1035 


Figure  5-16 


continued 


201 


1 

gbo  b 

^fq 

■q 

ll/3 


Pr 


1.7 


(5-38) 


where  1/Cf  is  a constant  depending  on  the  surface  condition 
and  b is  the  characteristic  diameter  of  bubbles  given  by 


b=  2 (5-39) 

N g(Pf-pg) 

The  above  equation  is  derived  as  a result  of  the  force  balance 
between  the  surface  tension  and  the  buoyancy  force,  a 
classical  model  of  bubble  detachment  in  pool  boiling.  The  use 
of  this  model  in  subcooled  boiling  heat  transfer  is  not 
physically  sound  and  the  predictions  grossly  underestimate  the 
wall  superheat  at  the  ONB  under  microgravity  conditions.  The 
force  balance  providing  the  detachment  criterion  should  be 
revised  for  microgravity  applications. 

Evaporation  heat  flux  a”„„.  In  order  to  investigate  the 
influence  of  gravity  on  the  governing  parameters  of  bubble 
growth,  it  is  necessary  to  develop  a model  to  describe  the 
heat  transfer  mechanism  for  a single  bubble  in  subcooled 
boiling.  The  mass  transport  model,  first  proposed  by 
Snyder[45],  was  adopted  in  the  present  analysis.  The 
importance  of  mass  transfer  inside  the  bubble  under 
microgravity  conditions  was  suggested  by  Zell  et  al.[21].  It 
is  reported  that  subcooled  pool  boiling  was  effectively 
continued  in  a reduced  gravity  environment  where  the  buoyancy 
force  was  absent.  It  is  interesting  to  note  that  heat  was 
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continuously  transferred  from  the  wall  even  if  the  bubble 
stays  on  the  surface.  If  a single  bubble  carries  the  enthalpy 
equivalent  to  one  bubble  volume,  heat  transfer  in  pool  boiling 
would  become  immediately  inefficient  under  microgravity 
conditions  where  bubble  detachment  no  longer  takes  place. 
This  is  a strong  evidence  to  support  the  mass  transport 
mechanism  for  their  subcooled  pool  boiling  experiment. 

The  mass  transport  model  was  then  incorporated  in  the 
bubble  growth  model  proposed  by  Unal[50].  Dynamics  of  bubble 
growth  in  subcooled  flow  is  determined  by  the  heat  balance. 
Upon  the  initiation  of  a bubble  at  the  activated  nucleation 
site,  heat  is  mainly  supplied  to  the  bottom  part  of  the  bubble 
by  both  heat  conduction  through  the  interface  and  evaporation 
from  the  microlayer  beneath  the  bubble.  Simultaneously,  heat 
is  removed  at  the  top  of  the  bubble  that  is  exposed  to  the 
low-temperature  liquid.  This  condensation  mechanism  is 
observed  only  for  subcooled  boiling  and  allows  an  extremely 
high  heat  flux  in  the  subcooled  boiling  region. 

Unal[48]  developed  a bubble  growth  model  based  on 
microlayer  evaporation.  The  heat  balance  between  evaporation 
from  the  microlayer  and  condensation  at  the  upper  part  of  the 
bubble  was  considered  and  the  conduction  through  the  bubble 
interface  was  neglected  since  the  conduction  heat  transfer  is 
small  compared  to  evaporation.  In  order  to  incorporate  the 
mass  transport  mechanism,  Unal's  model  was  modified  with  the 
following  assumptions: 
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a)  The  liquid  flow  was  fully  developed  and  the  velocity  and 
temperature  profiles  were  similar  to  those  of  single- 
phase flow.  This  assumption  was  justified  since  the 
number  of  nucleation  sites  was  relatively  small  for  the 
flow  boiling  system  with  large  subcooling  and  low  heat 
fluxes. 

b)  Heat  flux  was  constant  and  wall  temperatures  were  given 
as  input  data  to  compute  the  maximum  bubble  diameter  and 
growth  rate. 

c)  For  simplicity,  a constant  contact  angle  was  used. 
Deformation  of  the  bubble  during  the  growth  period  was 
neglected. 

d)  Inertia-controlled  bubble  growth  to  be  observed  at  the 
beginning  of  bubble  growth  was  neglected  since  this 
period  was  extremely  short  compared  to  the  period  of 
asymptotic  growth  mode  that  followed. 

e)  The  dry  area  in  the  microlayer  was  neglected.  As  van 
Stralen[49]  pointed  out,  the  dry  patch  in  the  microlayer 
was  hardly  observed  at  atmospheric  pressure. 

f)  The  interaction  between  bubbles  was  neglected. 

g)  The  thickness  of  the  superheated  layer  was  calculated  by 
the  temperature  profile  obtained  by  the  l/7th  law. 
Formulation  of  the  model.  A schematic  diagram  of  the 

proposed  model  is  shown  in  Figure  5-17.  In  this  figure,  b^ 

indicates  the  relative  height  of  the  superheated  layer 

thickness  to  the  bubble  height.  The  lower  part  of  the  bubble 
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Figure  5-17  Schematic  diagram  of  the  heat  balance 

for  a single  bubble  in  subcooled  flow. 
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up  to  y=2R*bh  is  immersed  in  the  superheated  layer,  while  the 
upper  part  from  y=2R*bh  is  subjected  to  the  subcooled  liquid. 
The  energy  balance  for  a single  bubble  on  the  heating  surface 
in  the  subcooled  flow  is  described  as 


g."  • A,  -gio.  • A A,  (5->o) 

where  q"bo  is  the  boiling  heat  flux  provided  from  the  heating 
surface,  q”con  the  condensation  heat  flux  removed  from  the 
bubble,  Ag  the  base  area  of  the  bubble  available  for 
evaporation,  Ac  the  area  exposed  to  subcooled  liquid,  V the 
vapor  volume,  Pg  the  vapor  density,  and  hfg  the  latent  heat  of 
vaporization.  The  transient  conduction  heat  flux  q”bo  is 
obtained  from  the  approximate  problem  in  which  a semi-infinite 
liquid  column  of  temperature  T^  makes  contact  on  the  heating 
surface  of  temperature  T„ 


-I 

<?ev  = - 


K6, 


f ''0 

2/  (TiaTp^ 


(5-41) 


where  kf  the  thermal  conductivity  of  the  liquid,  a the  thermal 
diffusivity  of  liquid,  and  0q  is  the  superheat  (T„-Ti)  in  the 
microlayer.  The  condensation  heat  flux  q"con  derived  by 
Unal[50]  for  subcooled  flow  at  O.lMPa  was  used  in  this 
analysis 


//  . . Chf(D/2) 

r —h  At  — fg  ' / / 

(con  ^^con“-^3ub  

'^fg 


(5-42) 


where  C is  a constant  obtained  experimentally  equal  to  63.4 
for  the  present  study,  and  ^Tsub  the  subcooling  at  the  top  of 
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the  bubble.  The  value  of  ATgui,  was  considered  zero  in  the 
superheated  layer,  while  outside  the  superheated  layer,  it  was 
calculated  from  the  previously  assumed  temperature  profile  of 
the  liquid.  Considering  the  geometric  factor  and  rearranging, 
Eq. (5-40)  is  expressed  as 


dR 

~Hi 

where  f ^ , f 2 , f 3 
contact  angle  (jjg 


CAT 


sub 


■'•2 


(l"Pg/Pf) 
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a Ja 


fA 


V ^ y 


\ 


t-1/2 


(5-43) 


are  geometric  factors  and  functions  of  the 


^ sin2<l)o 

(5-44) 

^ 2 

-f2  = 1 cos<j)o  - 2 jb^ 

(5-45) 

= 1 + 1 cos<j)o  - 1 COS^<J)o 

(5-46) 

When  ATgub  is  zero,  Eq.  (5-43)  reduces  to  a diffusion  bubble 
growth  in  saturated  boiling  with  ^Q=7i/2 


where  Ja  is  the  Jacob  number  defined  as 


(5-47) 


7^  _ ( ^ ^s)  ^-plPf 

■ 


(5-48) 


When  subcooling  is  not  negligible,  Eq. (5-43)  becomes  a first- 
order  non-linear  differential  equation  since  AT^ub  is  a 
function  of  the  bubble  radius  R.  Eq. (5-43)  was  solved 
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numerically  by  using  the  4th  order  Runge-Kutta  method  under 
the  following  initial  condition 

R=0  at  t=0  (5-49) 

The  initial  radius  of  the  bubble  is  on  the  same  order  of  the 
cavity  size  ( «1  /xm)  and  was  neglected  compared  to  that  of  the 
bubble  under  growth. 

To  calculate  the  superheat  of  the  evaporation  microlayer 
00,  the  surface  temperature  T^  needs  to  be  specified.  The 
initial  thickness  of  the  evaporation  microlayer [49]  do  is 
described  as 


(5-50) 


where  hf<-  denotes  the  heat  transfer  coefficient  in  the 
subcooled  region.  <p  denotes  the  ratio  of  wall  superheat  to 
the  total  temperature  difference  given  by 
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00+00 


(5-51) 


where  0o=(T„-Tsat)  and  0o*=(Tsat”Tf)  . It  is  assumed  that  this 
thickness  remained  almost  the  same  throughout  the  bubble  life 
since  the  evaporation  rate  was  small[49]. 

To  account  for  mass  transport  inside  the  bubble,  the  mass 
flux  due  to  evaporation  needs  to  be  determined.  The 
evaporation  rate  between  two  liquid  surfaces  at  different 
temperatures [45]  is  given  by 
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(5-52) 


where  m is  the  mass  transfer  rate  (kg/sec) , Ag  the  area  of 
evaporation  (m^) , the  surface  temperature  of  the  superheated 
microlayer  (K) , Tg  the  temperature  of  vapor  inside  the  bubble 
(K)  , R the  universal  gas  constant  (8 . 3143J/mol/K) , the 
molecular  weight  of  Freon  113  (187.376kg/mol) , and  a the 
accommodation  coefficient  which  accounts  for  freshness  of  the 
surface  and  considered  unity  for  such  a case  as  nucleate 
boiling.  The  vapor  equilibrium  temperature  Tg  was  determined 
by  Eg.  (5-19)  . The  heat  flux  due  to  evaporation  is,  therefore, 
given  as 


*fg^=giv  (5-53) 

The  above  amount  also  must  be  equal  to  the  heat  conduction 
through  the  microlayer  to  maintain  the  constant  heat  flux 


h 


fg 


Th 


(5-54) 


Here,  the  heat  capacity  of  the  evaporation  microlayer  is 
neglected  since  its  thickness  is  on  the  order  of  few  microns. 
The  only  unknown  parameter  in  the  above  equation  is  the 
interface  temperature  T^.  Eq. (5-52)  and  Eq. (5-54)  was  solved 
numerically  by  the  bi-section  method  to  obtain  the  interface 
temperature  Ti.  It  was  found  that  Ti  was  higher  than  the  vapor 
temperature  Tg  by  a fraction  of  degree  for  the  present  study. 


Rb  (m) 
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(a)  g/g0=1.0 


(b)  g/g0=0.01 


Typical  bubble  growth  predicted  by  the 
mass  transport  model  for  Case  F1035. 


Figure  5-18 
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The  measured  wall  temperatures  were  used  to  give  the 
superheat  0q  in  the  above  calculation.  Typical  bubble  growth 
for  Case  F1035  is  shown  in  Figure  5-18.  Within  the 
superheated  layer,  the  bubble  grows  proportional  to  Its 
growth  rate  decreases  after  the  top  portion  of  the  bubble  was 
exposed  to  the  subcooled  liquid.  Unlike  saturated  boiling, 
the  bubble  size  is  limited  by  the  thickness  of  the  superheated 
layer.  In  addition  to  the  wall  superheat,  the  growth  rate  is 
determined  by  the  growth  factor  2\/a*Ja.  The  variation  of  this 
factor  as  a function  of  pressure  is  shown  in  Appendix  C.  For 
example,  this  factor  increases  by  approximately  19.5%  for  the 
pressure  change  experienced  in  Case  F1020.  This  indicates 
that  a bubble  grows  faster  and  larger  within  the  same  period 
under  reduced  pressure  conditions,  enhancing  the  evaporation 
heat  transfer. 

In  the  mass  transport  model,  evaporation  and  condensation 
processes  continue  until  the  bubble  is  detached  from  the  wall 
either  by  the  drag  force  of  the  liquid  flow  or  by  the  buoyancy 
force.  At  the  moment  of  bubble  detachment,  the  growth  rate 
approaches  zero.  Therefore,  in  the  present  analysis,  it  is 
assumed  that  the  bubbles  leave  the  heating  surface  when  the 
maximum  diameter  is  reached. 

The  local  evaporation  heat  flux  q”ev  calculated  from 
Eg. (5-53)  was  extremely  high.  For  example,  the  calculated 
local  evaporation  heat  flux  at  the  inlet  and  the  exit  are 
86.8kW/m^  and  27.2kW/m^,  respectively  for  Case  F1015.  This  is 
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because  the  heat  resistance  becomes  extremely  low  in  the 
evaporation  microlayer.  The  boiling  heat  flux  in  Eq.  (5-38)  is 
the  averaged  heat  flux  over  the  entire  heat  transfer  area. 
The  relationship  between  the  local  evaporation  heat  flux  and 
the  average  evaporation  heat  flux  is  given  by 

giv,  avg  = T|  giv,  local  ( 5-55  ) 

where  i)  is  the  area  ratio  of  evaporation  area  to  the  total 
heat  transfer  area,  i]  is  defined  as 

= (5-56) 

/ij. 

where  is  the  base  radius  of  a single  bubble,  N the  total 
number  of  nucleation  sites  on  the  surface,  and  is  the  total 
heat  transfer  area.  Since  is  very  small  compared  to  1.0  in 
subcooled  boiling,  q”ev,  avg  becomes  a fraction  of  the  total 
heat  flux.  In  the  following  text,  q”e^  will  indicate  the 
average  value  unless  otherwise  specified. 

The  evaporation  process  is  sensitive  to  pressure  change 
but  hardly  affected  by  the  gravity  change  in  the  subcooled 
flow  boiling.  The  initial  thickness  of  the  evaporation 
microlayer  is  slightly  increased  when  bulk  convection  is 
deteriorated  in  microgravity,  but  the  influence  on  the  overall 
heat  transfer  is  insignificant.  The  growth  rate  is  indirectly 
increased  when  the  wall  superheat  is  elevated  as  a result  of 
deteriorated  bulk  convection  in  microgravity.  However,  under 
the  conditions  where  bubble  detachment  is  hydrodynamically 
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controlled,  the  effect  of  gravity  change  on  heat  transfer  will 
be  significant. 

Heat  flux  due  to  bubble  agitation,  For  a vertical 
subcooled  flow  boiling  system  under  normal  gravity  conditions, 
the  distributed  bubbles  rising  faster  than  the  surrounding 
liquid  phase  introduce  turbulent  motion  in  addition  to  eddies 
inherent  to  turbulent  flow.  It  is  considered  that  the  above 
turbulent  effect  due  to  bubble  motion  becomes  smaller  under 
microgravity  conditions  since  the  slip  velocity  of  the  vapor 
phase  approaches  zero.  However,  coalescence  of  the  detached 
bubbles  is  promoted  in  the  absence  of  relative  velocity  and 
resulting  bubble  size  grows  larger  than  that  of  normal 
gravity. 

As  is  discussed  in  Chapter  4 , the  turbulence  caused  by 
bubble  agitation  is  proportional  to  the  bubble  diameter  dg  and 
the  slip  velocity  Ug.  Under  extremely  small  gravity 
conditions,  Ug  approaches  zero  and  there  is  no  turbulence 
caused  by  bubble  agitation.  In  intermediate  low-gravity 
conditions,  however,  the  influence  of  dgXUg  on  eddy  generation 
is  unknown.  Besides,  the  decrease  in  static  pressure  changes 
the  shear  stress  distribution  in  the  liquid  flow.  These 
differences  are  especially  pronounced  in  vertical  flow  systems 
under  reduced  gravity  conditions  and  alter  the  velocity  and 
temperature  profiles. 

The  relative  magnitude  of  q"ey  to  q”3g  is  hard  to  define 
with  the  limited  data  available  from  the  present  experiments. 
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Bowring  proposed  an  empirical  relationships  for  the  ratio 
(q''ag/^3"ev)  analyzing  experimental  void  fraction  data  for 
water [46] 


Qag 

~Zi~i 

Qev 


= 1+3 


(5-57) 


The  above  equation  is  valid  for  l-9.5bar  (14.7-140psia) . 
Substituting  the  thermal  properties  of  Freon  113  at  the 
saturation  temperature,  the  ratio  becomes  4.1  for  Case  F1035. 
This  indicates  that  the  heat  flux  due  to  bubble  agitation  q”ag 
is  about  4.1  times  larger  than  the  evaporation  heat  flux  q”ev* 
However,  the  ratio  € obtained  from  the  mass  transport  model 
was  much  larger  than  4.1  as  shown  in  Figure  5-15.  This  is 
because  the  heat  flux  used  in  the  present  study  was  smaller 
than  that  of  Bowring's  experiment.  His  model  predicts  that 
q"ag  prevails  in  the  majority  of  the  subcooled  region  under 
normal  and  microgravity  conditions. 

The  empirical  factor  3.2  in  Eg. (5-57) , an  effective 
temperature  rise  of  the  water  pulled  in  and  pushed  out  by  the 
bubble,  is  considered  a function  of  bubble  size.  The  ratio  € 
will  be  larger  in  microgravity  because  the  bubble  agitation 
heat  flux  becomes  smaller.  Another  experiment  to  investigate 
the  influence  of  gravity  exclusively  on  q"ag  is  proposed  and 
described  in  Chapter  6 . 

g"„„  predicted  by  Larsen-Tona * s correlation.  A semi- 
analytic  model  to  account  for  the  evaporation  heat  flux  in 
subcooled  flow  boiling  was  proposed  by  Larsen  and  Tong[36]. 
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This  model  is  based  on  the  formation  and  growth  of  the  bubble 
boundary  layer  considering  a heat  balance  between  the  heat 
removed  by  the  liquid  core  and  the  external  heat  input.  The 
heat  flux  removed  by  liquid  flow  at  the  edge  of  the  bubble 
boundary  layer  q**2,  is  given  by 


16.5  l0.8\  (r2/-R)0c  ^^-0.3 

8 \ Pr  jTniT-r^/R) 


(5-58) 


where  r2  is  the  distance  from  the  center  line  to  the  edge  of 
the  bubble  boundary  layer  and  is  the  average  void  fraction 
in  the  bubble  boundary  layer.  Larsen  and  Tong  assumed  that  r2 
is  equal  to  the  laminar  sublayer  thickness  of  fully  developed 
turbulent  liquid  flow.  In  the  present  analysis,  r2  is 
computed  using  the  superheated  layer  thickness,  while  Ob  is 
obtained  from 


a^=a/(l- (r,/R)^)  (5-59) 

q”2  represents  the  single-phase  convection  in  a flow  area 
gradually  reduced  by  the  growing  bubble  boundary  layer.  Note 
that  this  theory  is  developed  for  the  horizontal  flow  and  the 
bubble  agitation  mechanism  normally  observed  in  a vertical 
flow  system  is  not  included.  Thus,  the  evaporation  heat  flux 
is  expressed  as  the  difference  between  the  total  heat  flux  and 
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(a)  Case  F1015 


(b)  Case  F1020 


(c)  Case  F1035 

Figure  5-19  Contribution  of  evaporation  heat  flux 

predicted  by  Larsen-Tong'  correlation. 
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qi^=qi'-{r,/R)q^'  (5-60) 

Eq. (5-58)  was  applied  to  Cases  F1015  through  F1035  and  the 
results  are  shown  in  Figure  5-19.  It  also  should  be  noted 
that  Eq.  (5-58)  does  not  include  gravity  effects  explicitly. The 
calculation  results  for  Cases  F1015,  F1020,  and  F1035  are, 

therefore,  regarded  as  those  for  hypothetical  experiments 
performed  under  normal  gravity  using  the  same  flow  conditions. 
The  predicted  evaporation  rate  was  enhanced  mainly  because  of 
reduced  pressure,  not  reduced  gravity. 

This  correlation  tends  to  overestimate  the  values  for 
evaporation  heat  flux  since  the  empirical  constants  are  chosen 
for  fully  developed  subcooled  flow  boiling  with  large  heat 
fluxes. 

predicted  bv  Chen's  correlation.  The  well-known 
Chen's  correlation  for  saturated  boiling  heat  transfer  was 
employed  to  estimate  the  subcooled  boiling  heat  transfer  for 
the  present  experiment.  Setting  the  suppression  factor  S 
equal  to  1.0,  the  nucleate  boiling  heat  transfer  coefficient 
hbo  for  subcooled  boiling  was  calculated  using  the  measured 
wall  temperature  data,  h^o  is  expressed  as 


=0.00122 


.0.79  „0. 45  _0.49 
■f  Cpi  Pf 


_0  5 . . 0 . 2 9 


0.24  ^0.24 
Ifg  Pg 


i0.24  A ^0.75 
sat  ^-^sat 


(5-61) 


where  AT^a^  denotes  the  wall  superheat,  while  AP^at  is  the 
saturation  pressure  difference  evaluated  at  wall  temperature 
and  at  the  bulk  temperature.  The  results  are  shown  in  Figure 


(a)  Case  F1015 


(b)  Case  F1020 


(c)  Case  F1035 

Figure  5-20  Contribution  of  boiling  heat  flux 

predicted  by  Chen's  correlation. 
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5-20.  It  should  be  noted  that  the  boiling  heat  flux  h^o 
consists  of  evaporation  and  bubble  agitation.  The  predicted 
values  tend  to  overestimate  the  heat  transfer  rate  because 
this  correlation  was  developed  for  saturated  flow  boiling  with 
high  heat  fluxes  where  the  evaporation  heat  flux  is  dominant 
over  the  bubble  agitation  heat  flux.  These  calculation 
results  are  used  to  estimate  the  pressure  effects  on  hto 
without  considering  the  influence  of  gravity  on  boiling. 

The  predicted  values  of  q''t,o  in  the  outlet  region  were  compared 
with  the  measurements  and  tabulated  in  Table  5-4.  Since  the 
Chen's  correlation  includes  no  effect  of  gravity,  these 
calculation  results  are  interpreted  as  the  data  that  would  be 
obtained  in  hypothetical  experiments  performed  under  normal 
gravity  conditions  using  the  same  flow  conditions.  Chen's 
correlation  predicts  an  increase  in  the  boiling  heat  flux 
under  microgravity  conditions,  which  is  opposite  to  the 
measurements.  This  is  because  qg^"  is  dominant  in  Chen's 
experiments,  while  q^g"  is  dominant  in  the  present  experiments. 
The  increase  in  q^v"  under  microgravity  conditions  as  predicted 
by  Larsen-Tong ' s model  agreed  with  the  present  measurements. 

5.5  Gravity  Effects  on  Boiling  Mechanisms 

The  total  heat  flux  consists  of  contributions  of 
evaporation,  bubble  agitation,  and  forced  convection  of  the 
liquid  phase.  For  flow  boiling  systems  provided  with  a 
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constant  heat  flux,  the  following  heat  balance  must  be 
satisfied 


q!,'  = const.  =qi'  +g/g  +g/c 
If  the  ratio  f is  defined  as  follows, 

f = 

( gig  + qlc ) 

Eq. (5-63)  is  rewritten  as 


1 

e 


+ 


\ 

qtc 


1 1 


j 


(5-62) 


(5-63) 


(5-64) 


where  e is  the  ratio  of  the  agitation  heat  flux  to  the 
evaporation  heat  flux  as  defined  in  Eq. (5-57) . The  value  of 
f near  the  test  section  exit  was  calculated  for  Case  F1035  for 
various  values  of  Ug  in  Section  4.3.  Under  the  conditions  of 
Case  F1035,  the  value  of  € was  equal  to  4.1  in  normal  gravity. 
For  simplicity,  the  same  value  for  e was  assumed  in  the 
microgravity  calculations.  Eq. (5-64)  yields  the  relative 
magnitudes  of  qev"/  qag"»  aud  qf^”  when  appropriate  values  are 
substituted  for  f and  € . 

The  resulting  fractions  of  each  contribution  in  the 
slightly  subcooled  region  for  Case  F1035  are  tabulated  in 
Table  5-5.  Note  that  the  method  described  here  is  applied 
only  to  bubbly  flow  in  the  subcooled  boiling  at  the  exit  of 
the  test  section.  The  selection  of  Case  F1035  was  made  based 
on  the  least  effect  of  pressure  change.  Under  normal  gravity 
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conditions,  qj^''  and  q^g"  account  for  54%  and  36%  of  the  total 
heat  flux,  respectively.  The  contribution  of  the  evaporation 
heat  flux  was  9%.  Under  microgravity  conditions,  the  above 
fractions  change  significantly,  mainly  due  to  the  increased 
void  fraction  and  smaller  slip  velocity.  In  Table  5-5,  three 
different  cases  were  calculated  by  changing  the  value  of  Ug. 
The  fraction  of  q^g''  increases  as  the  magnitude  of  the  eddy 
generated  by  bubbles  increases.  The  case  for  UB=0.015m/s 
appears  to  most  closely  simulate  the  actual  case  since  this 
value  is  close  to  the  bubble  rising  velocity  computed  by  the 
following  relationship  proposed  by  Peebles  and  Garber [51] 


= 0. 33  g 


0.76 


\0.52 

Pf 


iMf 


1.28 


/ 


for  2^i?e,^^4 . 02 


(5-65) 


where  the  bubble  Reynolds  number  Re,,  and  Y are  defined  by 
2Rt,UbPf/Mf  and  gMfV(Pf®^)»  respectively.  The  fractions  of  qf<.” 
and  qbo”  in  normal  and  microgravity  are  shown  Figure  5-21  and 
compared  with  the  measurements.  The  boiling  heat  flux  q”bo  in 
the  present  experiment  was  obtained  using  Sekoguchi's  method 
and  given  by 




Q 

1 

f 1 - xJ 

qj' 

j 

i"-aj 

( ^ - ) ONB 

(5-66) 


where  (T„-Tf)exit  and  (T„-Tf)oNB  were  the  temperature  differences 
measured  at  z=0.468m  from  the  test  section  inlet  and  at  the 
ONB,  respectively.  The  fraction  of  qf^”  in  microgravity 


(q'fc/q“w)  and  (q‘1>o/q’'w)  (q'Yc/q'Vr)  and  (q“bo/q‘‘w) 
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(a)  g/g0=1.0 


(b)  g/g0=0.01 


Figure  5-21 


Ratios  of  heat  flux  components  for  Case 
F1035  predicted  by  the  present  analysis 


224 


increased  slightly  compared  to  that  of  normal  gravity. 
However,  the  fraction  of  qag"  decreased  significantly  due  to 
loss  of  the  buoyancy  force.  The  forced  convection  heat  flux 
qfc"  accounted  for  70.7%  of  q„”,  while  the  bubble  agitation 
heat  flux  qag”  accounted  for  23.6%.  The  heat  flux  due  to 
evaporation  qg^”  was  5.7%. 

Based  on  the  fractions  computed  above,  the  ratio  of  the 
boiling  heat  flux  (qbo'V<3w”)  was  calculated  for  Case  F1035  to 
obtain  the  two-phase  heat  transfer  ratio  (hTp/hfo)  • For  the 
subcooled  flow  boiling  with  small  void  fractions,  the 
influence  of  gravity  on  the  forced  convection  term  [(1-Xa)/(1- 
jj)]0.33  negligible.  As  a result,  the  term  (l-qbo”/gw") 

became  dominant  in  determination  of  (hTp/hfo)  under  the 
conditions  of  Case  F1035.  This  result  is  plotted  in  Figure  5- 
22  and  indicates  that  the  predictions  qualitatively  agree  with 
the  measurements  under  the  conditions  of  F1035. 

The  two-phase  flow  heat  transfer  has  been  treated  as  a 
function  of  the  Lockhart-Martinelli  parameter  Xti  and  the 
Boiling  number  defined  by  Bo=q„"/ (G- hfg)  . The  ratios  hTp/hp, 
obtained  under  normal  and  microgravity  conditions  are  also 
plotted  in  Figure  5-22  as  a function  of  l/Xtp.  The  same 
experiments  were  repeated  12  times  on  the  ground  under  the 
conditions  similar  to  those  of  the  flight  experiment.  The 
atmospheric  pressure  in  the  laboratory  was  latm,  while  that  in 
the  Lear jet  cabin  during  the  experiment  was  0.85atm.  Although 
the  number  of  data  was  small,  the  ratios  hTp/hp,  obtained  in 
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Figure  5-22  Two-phase  flow  heat  transfer  rate  measured 

under  normal (l.Og)  and  microgravity (0 . Olg) 
conditions. 
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microgravity  tend  to  become  smaller  than  those  of  normal 
gravity.  A correlation  derived  by  Pujol  and  Stenning  for 
Freon  113  in  a vertical  up-flow  system  is  shown  as  well [52]. 
Their  correlation  is  given  by 


= 0.9 

BoxlO^  +4.45 

1 1 

0.37 

(5-67) 


Prediction  made  by  Sekoguchi's  method  without  modification  is 
also  plotted  in  Figure  2-22.  Both  correlations  overestimated 
the  ratio  hTp/hi,  mainly  because  their  correlations  were  derived 
for  higher  flow  rates  and  higher  heat  fluxes. 

The  boiling  heat  flux  qbo”  comprises  evaporation  q"ev 
bubble  agitation  q”ag*  It  is  observed  that  q”ev  became  larger 
in  the  inlet  region  where  the  superheat  is  large,  while  q”ag 
becomes  dominant  in  the  exit  region  where  the  small  wall 
superheat  reduces  the  contribution  of  q”ev* 

In  the  subcooled  region  near  the  test  section  outlet, 
bubble  agitation  is  significant  and  the  wall  superheat  is 
small.  To  explain  the  gravity  effects  on  subcooled  heat 
transfer,  the  heat  balance  in  the  boiling  test  section  was 
considered.  The  wall  superheat  for  the  subcooled  flow  boiling 
system  is  expressed  as 


(T  -T)  = 

V W S / 


q. 


1 1 


+hf^ 
ev  ag  rc 


+hf^ 

ag  ic  f ^ — 77  \ 

h +h  +h, 

ev  ag  fc 


(5-68) 


For  subcooled  flow  boiling  with  negligible  contribution  of 
qev”/  such  as  for  Case  F1035,  Eg.  (5-68)  is  further  reduced  to 


As  shown  in  Table  5-5,  hgg  decreases  significantly  in 
microgravity,  while  hf,,  increases  due  to  the  enlarged  void 
fraction.  As  a result,  the  first  term  on  the  right  hand  side 
of  Eg. (5-69)  is  increased.  The  second  term  on  the  right  hand 
side  of  Eg. (5-69)  remains  constant  unless  system  pressure  is 
changed.  The  above  discussion  suggests  an  increase  in  wall 
superheat  under  microgravity  conditions,  which  corroborates 
the  measurement  of  wall  superheats  shown  in  Figure  5-2 (f),  5- 
3 (f) , and  5-4 (f) . 

In  the  region  adjacent  to  the  ONB  point,  the  wall 
temperatures  are  higher  than  in  other  regions.  Although  the 
thickness  of  the  superheated  layer  is  thin,  the  number  of 
activated  nucleation  sites  is  considered  larger  than  in  other 
regions.  Under  microgravity  conditions,  the  wall  superheats 
decreased  as  predicted  by  the  ONB  model  previously  described. 
Although  the  decrease  in  pressure  as  well  as  saturation 
temperature  enhances  the  evaporation  heat  transfer,  the 
decrease  in  wall  superheat  suppresses  evaporation.  Since  the 
effect  of  wall  superheat  decrease  is  more  significant  than  the 
pressure  change,  the  resulting  heat  transfer  due  to 
evaporation  h^v  decreased. 

The  analysis  in  this  chapter  has  shown  that  gravity 
affects  the  flow  boiling  in  various  ways.  The  three  heat 
transfer  mechanisms  (evaporation,  bubble  agitation,  and  forced 
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convection)  were  discussed  in  view  of  their  dependence  on 
gravity.  Each  component  is  a function  of; 

h^v=-^i  (AT^atf  ^ '4  /Ys*  /«  fPr  ) (5-70) 

(5-71) 

4c=4  ,Pr  ,a  ,Xa)  (5-72) 

Among  these  parameters,  the  change  of  the  void  fraction 
affected  all  heat  transfer  mechanisms  of  flow  boiling  under 
variable  gravity  conditions.  The  detached  bubble  diameter  dg 
and  slip  velocity  Ug  were  found  to  be  significantly  affected 
by  gravity  change.  It  is  also  found  that  the  bubble  agitation 
heat  flux,  normally  included  in  the  boiling  heat  transfer  in 
the  existing  correlations,  plays  an  important  role.  However, 
Ub  and  dg  are  not  independent  of  each  other  and  a certain 
relationship  appears  to  exist  between  Ug  and  dg.  For  better 
understanding  of  physical  phenomena,  the  bubble  generation  and 
void  distribution  should  be  visually  observed  in  the  future 
experiments . 

The  bubble  detachment  mechanism  was  not  discussed  because 
the  point  of  the  net  vapor  generation  was  thermally  controlled 
under  the  conditions  used  in  the  present  study,  according  to 
Saha  and  Zuber's  criterion [ 54 ] . However,  gravity  effects  on 
the  bubble  detachment  will  be  significant  under  the  conditions 
where  the  detachment  criterion  is  hydrodynamically  controlled. 


CHAPTER  6 

SUMMARY  AND  RECOMMENDATIONS 
6.1  Discussions 

Existing  correlations.  Most  of  correlations  describing 
various  phases  of  flow  boiling  have  been  derived  under 
normal  gravity  conditions  assuming  a constant  value  for  the 
gravitational  acceleration.  The  exponents  of  Reynolds 
number,  Prandtl  number,  and  Grashof  number  obtained  through 
the  dimensional  analysis  are  expected  to  assume  different 
values  for  very  low  level  of  gravitational  acceleration.  It 
is  somehow  uncertain  that  the  same  dimensionless  groups  be 
involved  in  the  correlations  describing  flow  and  heat 
transfer  under  microgravity  conditions.  For  example,  the 
predictions  made  by  Rohsenow's  correlation  underestimated 
the  wall  superheats  measured  in  microgravity  as  shown  in 
Figure  5-16  (a)  - (c) . This  is  due  to  the  fact  that  the 
correlation  contains  the  force  balance  of  a detaching  bubble 
in  pool  boiling.  Numerous  observations  of  pool  boiling 
experiments  in  microgravity  show  that  bubble  detachment  is 
quite  different  from  that  in  normal  gravity.  In  addition, 
the  mechanism  of  bubble  detachment  in  flow  boiling  is  more 
complicated  than  in  pool  boiling.  According  to  Rohsenow's 
correlation,  the  predicted  heat  flux  is  proportional  to  g^^^. 
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However,  the  present  study  indicated  a weaker  dependence  on 
g in  the  range  of  microgravity  conditions.  In  the  absence 
of  gravity,  the  contribution  of  buoyancy-induced  heat 
transport  is  decreased  and  other  mechanisms  such  as  latent 
heat  transport  and  bulk  convection  of  the  liquid  phase 
become  dominant.  The  experimental  results  qualitatively 
reflect  the  above  tendency. 

Drop  tower  experiment.  Although  a fair  number  of 
experimental  investigations  of  adiabatic  two-phase  flow 
under  microgravity  conditions  have  been  conducted,  limited 
data  on  flow  boiling  with  heat  addition  is  available.  The 
drop  tower  experiment  of  the  present  study  was  designed  to 
observe  the  microgravity  flow  patterns  of  Freon  113  in  a 
vertical  channel  with  a constant  heat  flux.  The  use  of 
transparent  Cl2Sn  heating  elements  coated  inside  the  flow 
channel  provided  a uniform  heat  flux  distribution  along  the 
channel  and  simultaneously  allowed  taking  pictures  of  flow 
patterns  with  a 35mm  still  camera.  This  configuration 
allowed  for  an  unambiguous  observation  of  bubble  growth  on 
the  heating  surface  with  flow,  which  was  unavailable  in  the 
previous  adiabatic  two-phase  flow  experiments.  The 
difference  between  adiabatic  and  nonadiabatic  experiment  was 
quite  noticeable  in  the  bubbly  flow  region.  The  Cl2Sn 
elements  with  a thickness  of  50^lm  was  suitable  for  transient 
heat  transfer  experiments  within  a limited  period  of  time 
because  the  heat  capacity  of  the  heater  was  negligibly 
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small.  The  vertical  flow  was  chosen  since  the  transient 
time  from  normal  gravity  to  microgravity  in  the  vertical 
channel  is  shorter  than  the  horizontal  flow. 

Observed  flow  patterns.  A continuous  transition  from 
bubbly  to  annular  flow  patterns  was  also  observed  in  the 
vertical  flow  boiling  channel  under  microgravity  conditions. 
However,  the  detailed  structures  of  each  pattern  are 
different  from  those  in  normal  gravity. 

In  bubbly  flow,  two  types  of  bubbles  were  observed: 
nucleation  bubbles  on  the  heating  surface  and  coalesced 
bubbles  in  the  bulk  liquid  flow.  The  average  size  of 
nucleation  bubbles  under  growth  in  microgravity  was  larger 
than  those  in  normal  gravity,  indicating  that  the 
contribution  of  the  buoyancy  force  to  bubble  detachment  was 
nevertheless  significant  under  the  low  flow  conditions  of 
the  present  study.  It  was  also  noticed  that  the  number  of 
visible  nucleation  sites  was  increased  under  microgravity 
conditions.  This  is  due  to  an  increase  in  wall  temperature 
caused  by  reduced  bulk  convection.  The  main  reason  of  the 
reduction  in  bulk  convection  is  the  absence  of  buoyancy 
force,  which  otherwise  accelerates  bubbles  relative  to  the 
liquid  phase.  As  a result,  the  turbulence  generated  by  the 
bubbles  is  significantly  decreased  under  microgravity 
conditions. 

Coalescence  of  the  detached  bubbles  was  also 
significant  under  microgravity  conditions.  It  also  observed 
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that  the  bubble  coalescence  is  more  prominent  in 
microgravity  where  the  relative  velocity  between  bubbles  and 
the  liquid  phase  becomes  small.  This  assertion  is  supported 
by  the  observation  of  bubbly  flow  during  the  2 . 4g  period  in 
the  flight  experiment.  It  was  also  observed  that  the 
bubbles  with  smaller  diameter  moved  faster  in  2.4g  than 
those  in  normal  gravity.  The  size  of  coalesced  bubbles  is  a 
function  of  the  relative  velocity  between  the  phases,  which 
is  controlled  by  the  gravitational  acceleration. 

The  structure  of  annular  flow  in  microgravity  was  quite 
different  from  that  in  normal  gravity.  The  interface  was 
wavy  and  small  bubbles  were  trapped  in  the  thin  liquid 
films.  The  wavy  interface  was  formed  since  the  disturbance 
caused  by  the  vapor  motion  was  not  diminished  easily  in  the 
absence  of  gravity.  As  suggested  by  Klauser[25],  gravity 
affects  the  annular  flow  heat  transfer  through  the  liquid 
film  thickness.  The  formation  of  small  bubbles  seems  to  be 
due  to  a decrease  in  the  heat  transfer  rate  caused  by  the 
thickening  of  the  liquid  film  under  microgravity  conditions. 

In  the  slug  flow  region,  the  spherical  shape  at  the 
ends  of  the  vapor  slug  indicates  that  the  relative  velocity 
was  close  to  zero. 

Radial  distribution  of  the  vaoor  phase.  The  increase 
in  the  void  fraction  was  expected  by  the  theory  of  two-phase 
flow  since  the  slip  ratio  approaches  one  in  microgravity. 
However,  the  flow  pattern  observation  also  identified  a 
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significant  difference  in  the  radial  distribution  of  the 
vapor  phase  in  the  bubbly  flow  region  under  microgravity 
conditions.  This  is  due  to  the  existence  of  large  coalesced 
bubbles  moving  near  the  center  line.  After  detachment  from 
the  wall,  if  the  bubbles  move  faster  than  the  surrounding 
liquid  with  a velocity  gradient,  the  bubbles  are  driven 
toward  the  wall  by  the  lift  force.  The  peaks  of  the  vapor 
phase  distribution  near  the  wall  are  thus  formed  in  normal 
gravity.  However,  the  magnitude  of  the  lift  force  becomes 
smaller  in  microgravity  because  the  relative  velocity 
approaches  zero.  The  detached  bubbles  migrate  across  the 
cross  section  and  form  larger  bubbles.  As  a result,  a 
broader,  single  peak  is  formed  near  the  center  line  in 
microgravity . 

In  this  study,  a certain  dependence  between  the 
coalesced  bubble  diameter  and  the  gravitational  acceleration 
was  observed.  However,  a conclusive  explanation  of  the 
coalescing  process  of  multiple  bubbles  based  on  the  limited 
number  of  data  in  this  work  is  not  feasible. 

The  drift  flux  model  was  considered  and  the  values  of 
distribution  parameter  Cq  under  normal  and  microgravity 
conditions  were  calculated.  It  is  evident  that  the  value  of 
Cq  becomes  larger  than  one  in  the  bubbly  flow  region  under 
microgravity  conditions  because  of  the  vapor  phase 
concentration  near  the  center  line.  When  the  drift  flux 
model  is  used,  the  change  in  the  value  of  Cq  impacts  the 
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heat  transfer  calculation  which  is  dependent  on  the  average 
void  fraction.  It  was  found  that  the  gravity  change  affects 
the  void  fraction  through  both  the  radial  distribution  and 
the  relative  velocity  of  the  vapor  phase. 

Velocity  and  temperature  distributions.  One  of  the 
unique  features  of  Sato's  method  is  that  the  contribution  of 
bubble  agitation  to  the  total  heat  transfer  is  separately 
discussed.  Normally,  this  term  is  included  in  the  boiling 
heat  transfer  with  the  evaporation  heat  transfer.  Sato's 
approach  is  appropriate  for  the  present  analysis  since  it  is 
through  bubble  motion  that  gravity  affects  the  two-phase 
flow  heat  transfer.  The  calculated  velocity  profile  for 
0.02g  showed  an  increase  in  liquid  velocity  near  the  wall. 
The  calculated  wall  temperature  for  0.02g  was  higher  than 
that  for  l.Og,  indicating  a decrease  in  the  heat  transfer 
rate. 

Sato  proposed  that  the  eddy  diffusivity  due  to  bubble 
agitation  is  a product  of  the  void  fraction  a,  the  average 
bubble  diameter  dg,  and  the  relative  velocity  Ug.  It  was 
found  that  the  velocity  and  temperature  profiles  were 
sensitive  to  the  choice  of  the  above  variables.  Improper 
selection  of  the  value  for  Ug  lead  to  numerical  divergence. 
In  this  analysis,  the  measured  values  in  the  drop  tower 
experiment  were  used  as  input  data.  However,  the 
observation  of  bubbly  flow  under  variable  gravity  conditions 
suggested  that  a certain  relationship  exists  between  dg  and 
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Ug.  It  should  be  noted  that  the  values  of  these  parameters 
cannot  be  chosen  independently. 

Pressure  drop.  As  shown  in  Figure  4-15,  the 
accelerational  component  of  pressure  drop  was  increased  in 
microgravity  because  the  average  void  fraction  was 
increased.  Also,  the  frictional  pressure  drop  was  increased 
in  microgravity,  especially  in  the  low-cjuality  region.  This 
result  qualitatively  agrees  with  the  results  of  the  wall 
shear  stress  calculation  obtained  by  Sato's  method.  This  is 
because  the  liquid  phase  flowing  in  the  reduced  flow  area 
was  accelerated  and  the  velocity  gradient  became  larger  in 
microgravity.  The  increase  in  the  void  fraction  due  to  the 
reduced  gravity  affects  both  the  accelerational  and 
frictional  components  of  the  pressure  drop  in  the  two-phase 
flow. 

Subcooled  heat  transfer  experiment.  Observation  of  the 
flow  patterns  in  the  drop  tower  experiment  indicated  a 
significant  change  in  the  two-phase  flow  structures  in  the 
subcooled  region  under  microgravity  conditions.  The  second 
series  of  experiments,  the  reduced-gravity  airplane  tests, 
were  conducted  to  investigate  the  gravity  effects  on 
subcooled  heat  transfer.  A number  of  previously  performed 
two-phase  flow  boiling  experiments  using  pump-driven  closed 
loop  systems  had  transient  times  on  the  order  of  minutes 
upon  gravity  change.  These  transient  times  are  several 
times  longer  than  the  microgravity  period  available  on  the 
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microgravity  airplanes.  As  results  of  the  insufficient 
transient  times,  several  investigators  reported  conflicting 
results  on  microgravity  heat  transfer.  To  shorten  the 
transient  time,  the  present  system  employed  an  once-through, 
pressure-driven  piston  system  and  maintained  the  same 
conditions  at  the  inlet  of  the  boiling  test  section 
throughout  the  trajectories.  The  selection  of  vertical  flow 
also  helped  reduce  the  transient  time  because  the  flow 
pattern  was  already  axi-symmetric  in  normal  gravity. 

For  the  analysis  of  subcooled  heat  transfer,  the  flow 
rate,  the  void  fraction,  the  system  pressure,  bulk 
temperatures,  and  wall  temperatures  on  the  boiling  test 
section  were  measured  under  microgravity  conditions. 

The  measured  wall  temperature  profiles  indicated  that 
the  wall  superheat  at  the  onset  of  nucleate  boiling  was 
decreased,  while  that  in  the  subcooled  region  increased 
during  the  microgravity  trajectories.  The  analysis  of  the 
onset  of  nucleate  boiling  showed  that  the  former  change  was 
due  to  the  enhanced  bulk  convection  caused  by  an  increase  in 
the  flow  rate,  not  to  the  gravity  change.  However,  the 
latter  change  of  the  wall  superheat  in  the  subcooled  region 
was  induced  by  the  gravity  change.  It  was  observed  in  the 
drop  tower  experiment  that  bubble  motion  becomes  lethargic 
in  vertical  flow  under  microgravity  conditions.  Sato's 
calculation  showed  that  the  turbulence  generated  by  bubble 
motion  yields  a major  contribution  to  subcooled  heat 
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transfer.  It  is  concluded  that  the  reduced  gravity 
deteriorates  subcooled  heat  transfer  in  vertical  flow 
boiling  through  the  lethargic  motion  of  the  bubbles.  It  is 
interesting  to  note  that,  in  an  extreme  case  of  zero 
gravity,  the  velocity  and  temperature  profiles  in  vertical 
flow  approach  those  of  single  phase  flow  for  the  same  mass 
flux  because  the  turbulence  generation  due  to  bubbles  is  no 
longer  available. 

Gravity  effects  on  subcooled  boiling.  The  subcooled 
heat  transfer  mechanisms  consist  of  forced  convection  of  the 
liquid  phase,  bubble  agitation,  and  evaporation.  The 
contributions  of  each  heat  transfer  component  were 
calculated  using  Sekoguchi's  method.  The  bubble  agitation 
heat  flux  was  significantly  decreased  under  microgravity 
conditions.  On  the  other  hand,  the  forced  convective  heat 
flux  was  increased  due  to  the  increased  void  fraction. 

The  ratio  of  the  evaporation  heat  flux  to  the  bubble 
agitation  heat  flux  was  obtained  by  Bowring's  correlation. 
The  values  of  the  above  ratio  were  on  the  order  of  4.0. 

The  evaporation  heat  flux,  calculated  based  on  Bowring's 
correlation,  accounted  for  less  than  10%  of  the  total  heat 
flux  under  the  conditions  of  the  present  study.  This  ratio 
is  relatively  small  compared  to  that  in  saturated  boiling 
because  the  bubble  diameter  on  the  heating  surface  is 
smaller  in  the  subcooled  region.  According  to  the  analysis 
based  on  Sekoguchi's  method,  the  bubble  diameter  in 
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subcooled  flow  is  limited  by  the  thickness  of  the 
superheated  layer,  which  is  controlled  by  the  bulk 
convection. 

The  analysis  in  Section  5.4.3  showed  that  gravity 
affects  the  bubble  growth  rate  through  the  thickness  of 
superheated  layer  surrounding  the  lower  part  of  the  bubble. 
The  thickness  of  the  superheated  layer  estimated  by 
Sekoguchi's  method  was  slightly  increased  in  microgravity 
owing  to  the  deteriorated  bulk  convection. 

The  effect  of  gravity  on  the  bubble  detachment  was  not 
discussed  in  this  study.  The  flow  conditions  near  the 
bubble  detachment  point  were  classified  as  the  thermally 
controlled  region  defined  by  Saha  and  Zuber[54].  When  the 
inlet  subcoolings  are  much  smaller  than  those  of  the  present 
study,  the  bubble  detachment  will  be  controlled 
hydrodynamically . In  such  a case,  loss  of  the  buoyancy 
force  is  expected  to  have  a significant  impact  on  the 
evaporation  heat  flux. 


6.2  Conclusions 


a)  Experimental  rigs 

1.  Flow  patterns  and  bubble  growth  on  the  heating  surface 
were  observed  under  microgravity  conditions  in  a 15m 
drop  tower  facility.  Total  pressure  drops  along  the 
boiling  test  section  were  measured  and  the  results  were 
analyzed. 
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2.  The  second  series  of  flow  boiling  rig  was  designed, 
fabricated,  and  tested  aboard  the  NASA's  reduced- 
gravity  Learjet.  Wall  temperatures,  bulk  temperatures, 
and  pressure  drops  in  vertical  up-flow  were  measured. 
b)  Flow  analysis 

3.  Stable  boiling  in  the  subcooled  region  was  observed 
during  the  18  seconds  of  microgravity.  Flow  rates  were 
slightly  increased  due  to  the  static  pressure  balance 
of  the  vertical  flow  lines. 

4.  Bubbly,  Slug,  Annular  flow  patterns  were  observed  under 
microgravity.  Flow  characteristics  and  structures  of 
flow  regimes  were  significantly  different  from  those 
under  normal  gravity  conditions. 

5.  Under  similar  flow  conditions,  the  average  void 
fraction  in  vertical  up-flow  was  higher  due  to  the 
decreased  slip  velocity. 

6.  In  the  bubbly  flow  region,  radial  void  fraction  profile 
exhibited  two  peaks  near  the  heating  surface  under 
normal  gravity  conditions.  A broad,  single  peak  was 
observed  under  microgravity  conditions. 

7 . The  accelerational  and  frictional  components  of 
pressure  drop  increased  in  the  low-quality  region. 
Gravity  change  affects  the  pressure  drop  through  the 
flow  pattern  change  and  void  fraction. 
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c)  Heat  transfer  analysis 

8.  Gravity  effects  were  negligible  in  single-phase  flow 
and  at  the  onset  of  nucleate  boiling  under  the 
conditions  of  present  study.  The  decrease  in  wall 
temperature  at  the  onset  of  nucleate  boiling  was  caused 
by  the  decrease  in  system  pressure  and  the  increase  in 
flow  rate. 

9.  The  wall  superheat  in  the  subcooled  region  was 
increased  under  microgravity  conditions.  This  is  due 
to  the  decrease  in  bubble-induced  turbulence  in  the 
absence  of  gravity.  The  above  fact  indicates  an 
increase  in  evaporation  heat  flux. 

10.  Enlarged  bubble  size  and  the  reduced  slip  velocity 
oppositely  contributed  to  the  bubble- induced 
turbulence.  The  total  bubble-induced  turbulence 
diminished  as  the  gravity  level  became  smaller.  As  a 
result,  there  was  a reduction  in  the  heat  transfer  rate 
in  vertical  up-flow  under  microgravity  conditions. 

11.  The  influence  of  gravity  on  turbulent  flow  heat 
transfer  was  smaller  than  that  of  laminar  flow  due  to 
inherent  turbulence  in  turbulent  flow. 

12.  Total  two-phase  heat  flux  ratio  (hip/hfo)  of  subcooled 
boiling  in  vertical  up-flow  was  decreased  under 
microgravity  conditions. 
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6 . 3 Recommendations 

a)  Variable  heat  flux.  In  the  present  experiment,  the 
flow  rate  was  slightly  changed  while  using  the  same  value  of 
uniform  heat  flux.  However,  wall  superheats  for  various 
heat  fluxes  are  needed  to  develop  the  boiling  curve  for 
microgravity  applications.  In  particular,  the  data  near  the 
onset  of  nucleate  boiling  followed  by  the  nucleate  boiling 
region  should  be  gathered  and  analyzed  since  Rohsenow's 
correlation  fails  to  predict  the  boiling  heat  flux  in  this 
region. 

The  heat  flux  used  in  the  present  study  was  about  10 
times  smaller  than  other  studies.  Thus,  we  could  not 
compare  our  results  with  the  published  results  because  the 
number  of  activated  nucleation  sites  was  by  far  less  than 
others.  Besides,  most  correlations  for  the  two-phase  heat 
transfer  rate  or  void  fraction  tend  to  overestimate  our 
results.  To  perform  a reproducible  and  accurate 
microgravity  experiment,  the  first  task  should  be  directed 
toward  generating  the  data  which  can  be  analyzed  by  existing 
correlations  under  normal  gravity  conditions.  In  the 
Lear jet  cabin,  a maximum  of  80A  of  28V  dc  is  available  from 
the  on-board  power  supply.  The  current  of  14. 7A  was  used  to 
provide  a heat  flux  of  23.9kW/m^  into  the  0.5m  test  section 
in  the  present  study.  Since  electrical  power  available  on 
the  Lear jet  is  limited,  a maximum  of  approximately  240kW/m^ 
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will  be  the  upper  limit  of  heat  flux  for  the  present 
equipment . 

\ 

b)  Inlet  section  length  and  flow  rate.  Heat  flux  near 
the  inlet  region  is  normally  higher  than  that  of  the  fully 
developed  region.  As  shown  in  Sekoguchi's  method,  single- 
phase heat  transfer  prevails  upon  the  entire  region  of  the 
test  section.  Successful  prediction  of  the  boiling 
component  of  heat  flux  depends  on  an  accurate  measurement  of 
single-phase  heat  transfer.  The  inlet  section  of  the  test 
tube  should  be  long  enough  for  the  flow  to  fully  develop. 

In  the  present  study,  the  liquid  Reynolds  number  ranged 
from  4000  to  9000,  which  is  larger  than  the  laminar  flow 
criterion  of  2000  but  not  large  enough  for  fully  developed 
turbulent  flow.  In  this  intermediate  range,  the  effects  of 
laminar  flow  and  natural  convection  still  remain.  Higher 
flow  rates  should  be  used  since  most  existing  models  were 
developed  for  fully  developed  turbulent  flow.  The  pressure- 
driven  piston  system  used  in  this  investigation  provides  a 
maximum  mass  flux  of  500kg/m^/s  for  about  one  minute  and  the 
corresponding  Re  is  over  13,000.  In  such  a case,  the  heat 
flux  is  not  large  enough  to  achieve  saturation  conditions  at 
the  exit  of  the  test  section.  A new  heater  providing  higher 
heat  fluxes  needs  to  be  installed. 

We  experienced  a flow  interruption  during  the  2.4g 
period  of  flight  experiment  under  low  flow  rate  conditions. 
This  is  because  the  static  pressure  loss  exceeded  the 
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driving  pressure  in  the  piston.  To  prevent  flow 
interruption,  pressure  should  be  large  enough  compared  to 
the  static  pressure  loss  in  a vertical  tube.  This  can  be 
done  by  adjusting  the  pressure  regulator  and  replacing  check 
valves  on  the  pistons  preset  at  llpsig  for  the  present 
experiment . 

c)  Observation  of  vapor  phase  growth.  Void  fraction 
profiles  in  the  radial  and  axial  direction  convey  important 
information  for  analysis  of  flow  structure  and  pressure  drop 
mechanism  under  microgravity  conditions.  More  sophisticated 
models  for  pressure  drop  calculation  requires  the  growth  of 
the  void  fraction  in  the  axial  direction  must  be  accurately 
accounted  for.  In  addition,  bubble  growth  at  the  onset  of 
nucleate  boiling  and  bubble  size  at  the  moment  of  detachment 
needs  to  be  observed  to  understand  the  gravity  effects  on 
subcooled  boiling.  In  this  work,  the  void  fraction  was 
measured  at  the  exit  of  the  test  section.  To  observe  the 
evolution  of  the  vapor  phase  in  a boiling  channel,  Hino  and 
Ueda[41]  used  a transparent  plastic  annular  flow  channel 
with  a cylindrical  heater  at  the  center  (Figure  6-1) . Their 
system  gives  a clear  view  of  vapor  phase  evolution  all  along 
the  tube.  The  use  of  similar  geometry  is  recommended  for 
the  future  studies. 

d)  Contribution  of  bubble  agitation.  One  of  the  main 
mechanisms  influenced  by  gravity  is  the  heat  transfer  due 
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Figure  6-1  Transparent  annular  flow  channel 

for  observation  of  vapor  phase 
evolution . 
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to  bubble  agitation.  Sato  et  al.  described  the  eddy 
diffusivity  due  to  bubble  agitation  as  a function  of  the 
average  bubble  diameter  and  the  relative  velocity  of  vapor 
phase  to  the  liquid  phase.  The  average  bubble  diameter  was 
obtained  from  the  photographs  taken  in  our  experiment,  while 
the  relative  velocity  of  vapor  phase  was  still  unknown  under 
reduced  gravity  conditions.  The  drift  flux  model  was  used 
to  obtain  the  relative  velocity.  Under  an  ideal  case  of 
zero-gravity,  the  relative  velocity  disappears.  For  the 
practical  cases,  however,  the  residual  gravity  on  the  order 
of  lO.gg  always  exists  even  in  spacecraft  in  orbit.  In  the 
present  study,  it  is  expected  that  small  magnitudes  of 
relative  velocity  still  remain  under  the  gravity  levels  of 
0.02g  and  O.Olg  for  the  drop  tower  and  the  Learjet  flight 
experiments,  respectively.  Thus,  the  competitive  effects  of 
enlarged  bubbles  and  the  reduced  relative  velocity  on  the 
eddy  diffusivity  under  microgravity  conditions  needs  to  be 
further  examined. 

An  experimental  apparatus  is  proposed  to  investigate 
the  contribution  of  bubble  agitation  to  the  total  heat  flux. 
Figure  6-2  shows  a schematic  of  the  air-water  two-phase  flow 
system.  At  the  bottom  of  the  test  section,  air  is  supplied 
through  orifices  to  simulate  the  vapor  bubbles.  A specific 
heat  flux,  small  enough  not  to  initiate  boiling,  is  supplied 
from  outside  the  tube.  An  array  of  sensitive  thermocouples 
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Figure  6-2  A proposed  two-phase  flow  channel 

for  measurement  of  bubble-induced 
heat  transfer. 
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moun'ted  on  the  tube  surface  will  detect  the  wall  temperature 
change  caused  by  bubble  agitation.  In  the  above  system,  the 
heat  transfer  mechanism  involves  forced  convection  by  liquid 
and  turbulence  caused  by  bubble  agitation.  Since  nucleate 
boiling  is  excluded  here,  the  contribution  of  bubble-induced 
heat  transfer  will  be  measured  with  high  accuracy.  The 
influence  of  gravity  on  the  bubble-induced  heat  transfer 
will  be  understood  by  conducting  the  above  experiment  and 
observing  bubble  behavior  under  microgravity  conditions. 


APPENDIX  A 

STATIC  PRESSURE  CALCULATION  IN  THE  VERTICAL  FLOW  SYSTEM 

The  static  pressures  in  the  vertical  tubes  in  the  Lear jet 
experiment  changed  significantly  during  the  gravity 
transition.  This  change  caused  a decrease  in  system  pressure, 
which  would  not  occur  in  the  single-phase  flow  system.  During 
the  flow  boiling  experiment  with  low  flow  rates  in  high 
gravity  levels,  flow  interruption  occurred  because  of  the 
insufficient  driving  pressure  in  the  piston  system.  To 
prevent  the  flow  interruption  for  future  experiments,  the 
static  pressure  balance  was  calculated  for  the  present  flow 
boiling  system. 

Figure  A-1  shows  a schematic  diagram  of  the  flow  lines 
and  piston  system  used  in  the  Lear jet  experiment.  Assuming 
that  accelerational  and  frictional  pressure  drops  are 
negligible  compared  to  static  pressure  drop,  the  static 
pressures  at  points  PI  through  PIO  are  computed  as  follows; 

PlO  ~ Pq  “ (20”-ho)  Pf  ' g + Ploss,C2 

P9  = Pio  + 20»  Pf'  g 

Pg  = P9  - 23"  pj  g 

P,  = Pg  + 8.5"  pj  g 

Pg  = P,  - 17"  pj  g 

P5  = Pg  + 29"  p„  g + 5.5"  Pf  g 
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igure  A-1  Static  pressure  balance  in  a varying  gravity  field. 
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P4  = Ps  - 17.5”  Pf  g 
P3  = P4  + 7”  Pf  g 
^2  ~ P3  ” 15.5”  Pf  g + Ploss,V2 
P = P2  + (23”-h,)  Pf  g + P.oss.ci 

where  Pf,  Pf',  p„,  and  p„'  are  the  average  densities  of  liquid 
at  the  test  section  inlet,  condenser  exit,  inside  the  test 
section,  and  condenser  exit,  respectively.  Pioss,v2f  Pioss.ci/  arid 
Pio3s,c2  are  the  pressure  losses  due  to  the  throttle  valve,  the 
piston  head  in  Cylinder  1,  and  the  piston  head  in  Cylinder  2, 
respectively.  Assuming  further  Pf  « Pf'  and  p„  « p„',  the 
summation  of  the  above  equations  yields 

^drv^^io33+[31.5"pfa-(h,-ho)  Pf]g  (A-1) 

where  P^j-v  is  the  driving  pressure  supplied  from  a high- 
pressure  air  tank,  Pioss  is  the  total  pressure  loss  due  to  the 
throttle  valve  and  friction  of  cylinder  piston  heads,  and  a is 
the  void  fraction.  In  order  to  maintain  flow  in  the  system, 
Pdrv  must  exceed  the  pressure  loss  and  the  static  pressure  head 
described  in  the  right  hand  side  of  Eg. (A-1) . It  indicates 
that  the  static  pressure  head  increases  as  the  void  fraction 
increases.  The  value  of  (hf-ho)  is  at  most  10”.  It  is 
considered  that  the  total  of  the  r.h.s.  of  Eq. (A-1)  exceeded 
the  driving  pressure  under  the  conditions  of  Case  F1015,  where 
the  void  fraction  was  large.  In  addition,  the  acceleration 
level  g became  2 . 4 times  larger  than  that  of  normal  gravity 
during  the  pushover  maneuver  before  entering  the  microgravity 
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trajectories,  only  to  increase  the  value  of  the  right  hand 
side  of  Eq. (A-1) . 

To  prevent  the  flow  interruption,  it  is  suggested  to; 

a)  increase  the  driving  pressure  P<irv  large  enough  to 
overcome  the  pressure  loss,  or 

b)  confine  the  experimental  conditions  to  the  small 
void  fraction  range. 

In  the  present  experiment,  the  latter  approach  was  taken  since 
the  driving  pressure  was  limited  to  llpsi  by  pre-set  check 
valves  on  the  cylinders.  Future  experiments  will  reinforce 
the  cylinder  structure  to  sustain  higher  driving  pressures. 


APPENDIX  B 
ERROR  ANALYSIS 


The  data  from  the  experiment  is  meaningful  only  within 
the  limits  of  the  accuracy.  An  analysis  on  the  errors  was 
conducted  using  the  error  propagation  formula. 

The  basic  measured  quantities  and  their  respective 
uncertainty  intervals  are  summarized  in  Table  B-1.  The 
confidence  interval  of  the  derived  quantities  are  calculated 
in  the  following  section. 

a^  Exit  quality.  The  exit  quality  is  given  by 


Ah  +h 


in 


h 


fg 


(B-1) 


where  the  increase  in  the  enthalpy  Ah  is  expressed  as 
Q/ (G- p)  [kJ/kg] . The  heat  input  Q is  obtained  from  the 
measured  voltage  V and  current  I.  The  error  associated  with 
Q is  obtained  by  referring  the  individual  uncertainties  of  V 
and  I in  Table  B-1.  The  resulting  relative  error  in  Q is 

W 

-J-  = ±0.0345 

0 

The  error  involved  in  the  mass  flux  G stems  from  the 
uncertainty  in  the  turbine  flow  meter  and  the  resolution  of 
the  8-bit  data  acquisition  system.  The  combined  accuracy  of 
G becomes 
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Table  B-1  Values  and 

quantities. 

confidence  intervals  of  measrued 

Quantity 

Confidence  interval 

1 . Pressure 

5.50±0.05  psig 

2 . Mass  flux 

0.370±0.002  GPM 

3 . Temperature 

65±1  'C 

4 . Voltage 

28.0±0.25  V 

5 . Current 

15.0±0.5  A 

6.  Length  of  tube 

500±2  mm 

7.  Radius  of  tube 

12.810.1  mm 
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-^=±0.0341 

The  values  of  enthalpy  at  various  temperatures  were  read 
from  the  ASME  steam  table.  The  accuracy  of  these  values  are 
conservatively  assumed  ±0.02.  Using  the  values  obtained  in 
Case  F1035,  the  interval  of  uncertainty  for  the  exit  quality 
is 

IV 

-^  = ±0.133 
A 

It  is  noted  that  the  addition  and  subtraction  of  enthalpy 
cause  an  increase  in  uncertainty  although  those  of  individual 
values  are  small. 

b)  Pressure  measurement.  The  value  of  uncertainty  for 
the  high-accuracy  pressure  transducers  provided  by  the 
manufacturer  is  ±0.1%,  including  thermal  and  hysteresis 
effects.  The  signal  from  the  transducers  are  introduced  to 
the  data  acquisition  system  with  a resolution  of  2.44mV  in  the 
range  of  ±5V.  The  uncertainty  from  the  data  acquisition 
system  alone  is  ±0.0049.  The  total  error  of  the  pressure 
measurement  is 

IV 

-±=±0.0011 

This  error  is  to  be  applied  to  the  actual  readings  at  the 
inlet  and  outlet  of  the  test  section. 

The  pressures  between  the  inlet  and  outlet  Pj,  were 
estimated  by  the  following  equation 
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W 


(B-2) 


where  z is  the  distance  from  the  inlet,  L is  the  total  length. 
Pin  and  Pout  are  the  inlet  and  outlet  pressures,  respectively. 
Using  the  value  of  the  pressure  in  the  middle  of  the  test 
section,  the  uncertainty  of  P^  yields 


The  errors  involved  in  the  frictional  pressure  drop 
measurement  are  larger  than  those  in  the  pressure  measurement. 
The  major  sources  of  the  error  are  the  vapor  quality  and  the 
Reynolds  number.  The  uncertainty  of  the  vapor  quality  was 
discussed  in  the  previous  section.  The  uncertainty  of  the 
Reynolds  number  is  estimated  ±0.12  in  the  present  study.  The 
relatively  large  error  of  the  Reynolds  number  comes  from  the 
uncertainty  of  the  viscosity.  The  pressure  drop  in  the 
single-phase  region  is  calculated  by  Eq. (4-41)  and  Eq. (4-42) . 
The  estimated  uncertainties  in  this  region  are 


^ =±0.0016 

K. 


=±0.138 


For  the  two-phase  flow  region,  the  individual  errors  from 
accelerational  and  gravitational  components,  calculated  by 
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Eq. (4-36)  and  (4-37) , are  accumulated.  The  uncertainties  of 
each  component  in  the  two-phase  region  yield 


=±0.192 


=±0.118 


The  frictional  component  is  obtained  by  subtracting  the  above 
components  from  the  measured  total  pressure  drop.  Thus,  the 
resulting  total  uncertainty  for  the  two-phase  frictional 
pressure  drop  becomes 


=±0.340 


Significant  errors  tend  to  build  up  during  the  subtraction  of 
the  pressure  components. 

c)  Temperature  measurement.  The  voltage  signal  from  the 
OMEGA'S  E-type  thermocouple  is  converted  into  the  unit  of 
temperature  using  a 9th-order  polynomial  provided  by  the 
manufacturer.  The  nominal  accuracy  of  the  polynomial  in  the 
manual  is  ±0.5°C.  Since  the  relative  error  of  the  temperature 
measurement  is  larger  than  that  of  the  data  acquisition 
system,  the  latter  was  neglected  in  this  calculation. 

The  bulk  temperature  Tb  is  obtained  by  the  following  equation 

31  = (T„„. ^ +T„  (B-3) 

b ' out  in/  L ^ 

Using  the  error  propagation  formula,  the  error  in  the  bulk 
temperature  yields 
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=±0.025 


d)  Void  fraction  measurement.  The  diameters  of 
individual  babbles  were  automatically  measured  by  Particle 
Analysis  in  Macintosh's  Image  vl.31  by  choosing  an  appropriate 
scale.  The  error  of  this  visual  measurement  stems  from  the 
resolution  of  the  images.  The  scale  was  calibrated  so  that 
9.40  pixels  exist  in  1.0mm.  This  introduces  an  error  of 
±0. 106mm  in  the  measurement  of  bubbles.  The  measured  diameter 
of  the  bubbles  ranged 

0 . 19mm  ^ ^ 4 . 2mm 

As  is  noticed,  the  uncertainty  increases  as  the  bubble 
diameter  becomes  smaller.  The  resulting  accuracy  of  the  void 
fraction  measurement  ranges 

W 

±0.043  ^ ^ ±0.426 

a 

The  influence  of  the  error  occurred  on  measuring  the  diameters 
of  small  bubbles  is  considered  insignificant  since  the 
accumulated  volume  of  large  bubbles  are  dominant  in  the  void 
fraction  measurement. 


APPENDIX  C 

THERMAL  PROPERTIES  OF  FREON  113 

Variation  of  thermal  properties  of  Freon  113  was  not 
negligible  for  the  change  of  system  pressures  encountered  in 
the  Lear jet  experiment.  This  change  in  the  properties  affects 
thermodynamic  calculation,  bubble  growth,  and  heat  transfer. 
The  properties  covered  in  this  chapter  are  liquid  and  vapor 
densities,  latent  heat  of  vaporization,  specific  heat, 
viscosity,  thermal  conductivity,  thermal  diffusivity,  Prandtl 
number,  Jakob  number,  and  Ja*a°'^.  The  last  properties  are 
important  to  assess  the  influence  of  pressure  change  on  bubble 
growth . 
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Figure  B-1  Variation  of  liquid  density  as  a function  of 

pressure • 
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Figure  B-2  Variation  of  vapor  density  as  a function  of 

pressure . 


0.2 


Latent  heat,  (J/Kg) 
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Figure  B-3  Variation  of  latent  heat  as  a function  of 

pressure . 
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Figure  B-4  Variation  of  liquid  specific  heat  as  a 

function  of  pressure. 


0.2 


Liquid  viscosity,  (Pa  sec) 
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Figure  B-5  Variation  of  liquid  viscosity  as  a function  of 

pressure . 


Liquid  conductivity,  (W/m/K) 
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Figure  B-6  Variation  of  liquid  thermal  conductivity 

as  a function  of  pressure. 


Thermal  diffusivity,  (m'^2/sec) 
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Figure  B-7  Variation  of  thermal  diffusivity  as  a 

function  of  pressure. 


Prandtl  number 
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Figure  B-8  Variation  of  Prandtl  number  as  a function  of 

pressure . 
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Variation  of  Jakob  number  as  a function  of 
pressure . 
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Figure  B-10  Variation  of  Ja"*"a^0.5  as  a function  of 

pressure . 
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